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Abstract 



Quantitatively understanding hadronic physics from first principles requires numerical so- 
lutions to QCD, known as lattice QCD, which are necessarily performed with a finite spatial 
extent at finite lattice spacing. Current simulations employ quark masses that are heavier than 
those in nature, as lighter quarks arc computationally too costly, the physical systems simu- 
lated are not significantly larger than the lattice spacing nor are they significantly smaller than 
the lattice size. These approximations, or lattice artifacts, modify the observable quantities of 
interest. To make a rigorous connection between the physical world and current lattice QCD 
simulations, these lattice artifacts must be understood. The tool to systematically understand 
these effects is effective field theory. 

In this work, we extend and apply effective field theory techniques to systematically un- 
derstand a subset of these lattice artifacts in addition to the underlying physics of interest. 
Where possible, we compare to existing lattice QCD. In particular, wc extend the heavy baryon 
Lagrangian to the next order in partially quenched chiral perturbation theory and use it to 
compute the masses of the lightest spin- 1/2 and spin-3/2 baryons to next-to-next-to leading or- 
der. We then construct the twisted mass chiral Lagrangian for baryons and apply it to compute 
the lattice spacing corrections to the baryon masses simulated with twisted mass lattice QCD. 

We extend computations of the nuclcon electromagnetic structure to account for finite vol- 
ume effects, as these observables are particularly sensitive to the finite extent of the lattice. 
We resolve subtle peculiarities for lattice QCD simulations of polarizabilities and we show that 
using background field techniques, one can make predictions for the 4 spin-dependent nucleon 
polarizabilities, quantities which arc difficult to access experimentally. 

We then discuss the two-pion system in finite volume, determining the exponentially small 
volume corrections necessary for lattice determinations of the scattering parameters. We also 
determine the lattice spacing artifacts that arise for a mixed-action lattice simulation of the 
two-pion system with Ginsparg- Wilson valence quarks and staggered sea quarks. Wc show that 
the isospin 2 scattering length has a near continuum like behavior, differing from the chiral 
perturbation theory calculation by a computable difference. 
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Chapter 1 

LATTICE QCD AND EFFECTIVE FIELD THEORY 

Quantum chromodynamics (QCD) is the fundamental theory of the nuclear strong force, 
describing the interactions between quarks and gluons, which bind into the hadrons we observe 
in nature; the proton, the neutron, the pions etc. QCD is a non-Abelian gauge theory ^ based 
upon the SU{3) color group j2|Ej- The non-Abelian nature of the theory gives it the property 
of asymptotic freedom [3] , which is necessary to describe the observed Bjorken scaling |^ [7] in 
deep-inelastic scattering processes, in which at high momentum transfer the individual partons, 
or hadron constituents, act like nearly free, point-like objects ^8^. Conversely, at lower energy 
scales relevant to nuclear physics, ;U ~ 1 GeV, the coupling between the quarks and gluons 
becomes strong, as{fi) ~ C(l); and our well established perturbative treatment of a quantum 
field theory (QFT) breaks down. Currently, we still lack a closed- form solution to QCD, making 
a rigorous description of the nuclear physics governing the world around us directly from QCD 
difficult, to say the least. 

A conceptual breakthrough in our understanding of field theories .9' has lead to a lattice 
formulation of quantum field theory, in which spacetime is treated discretely instead of con- 
tinuously. This technique naturally lends itself to numerical solutions which can be applied 
to the non-perturbative regime of QCD, and in fact, lattice QCD (LQCD) is currently our 
only known solution to this non-perturbative physics. With this powerful technique, we can 
ask a number of previously unanswerable questions about the world around us. It is now well 
accepted that QCD gives rise to quark confinement, as evidenced by the lack of any direct 
observation of isolated quarks and also by the quark potential plots determined with lattice 
QCD. Lattice techniques have also allowed us to gain at least a qualitative understanding of 
the hadron spectrum we observe experimentally. However, we would also like to know if we 
can use QCD to explain the structure of hadrons? Does QCD contribute to the observed 
CP violations in the universe (in addition to the known CP violations of the electroweak the- 
ory [iniiniiiiiiaiiiiiaiiainiiiaiiaiiniiiiiiiiiisiiMiiis usi), or do we need a cp 

violating extension to the Standard Model? Can we understand the formation and structure of 
the nuclei, deuterium to uranium, from the fundamental parameters of the QCD Lagrangian? 
These are just a few examples of many interesting questions which are actively being pursued 
with LQCD. It has only been in these last few years that we have seen serious quantitative 
progress being made to answer these questions, which has largely been possible through the 
advances in available computing power. 

The QCD Lagrangian, along with the electroweak interactions, is responsible for all of 
hadronic physics, and thus in principle expressible as a function of only 4 fundamental QCD 
parameters, /(Agc/j, m„, m^, m^), along with the fundamental electroweak parameters. Here, 
^QCD is the dynamically generated QCD scale where roughly speaking the QCD coupling 
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becomes strong. The light quark masses, which dominate the infrared hadronic physics, are 
given by m„, rrid and rus (the other three quarks are heavy enough that they can be integrated 
out of the theory). Lattice QCD, in principle, is a tool we can use to numerically determine these 
functions, at least for systems with a small number of hadrons. LQCD simulations turn out 
to be very computationally demanding, requiring some of the world's fastest super computers. 
Even still, approximations must be made to perform simulations in a finite amount of time. 

All LQCD simulations are necessarily performed with a finite spatial extent (finite box) and 
at finite lattice spacing. Typically, the box sizes are not significantly larger than the physical 
system of interest, the lattice spacing is not significantly smaller than the physical system, and 
currently, the quark masses are larger than those of nature. These approximations, or lattice 
artifacts, modify the observable quantities of interest. Therefore, to make a rigorous connection 
between the physical world and current lattice QCD simulations, these lattice artifacts must 
be understood. The tool to systematically understand these effects is effective field theory 
(EFT). In this thesis, we contribute formal steps in the direction of understanding hadronic 
physics from lattice QCD. We accomplish this by developing, extending and applying various 
effective field theory techniques to hadronic observables. Much of this work has been published 
previously and we collect it here in three categories. 

Baryon Masses 

Using what is known as heavy baryon chiral perturbation theory (HBxPT) |271 128j , one can 
compute the masses of the lowest lying spin-^ and spin-| baryons, as functions of the light 
quark masses. In this way, one can fit the derived mass formulae to lattice computations of 
the baryon spectrum, and determine the mass of the nucleon as an explicit function of the 
quark masses. Why do we want to know the mass of the nucleon from LQCD? In addition 
to providing information about HBxPT, it turns out that knowing the mass of the nucleon as 
a function of the quark masses, in particular the strange quark mass, allows one to constrain 
the possibility of kaon or hyperon condensation in neutron stars |29l I30j . Additionally, the 
expansion parameter of HBxPT is not as small as that of purely mesonic chiral perturbation 
theory (xPT), and it is therefore important to push the determination of the baryon masses 
to higher orders so one can test the convergence of HBxPT. To this end, we extend the heavy 
baryon Lagrangian to next-to-next-to leading order (NNLO) in HBxPT as well as extensions 
of this theory for LQCD, and compute the quark mass dependence of the spin-^ and spin-| 
baryons to 0{m^). This work has been published; 

• Andre Walker-Loud, Octet baryon masses in partially quenched chiral perturbation theory, 
Nucl. Phys. A 747 (2005) 476. 

• Brian C. Tiburzi and Andre Walker-Loud, Decuplet baryon masses in partially quenched 
chiral perturbation theory, Nucl. Phys. A 748 (2005) 513. 

• Brian C. Tiburzi and Andre Walker-Loud, Strong isospin breaking of the nucleon and 
delta masses on the lattice, Nucl. Phys. A 764 (2006) 274. 
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One exciting new lattice discretization technique is known as twisted mass LQCD 2>1 . We 
construct the twisted mass heavy baryon Lagrangian which incorporates the lattice artifacts 
arising from the twisted mass lattice action. We then apply this Lagrangian to compute the 
effects on the nucleon and delta spectrum arising from the twisting. We also discuss the subtle 
but very interesting consequences of performing twisted mass LQCD simulations which incor- 
porate both isospin violation in the up and down quark masses as well as the twisted mass 
effects. This work has been published; 

• Andre Walker-Loud and Jackson M. S. Wu, Nucleon and delta masses in twisted mass 
chiral perturbation theory, Phys. Rev. D 72, 014506 (2005). 

Electromagnetic Structure of the Nucleon 

The structure of the nucleon is interesting in its own right, from the nucleon parton •8_ dis- 
tributions, to the spin content of the nucleon |32| I33j to the electromagnetic moments and 
polarizabilities of the nucleon. Each of these is a measure of the content and structure of the 
nucleon, as well as a necessary building block to understanding the structure of nuclei. In this 
work we show how to determine the polarizabilties of the nucleon from lattice QCD calculations 
using background field techniques. We show that the polarizabilities are particularly sensitive 
to the finite size of the box used in the simulation. We address other subtleties which arise 
for LQCD determinations of the nucleon polarizabilities and additionally, we show how LQCD 
in conjunction with effective field theory can be used to predict the spin polarizabilities of the 
nucleon |34l 135115^ . which are difficult to access experimentally (in fact, there are 4 spin polar- 
izabilities, and currently only two linear combinations of them are measured). This work has 
been published; 

• William Detmold, Brian C. Tiburzi and Andre Walker-Loud, Electromagnetic and spin 
polarisabilities in lattice QCD, to be published in Phys. Rev. D (2006). 

Two-hadron interactions 

To understand hadronic physics from LQCD, in particular nuclear physics, one must understand 
multi-hadron systems with LQCD. A first step in this direction is the two-hadron system and in 
this work we focus on the two-pion system, which is both formally and numerically easier than 
the two-nucleon system. We study both the volume dependence as well as the lattice spacing 
dependence of the two-pion system. An elegant method has been developed for extracting 
two-particle infinite volume scattering parameters from finite volume energy levels of the two- 
particle system jS3 |SH1 [HHl HOI ISl ^21 HSl ■ We compute the exponentially suppressed volume 
corrections to this relation for the two-pion system and show that these effects will become 
important in the next generation of lighter pion masses simulated with lattice QCD. We also 
compute the lattice spacing and partial quenching corrections to the isospin 2 (7 = 2), two- 
pion scattering length, in particular for a mixed-action simulation employing Ginsparg- Wilson 
valence quarks jlS] and staggered sea quarks |46j . We demonstrate that a very nice cancellation 
of lattice spacing and partial quenching effects occurs if one expresses the scattering length 
in the effective field theory in terms of correlation functions one measures in the lattice QCD 
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simulations instead of in terms of the bare parameters of the theory. These results explain 
the remarkable success of the recent determination of the 1 = 2 tttt scattering length from a 
dynamical mixed action lattice QCD simulation j47j . 

Finally, based upon the insights gained from determining the lattice spacing and partial 
quenching effects mentioned above, we demonstrate the utility of mixed-action simulations 
employing Ginsparg- Wilson valence quarks. In particular, we show that through the one- loop 
order in the effective theory, the lattice spacing effects which arise can all be absorbed by 
multiplicative redefinitions of the coefficients of the continuum theory. We also demonstrate 
how at the next order in the chiral expansion, this renormalization scheme breaks down, and in 
particular, for staggered sea quarks, we show how the taste-breaking enters these mixed action 
simulations. Most of this work has been previously published, however the final chapter of this 
thesis is work in progress; 

• Paulo F. Bedaque, Ikuro Sato and Andre Walker-Loud, Finite volume corrections to tttt 
scattering, Phys. Rev. D 73, 074501 (2006). 

• Jiunn-Wei Chen, Donal O'Connell, Ruth S. Van de Water and Andre Walker-Loud, 
Ginsparg- Wilson pious scattering in a sea of staggered quarks, Phys. Rev. D 73, 074510 
(2006). 

• Jiunn-Wei Chen, Donal O'Connell and Andre Walker-Loud, Mixed-action simulations with 
Ginsparg-Wilson valence quarks, to be published. 



1.1 Lattice QCD 

Here we do not pretend to give a thorough introduction or review of lattice QCD. We aim 
only to review some key concepts and introduce common jargon useful to understanding this 
work. There are numerous reviews and books describing lattice QCD (and lattice QFT) in the 
literature. A few references the author finds particularly useful are j48| I49j . 

As we mentioned in the previous section, lattice QCD (or more generally lattice field theory), 
is a powerful numerical technique which allows us to study the non-perturbative regime of field 
theories, in particular, low-energy QCD. Lattice theories are formulated with the path integral 
formulation of quantum field theories. The generating function for QCD is 

Z = j DA^DqDqeicp{iSG[A^\+i j d^xq{ilp-mg)q^ , (1.1) 

where g is a quark field, rUg is the quark mass matrix, Af^ are the gluons, S'f^^] is the standard 
Yang- Mills actions, qlpq is the gauge-interaction between the quarks and gluons, and the func- 
tional integral is performed over the infinite set of possible field configurations of the quarks 
and gluons. The idea is then to discretize spacetime, which regularizes the theory by providing 
a natural ultraviolet (UV) cut-off set by the inverse lattice spacing, a~^. Additionally, one 
provides a finite spatial extent to the universe, with some user specified boundary conditions 
(most commonly periodic) , and this then provides a natural infrared (IR) cut-off of the theory, 
set by the inverse box size, L~^. In this way, the generating function is now a finite dimensional 
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integral and one can think about using computers to evaluate it. Unfortunately, this is still not 
possible as there are an infinite set of field configurations, and in Minkowski space, the integral 
is performed over the phase of an exponential, and thus every configuration is important, with 
large cancellations happening amongst configurations which are far from the classical action. 

To get around this problem, we can Wick rotate to Euclidean space, t — > —it, such that the 
functional integral becomes an integral over a real valued function, which is now weighted by 
exponentially damped field configurations away from the field configurations which give rise to 
the minimum of the action, 

2e = j DA^DqDqexp(^-SG[A^]- j d''xq{Ip+m,)(^ . (1.2) 

One can then apply Monte Carlo techniques to determine a set of field configurations which give 
the dominant contributions to Eq. (|1.2() , and in this way limit the number of field configurations 
needed to estimate physical quantities. To do this, one first does the integral over the quark 
degrees of freedom and uses the resulting fermionic determinant as part of the weighting in the 
Monte Carlo updating routines. In this way, the expectation value of some physical observable 
is given by 



N, 



configs 



{o) = —^ om, (1.3) 

configs ^ 

where 0[Ui] is a functional of the latticized gluon gauge fields, C/j, and the selection of these 
gauge fields is weighted by 

Det[^+mq] expi-Scp]) . 

In the limit of an infinite set of gauge configurations. Neon figs oo, this Monte Carlo averaging 
converges to the exact answer. To make these ideas more concrete, and to introduce some lattice 
jargon, we shall discuss these ideas in more detail with a specific example, the pion two-point 
correlation function. 

The first thing to do is construct and interpolating field for the pion which has an overlap 
with the physical pion in the lattice simulation. This interpolating field will excite all the 
states of the theory which have the same quantum numbers as the pion, but as the pion is the 
lightest mode of the theory, it will dominate the correlator. Let us consider the vr"*", for which 
an interpolating field is given by 

7r+(f,r) = d(f,r)75n(f,T) , (1.4) 

where d and u are the down and up quarks respectively, and 75 is used as the pions are pseudo- 
scalars. In Figure 11.11 we provide a cartoon of the pion two-point function, for which the 
correlator is given by 

( TT-ixE) vr+(0) ) = ^ j DA Det + m,] tr(g«(0, x^j) 75 ^^(xij, 0) 75 )exp (-5g[A]) 



N. 



configs 



-^configs 



tr(gu{0,XE)75GdixE,0)J5)[U^], (1.5) 
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Figure 1.1: The pion two-point correlation function. In this figure the curly lines represent 
gluons, and the grey blobs are sources/sinks for the tt"'". The outer solid lines (black online), 
which are connected to the source and sink are valence quarks. The closed inner grey lines (red 
online) represent dynamical sea quark - antiquark pairs. The fermionic determinant in Eq. 1)1.5(1 
involves the dynamical sea quarks, while the operators which make up the sources and sinks 
involve valence quarks. 



where the up and down quark propagators are given by 

g{x,y) = []p, + m,]-; . (1.6) 

The labels, 'u' and 's' stand for valence and sea respectively. The valence quarks are those which 
are tied to the source and sink, and thus their masses appear in the propagators, Eq. ()1.6|1 . 
arising from the sink and source operators. The sea quarks are the dynamical quarks which 
appear in quantum loop fluctuations and thus their masses appear in the fermion determinant 
of Eq. (|1.5() . The most expensive part of lattice QCD simulations is computing the fermion 
determinant, as it is a non-local operator, for which the computation time scales roughly as 
t ~ This is in contrast to the computation time of the quark propagators, which scale 

roughly as t ~ m~^. 

The process of generating of gauge configurations, in which the sea quarks are simulated, 
is independent of the process of constructing interpolating fields and computing correlation 
functions, in which the valence quarks are simulated. There is therefore no technical reason that 
masses of the valence and sea quarks must be the same. Moreover, the discretization methods 
(or Dirac operators) used in the valence and sea sectors can also be different. However, there 
are unphysical consequences for using different masses or Dirac operators, the most notable 
being a loss of unitarity. This is simple to understand as with either different valence and sea 
masses or Dirac operators, the internal quarks of some correlation function (the sea quarks) are 
not the same as the external quarks (the valence quarks) as is required by unitarity and the 
optical theorem. 

Before discussing the motivations of using these unitarity violating theories, we will first 
define some common lattice jargon, which for convenience we list in Table ITTTl When the sea 
and valence Dirac operators are not equal, this is known as mixed action (MA) or hybrid lattice 
QCD. Mixed action lattice simulations can never be tuned exactly to the QCD point as in any 
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Table 1.1: Common lattice QCD jargon. 





■nis / "Id 

TTls = 


Mixed Action (MA) or hybrid lattice QCD 




rUs ^ ruy 


partially quenched QCD (PQQCD) 




nig = rriy 


QCD 


NA 


TJiy / , nis = oo 


quenched QCD (QQCD) 



simulation, there will be different lattice spacing corrections to the valence and sea quarks which 
will result in different quark masses and dispersion relations even if the mass parameters are 
tuned equal. However, for sufficiently small lattice spacings, one can incorporate these lattice 
artifacts into the effective field theory (EFT) description of the mixed action LQCD simulation,^ 
and use these expressions to remove the lattice spacing effects, allowing an extrapolation to the 
QCD point. We shall discuss a particular example in some detail in Chapter 1^1 

A simulation with the same Dirac operator but different quark masses in the sea and valence 
sectors is known as partially quenched QCD (PQQCD). This theory is also unitarity violating, 
but has the advantage over MA simulations that at the level of the simulations, one can always 
set the valence and sea quark masses equal, and thus perform numerical simulations without 
any sicknesses. The important point here is that at the level of simulations, one knows that 
there is a limit, = m„ in which one is simulating QCD, the unitarity conserving theory. 
This also implies that the low energy effective theory of PQQCD must also contain this limit, 
for the same tuning of valence and sea quark masses. One can then perform simulations away 
from this limit, match the EFT description of observable quantities to the lattice correlation 
functions, and later by hand take the QCD limit, thus recovering the physics of interest.^ 

Lastly, one can set the sea quark masses to infinity, such that they can be integrated out 
of the theory. This is known as quenched QCD (QQCD). The motivation of this theory is 
that up until a few years ago, it was numerically too expensive to simulate dynamical quarks 
in a reasonable amount of time. It has also been found that QQCD typically introduces an 
uncontrolled error of < 20%. Thus, QQCD provides an inexpensive check of methodology 
and can provide approximate predictions for physical observables. We stress that there is no 
rigorous connection of QQCD to QCD, meaning there is no systematic way to correct the 
quenched approximation and include the effects of the neglected sea quarks. 

The motivation to consider these partially quenched simulations is related to the extreme 



^Of course, for arbitrarily small lattice spacings, these differences arising from the different lattice actions 
should become negligible, which practically means smaller than the statistical uncertainty of a given quantity. 
This would also imply that the need for a mixed action effective field theory description would not be necessary, 
but today and for the foreseeable future, these differences must be incorporated into the MA effective theory 
relevant for a given MA lattice simulation. 

■^The belief of this author is that this same methodology works for the MA theories, with the combined 
continuum limit (in which the Dirac operators become equivalent), despite not being able to tune to the QCD 
point in the simulations. This can only be verified or disproven by a matching of the MA effective field theory 
to MA LQCD simulations, and then performing the continuum, QCD limits. 
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cost of full scale lattice QCD computations (which are presently measured in tera-flop years). 
The first reason is that as stated above, it is numerically much cheaper to simulate valence 
quarks as compared to sea quarks. If one can make more use of the expensive gauge configura- 
tions containing the sea quarks, by using a number of different valence quark masses for each 
configuration, then one obviously gains much in terms of computing resources. Additionally, it 
is important to reduce the masses of the quarks used in the simulations, as it has only been in 
these last few years that the quarks masses have been in the chiral regime [511 152j such that 
one can make use of effective field theory techniques to compare the lattice QCD simulations 
to the physical world. By reducing the valence quark masses more than the sea quark masses, 
one can further probe the chiral regime, increasing the region of overlap between the parameter 
space where lattice simulations are performed and where one has confidence in the effective field 
theories (53j . 

The motivation to consider mixed action simulations, which are generalizations of partially 
quenched simulations, is also related to the numerical cost of full scale lattice QCD simulations. 
Additionally, the MIMD Lattice Computation (MILC) collaboration has made there lattice 
source code and gauge configurations publicly available [51]. For reasons which we will discuss 
in some detail in Chapters 101 and it is beneficial to use a different discretization technique for 
the valence quarks, than the technique used to generate the MILC configurations. Two groups 
in particular have determined some very impressive physical results using these mixed action 
techniques [Hi EHl EH EH EHl EHl EOI • 

There are good reasons to believe that the low energy effective theories of these unitarity 
violating lattice theories fully capture the unitarity violations introduced by either different 
masses for the sea and valence quarks, and also for the use of different discretization techniques 
in the sea and valence sectors. Provided this is true,^then the EFTs for these lattice simulations 
provide a rigorous means to describe the unitarity violations in the region where the quark 
mass and lattice spacing corrections can all be treated perturbatively, and thus remove their 
contribution to the lattice correlation functions and extract the physics of interest. We will 
come back to this issue below in Section 11.2.1^1 as well as Chapters IHl and d We now move on 
to describe some of the effective field theories necessary to understand lattice QCD simulations 
today, and extrapolate these LQCD computations to the physical world. 

1.2 Effective Field Theory 

What is an effective field theory (EFT)?^ Effective field theories are based upon simple ideas 
with powerful consequences. The underlying idea is that physical phenomena at low-energies 
(long wave-lengths) are not sensitive to the details of the high-energy (short wave-length) struc- 
ture of particles or fields. We use these ideas everyday in physics. In freshman undergraduate 
physics courses we teach that the change in gravitational potential energy of a particle of mass 

•^We stress again, that this can only be tested by performing MA simulations, and comparing them to the 
EFT predicted form of correlation functions. 

■^As with our introduction to LQCD, we do not pretend to give a comprehensive review of effective field 
theories. We merely summarize the main ideas and introduce the EFTs necessary to understand the work in 
this thesis. The following is a list of EFT references the author has found invaluable in the course of pursuing 
his degree in physics, [CTIIS^ 0151 1551 15^ . 
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m raised a height h above the surface of the Earth is given by 

AU = mgh, (1.7) 

where g is the acceleration of a massive particle induced by the Earth's gravitational field. This 
is actually the leading term in an effective expansion of the height of the object over the radius 
of the Earth, R. The actual change in potential energy has the form 

\ TT r ^ 

l\U = mg a 



mgh[l-'^ + o(^A] . (1. 



R + h 

R ' ^ \W 

This is a simple example of a classical effective theory (no quantum loop effects), of which 
there are many. Classical Newtonian mechanics is an effective theory of quantum mechanics for 
distances large compared to quantum length scales, set by h. The pattern we see from these two 
examples is that there is a separation of scales, in these cases length scales, which sets a range of 
validity of the theory. For heights comparable to the radius of the earth, Eq. (|1.8)) is no longer 
a good approximation of the true change in potential, as the neglected terms, Qi/R)"^ ~ 1 
are no longer small, and must be included. However, for /i << i?, Eq. (|1.8|1 is a very good 
approximation to the change in potential, and moreover, the neglected terms give us a rigorous 
estimate upon the theoretical uncertainty we have introduced by not including them.^ Let us 
define a small scale that serves as an expansion parameter for the effective theory, 

eG^\. (1.9) 

We see that e both acts as a small parameter, allowing computations to arbitrary precision 
(powers of e), and also is an indicator of the range of validity of the theory. These are charac- 
teristic features of effective theories, which also apply to quantum field theories. 

When we apply these ideas to quantum field theories, there are a few more steps to the 
process. The first task is to identify the appropriate degrees of freedom relevant to the physics 
questions we are interested in asking, i.e. what are the particles one would like to describe? 
Next, we must identify all the symmetries (and approximate symmetries) the underlying, or 
fundamental theory exhibits. We then construct a Lagrangian which includes all linearly inde- 
pendent operators consistent with the underlying symmetries (this is often an infinite number 
of operators). The power of the EFT treatment is the separation of scales, or the expansion 
parameters, such as Eq. (|1.9|) in the gravitation potential example. Often, there will be more 
than one ratio of relevant scales, let us call them £a and £5. All the operators in the effective 
Lagrangian will have associated powers of these expansion parameters. One then decides to 
which precision in 6^ and £^ we wish to know the physics, and computes all processes which 
contribute to this order of precision. For instance, in our gravitational potential energy exam- 
ple, the expansion parameter for heights on the order of 10 meters is eg ~ 1.6 x 10~^, and thus 
the first term in the series is a very good approximation. The (often infinite) set of operators 



^Actually, in this example, it is an exact uncertainty as the coefficient of the neglected term is exactly known. 
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in the effective Lagrangian we need are thus truncated by this desired precision, and we have a 
theoretical estimate on the neglected physics, the first order in and Sb not included. 

If we know the underlying theory, as in the case of the gravitational potential example, then 
the coefficients of the operators in the effective theory can be exactly determined by matching 
the effective theory to the underlying theory at an energy scale both theories converge in, and 
requiring the physics to be the same. If, however the underlying theory is unknown, as is the 
case with extensions to the Standard Model, or if the underlying theory is not computable at 
energy scales of interest, as with QCD, then there is no way a priori to know the coefficients 
of the operators in the effective theory, and they must be determined by fitting expressions 
of physical observables computed with the EFT either to experimental data or in the case of 
QCD, to lattice QCD simulations. The predictive power of an EFT is that these a priori unde- 
termined coefficients of the effective operators are universal. Once they have been determined 
by comparing to one set of experiments, they can be used to make predictions about other 
experimental quantities not used to constrain them. 

This is a crucial component of EFTs. As one has included all operators consistent with 
the symmetries in the effective Lagrangian, the EFT provides a model independent description 
of the physics in the region of applicability of the theory. The quantum loop contributions 
arising from the effective operators encode the long distance, or infrared (IR) physics, while the 
coefficients of the higher dimensional operators in the theory, known as low-energy-constants 
(LECs), encode the short distance, or ultra-violet (UV) physics. The IR physics is dominated 
by the lightest degrees of freedom in the theory, and at a given order are generally known, while 
the UV physics at a given order in and is dominated by the highest dimensional operators 
included in the computation. The LECs of these operators are generally unknown and must be 
determined to have predictive power. We now move on to introduce the low energy effective 
theory of QCD, chiral perturbation theory. 



1.2.1 Chiral Perturbation Theory 



The most well known effective field theory is the Fermi theory of weak interactions [lOl lllj . 
Perhaps the next most well known EFT is chiral perturbation theory (xPT) which was intro- 
duced by Weinberg [HZllEHllEni and systematically developed by Gasser and Leutwyler [701 [7T|. 
One can not derive from QCD as that would require an exact closed-form solution of the 
QCD equations of motion. However, with one phenomenologically based assumption, one can 
construct the most general Lagrangian consistent with the symmetries of the QCD Lagrangian. 
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The quark content of the QCD Lagrangian is^ 

.9i 



^QCD = ^ Qi ilp-mQ 

i,j=u,d,s 

= X] [^^'^ ^^^i ^i'^ + ^"-'^ ^^^i ^^'^ + ^^'^ "^<3'j ^J'-f^ + ^^'^ "^Q'j ^J'-^ 
where we use the quark mass matrix 



(1.10) 



irid 



(1.11) 



and the left and right handed quark fields are defined as 



1 



75 



1 



QR 



■75 



(1.12) 



We can see from Eq. (|l.lUj) that the kinetic term of the QCD Lagrangian is invariant under 
independent L and R chiral transformations, where L,R € SU{3)L^fi. The quark mass term 
explicitly breaks this chiral symmetry of SU{3)l (gi5C/(3)/j down to the vector subgroup. How- 
ever, for sufficiently small quark masses, as compared to Aqcd, this explicit breaking can be 
treated perturbatively. In the zero quark mass limit, the QCD Lagrangian exhibits an exact 
chiral symmetry. However, in nature we observe that the parity partner of the nucleon is sig- 
nificantly heavier than the nucleon. We also observe that all the observed hadrons composed of 
the three light quarks, H, have masses mn ~ 0(1 GeV), except for a handful of pseudo scalars, 
the pions (m^ ~ 135 MeV), the kaons {niK ^ 498 MeV) and the eta (M^ ~ 548 MeV). This 
phenomenological evidence leads us to postulate that the vacuum of QCD spontaneously breaks 
the global chiral symmetry of the QCD Lagrangian, Eq. (|1.1U|) down to the vector subgroup (83) . 



SU{Nf)L(^SU{Nf)R^ SU{Nf)v, (1.13) 

creating the 8 pseudo-Nambu-Goldstone modes jMJ |SS1 EH] which encode the long distance 
fluctuations about the many possible minima of the QCD vacuum. In the zero quark mass limit, 
these 8 mesons would be exact Nambu-Goldstone modes, and thus massless. The inclusion of 
the quark mass parameters, however, gives rise to the non-vanishing mass of these mesons. 

In the following, we will consider the theory with two light quarks, the up and down, as 
well as the theory with three quarks, the up, down and strange. The mass of the strange 
quark is significantly heavier than the up and down quarks but still lighter than Aqcd- The 
inclusion of the strange quark in the theory is thus necessary, but the convergence of the chiral 
expansion is not as good. This is especially true for the theory including heavy baryons, as 

We should also include the theta term in the Lagrangian, d-^^G^vC^" , as this term can be rotated into a 
phase of the quark masses by a chiral transformation. This theta term gives rise to CP violations in QCD and 
also to the neutron electric dipole moment (EDM), which has yet to be found experimentally |72|. One can 
also theoretically understand the neutron EDM using effective field theory methods Ellll EH Esl E] and 
lattice QCD 78 79 80 8f 82 . However, in this thesis, we do not explore any CP violating processes, and 
thus shall not include this term. 
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the expansion parameter of heavy baryon (HBxPT) scales with the meson mass, and 

not the meson mass squared as in the purely mesonic theory. Therefore, sometimes it will be 
prudent to explicitly include the strange quark, and other times it will be more beneficial to 
include its effects indirectly through modifications of the LECs. This can only be determined 
phenomenologically through a comparison of the two theories with their respective LECs fit 
either to experiment or lattice QCD simulations. We will show this explicitly in the next two 
sections in which we first focus on the two-flavor theory, and then the three-flavor theory. 



Two Flavors: The Pions 

The hypothesis is that quarks condense in the QCD vacuum, 

(0|gi,gi|0) = A5*^ (1.14) 

where i,j are flavor indices and A is a dimension-3 quantity whose scale is set by KqcD) 
A ~ [A-qcdY- Under a chiral transformation, 

(0| in ql |0) (0| tn RV L\ q\ |0) = A {LR^T 

= AS*^ (1.15) 

and we see that this quark condensate is only invariant under the vector subgroup, L = R = V ^ 
for which Y^^ = 5*-' . Otherwise, for L ^ R, S*-' corresponds to a different vacuum of QCD than 
that of Eq. (|1.14|1 . and in the massless quark limit, these different vacua would be degenerate. 
The spontaneous breaking of the chiral group then would give rise to exact Nambu-Goldstone 
bosons representing the long wavelength fluctuations of the QCD vacuum, S. In the case of 
SU{2), for example, three generators are broken by the vacuum, giving rise to 3 pions. We then 
wish to parameterize these bosons with a linear combination of the broken generators |87| I88j . 
The choice is not unique, (although all forms must give the same physics) and here we choose 
the non- linear sigma model form. We replace S with a non- linear realization of the pions, 

S = exp (^^) with ^=(^ , (1.16) 

where we use the normalization / ~ 132 MeV. From Eq. (jLlSf) . one can see that the S-fleld 
transforms under the chiral group as 

E — >LT,R^. (1.17) 

To construct the most general chiral Lagrangian, we must first include the quark mass param- 
eters by utilizing spurion analysis, in which one first treats the quark mass matrix as a spurion 
field which transforms under the chiral group as 



LmqRK (1.18) 
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Setting mq = diag(mti, m^),'' explicitly breaks the chiral symmetry down to the vector sub- 
group, SU{2)y. The most general Lagrangian consistent with chiral symmetry at leading order 
is then ^ ^ 

£ = ^tr (a^SS'^S^) + ^^^tr (Sm{^ + mgSt) , (1.19) 
where the LEG Bq is defined by 



So = Hm^^^o . (1.20) 



The normalization of the first term in Eq. 1)1.19(1 is fixed by the canonical normalization for a 
scalar kinetic operator, which is seen by expanding the S-field in powers of (j)/ f ■ To quadratic 
order in cf), Eq. (|1.19|) gives, ^ 

C = l^vr+p + i(97r0)2 + i?o(m„ + m,) ^+7r- + iyrVO j + . . . (1.21) 

and thus the pion masses are given at leading order by 

ml = BQ{mu + md). (1.22) 

The theory gets much more interesting when one includes the quantum loops arising from this 
Lagrangian. To be consistent with EFT techniques, we then need to include tree-level terms 
arising at the next order, treating them as the same order as the loop corrections from the LO 
Lagrangian. The most general Lagrangian at the next order is given by 



2 p 

£(4) =5i tr { a^Sta'^s)] + ^tr { djl'^d^^] tr { d^T^'^d^'Y. 



4 



+ 



2 r 



tr (mQSt + SmQ^J ^ + ^^tr (^S^SS'^St) tr [rnqT^ + Sm^ , (1.23) 



4 

where the coefficients jBSj are related to the original LECs of Gasser and Leutwyler [70] by 

4=^f^ , ^2 = (-2^ , 4=^P+^?^ , h=i4^- (1.24) 



n 



The loop effects typically give contributions proportional to y (4„-/)^ ) ' '^tiere n is the number 
of loops. Naive dimensional analysis and phenomenological evidence suggest that the chiral 

^The choice of the real, diagonal quark mass matrix is not necessary but is the most natural choice for obvious 
reasons. It is this choice, which then leads dynamically to the vacuum energy being minimized for E oc 1. 
However, any complex choice is allowed, and explicitly breaks the chiral symmetry group. The different choices 
are all related by chiral transformations. 

*We see that the square of the meson mass is related linearly to the quark masses, and the chiral condensate, 
in agreement with the Gell-Mann-Oakes-Renner prediction |89| . However, it is only recently that this relation 
has been phenomenologically confirmed by precisely determining the LECs of the corrections to the pion mass 
which arise at the next order in the chiral expansion 

|M| HH El EH El EH- This rules out the alternative 
description of nature using a generalized chiral perturbation theory in which the chiral condensate is not the 
leading order parameter of chiral symmetry breaking 
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symmetry breaking scale is A^^ ~ 47r/, and thus the chiral corrections are expected to scale as 
powers of the expansion parameter, 

^ ^ ^ (1 25) 

(4^/)2 (47r/)2 ^ 

This is the promised EFT expansion parameter, and we can see that in the real world with 
~ 135 MeV, this is expected to be a good expansion, with corrections being on the order of 

1% for SU (2) xPT. This will not be as good with the explicit inclusion of the strange mesons, as 

the kaon mass is significantly heavier than the pion masses. We will come back to this shortly. 
The Lagrangian H1.19() can also be used to compute the leading nn scattering amplitudes, 

which give the prediction of the scattering lengths by Weinberg using current algebra llODj . 

with the 1 = 2 vrvr scattering amplitude, for example, given by 

2 3s - 6ml 

(1.26) 

We shall discuss the tttt scattering in much more detail in Chapters [5] and IHl 



Three Flavors 



As will be a recurring theme in this Chapter, the formalism for the three flavor theory is identical 
to the two-flavor theory, the difference being the group algebra. There are some simpliflcations 
which occur in the two-flavor theory, which do not hold for the three flavor theory. These will 
become clear in this section. 

The LO Lagrangian, Eq. p.l9|) . takes the same form as in the two flavor theory, with the (p 
fleld now defined as 

K+ \ 



V2 V(i ' 



V 



vr 



K- 



(1.27) 



In the isospin limit {rriu = rrid), the vr" and rj fields do not mix, and the masses of the 
pseudo-Nambu-Goldstone mesons are given by 



rriK = Bo{mu4 + m^) 

2 Bq 
m, = -(m„ + m, + 4m. 



(1.28) 



The next difference from the two-fiavor theory occurs when constructing the NLO chiral 
Lagrangian. It has the same general form as Eq. H1.23() . but there are more simplifications 
one can make in the two-fiavor theory, as SU{2) has a real representation. The most general 
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Lagrangian consistent with the symmetries is given by |71j 



C^^^ = Li [tr ( d.^d^'T)'] 



+ L2 tr ( daT.duT.'^ ) tr ( d^'^d^'T) 



+ L3 tr + 2Bq L4 tr (^S^SS^S^) tr ^TTlgS^ + StTIq 

tr (m-gS^ + Sm^ 



+ 45^ L7 



+ 4Bo' ^6 



tr (^mqYj^ — Tim}(^ + 4i?g Lg tr ('rngS^rngS^ + SmQSm^ 



(1.29) 



One can match this SU{'i) Lagrangian onto the SU{2) Lagrangian of Eqs. (|1.19|) and ()1.23p . 
and learn that the SU{2) LECs depend logarithmically upon the strange quark mass. This is 
important for comparing SU{2) xPT calculations to lattice QCD simulations which include the 
strange quark in the simulation. We now move on to show how to include heavy baryons into 
the chiral Lagrangian. 



1.2.2 Heavy Baryon Chiral Perturbation Theory 

To systematically include the baryons into the chiral Lagrangian, we use heavy baryon chiral 
perturbation theory (HB^PT) first formulated in Refs. 27, 221. The idea is based upon heavy 
quark effective theory (HQET) jlOU ll()2j in which ones performs a velocity dependent phase 
transformation of the Lagrangian to separate the hard and soft modes of the theory. This 
amounts to studying the fluctuations about the heavy mass parameter. The velocity dependent 
fields are 



iMv-x 



■■Bix), 



(1.30) 



where M is the mass of the heavy field and is the four- velocity of the baryon, B. The 
momentum of the baryon can be parameterized as 



(1.31) 



where /c^ is the off-shell momentum of the baryon. The effect of this parameterization is to 
transform the standard Dirac Lagrangian for the heavy baryon of spin-^ in the following way. 



B {i^ - M) B 



B,i^B^ + 0[ ^ 



(1.32) 



Prom Eq. (|1.3U|) . it is easy to verify that derivatives acting on B^ bring down powers of the off- 
shell momentum, k, rather than the full baryon momentum, p. Thus, higher derivative operators 
in the EFT are suppressed by powers of the heavy mass M, and a consistent power counting 
emerges.^ Heavy baryon is applicable in the limit that the pion masses, momentum 



With the standard Dirac Lagrangian for the baryons, the power counting problem arises when one considers 
loop graphs and higher dimensional operators of the effective theory. The baryon mass is not small compared to 
the chiral symmetry breaking scale and thus one can not power count loop graphs with internal baryons, as they 
will be proportional to powers of M/A^ ~ The same problem occurs for higher dimensional operators, 

as one must include all operators with arbitrary powers of A'l/A^, which are not suppressed compared to the 
leading mass operator, B M B. For an alternative approach to resolving this problem, see Refs. |1U3I IIU^ . 
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and the off-shellness of the baryon are small compared to the heavy baryon mass (which is 
approximately the chiral symmetry breaking scale), 

m,r q V ■ k , ^ 

1-33 

M M M ^ ^ 

We also must include the lowest lying spin-| baryons in the theory as they play a very 
important phenomenological role in nucleon (or octet baryon) properties [23 EHl 11051 11061 

una hhhi unni ima cni ina inni uni insi inni nia ehi inni • However, there are questions 

which arise because these states are unstable, undergoing strong decay with a lifetime of r ~ 
10~^^ seconds (except for the Q,~ which decays weakly), and thus one can not define unique 
S-matrix elements for them. The effects of these states are not appropriately included by only 
absorbing their contributions to observable quantities into the LECs of the nucleon (or octet 
baryon) chiral Lagrangian for energies on the order of the pion mass. This is because the mass 
splitting between these resonant states and the spin-^ baryons is small compared to A-^, and 
thus they are easily excited. This leads to large contributions to the spin-^ baryon observables 
from these spin-| baryons, and more importantly, in scattering processes, a rapid oscillation of 
the scattering phase shift for energies near the mass of the resonance. Moreover, for the vast 
majority of lattice QCD simulations to date, the pion masses are heavy enough that the spin-| 
states are actually stable. We will come back to this issue in Section [2.1.21 

The inclusion of the spin-| fields into HBxPT requires the addition of an additional mass 
parameter, the mass splitting between the T and B fields in the chiral limit, A = Mt—Mb- This 
mass parameter is dynamically generated from the interaction of the quarks with the gluons, 
and is independent of the quark mass. The need for A arises from the velocity transformation 
of Eq. ()1.30|) . as one can not simultaneously remove both the nucleon and the delta masses (or 
octet and decuplet baryon masses).^'' To see the explicit effect upon the HBxPT Lagrangian, we 
will now focus on the 2-flavor theory, then the 3-flavor theory and lastly the partially quenched 
versions of these theories. 



Nucleons and Deltas 

Heavy matter fields (nucleons, deltas, heavy mesons etc.) transform under the vector subgroup 
(parity even subgroup) of the full chiral group. For example the nucleon fields, which are a 
doublet under SU{2)v^^ 




(1.34) 



transform as 

Ni ^ Uij Nj . (1.35) 



^"This can be accomplished for the full octet or decuplet masses, as in each case there is a symmetry protecting 
their common masses in the chiral limit. However, with the inclusion of both sets of fields simultaneously, 
there is no symmetry relating Mb to Mt, and thus one must choose a mass to phase away, the simplest choice 
being the lightest mass. We note that in the large Nc limit, A 0, and thus the B and T fields become 
degenerate |120| . This is not relevant to our current discussion however. 

^^Here, and in the rest of this work, we always drop the subscript, "v" on the heavy baryon fields, as we never 
use the fuU field defined in Eq. lll.3Ut . 
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When including heavy matter fields in xPT it is useful to introduce the field 

e = ^, (1.36) 

where S is defined in Eq. ()1.16() . Using the chiral transformation properties of S, Eq. (|1.17|) . 
one can deduce that under an arbitrary SU{2)l SU{2)ji transformation, 

LiU^ = U iR\ (1.37) 

where L,R are elements of SU{2)l,r respectively and U is an element of the chiral group. In 
the case of a vector transformation, L = R = V , then U = V as well. Otherwise, [/ is a 
complicated function of the pion fields in T, as well as space-time. It is also now clear that the 
choice of nucleon (or other matter field) is not unique, as the field 

N' = ^N^LN', (1.38) 

transforms identically to the field in Eq. (|1.35|) under the vector subgroup [HZ[ IHHl 1^ ESI • Thus, 
any choice of nucleon field which transforms as Eq. H1.35() under the vector subgroup is a good 
choice. The difference between these choices is simply a field redefinition of the nucleon field, 
and must give the same S-matrix elements. When constructing the most general Lagrangian, 
it is useful to introduce the chiral vector and axial-vector fields, 

A^=l{Cd^^^ -C^d^i) ^UAf^uK (1.39) 

The most general nucleon pion Lagrangian is then constructed from the and fields. At LO 
this results in 

C = NivDN + 2aMNM+N + 2 aMNNii{M+), (1.40) 
where the chiral covariant derivative is given by 

{D^N), = d^Ni + v/ Nj. (1.41) 

In the isospin limit, oc lsu{2) aiid thus the Lagrangian reduces to 

C = Niv ■ DN + 2dM'NNti[M+) (1.42) 

where ctm = ^oa/ + cm- For generality, we will always use Eq. (|1.4U|) when computing the 
masses but will often specialize to the isospin limit. 

We also want to include the delta fields, which are spin-| and isospin-|. To incorporate 
the spin degrees of freedom we use a Rarita-Schwinger field jl21| . T^, satisfying the constraint 
'jfiT^ = 0. To include the flavor degrees of freedom, we embed the delta fields in a rank-3 
totally symmetric flavor tensor, which transforms under the chiral group as 

Tijk^U,'U."'Uk''Tiran, (1-43) 
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with the normalization, 

rni = A++ , Tu2 = ^A+ , ri22 = ^A0 , Tasa = A" . (1.44) 

The full LO Lagrangian including the delta fields is 

C=NivDN + 2 au'NM+N + 2 cjA/lViV tr(7W+) 

-iT^'v DT^ + Ao T^T^ + 2-fM T^" M+ - 2a m T^T^ tr(A^+) (1.45) 

Here, Aq is leading order mass splitting between the nucleons and deltas in the chiral limit. 
There is some ambiguity in the definition of Aq and what we denote as the renormalized mass 
splitting in the chiral limit, A. We will clarify this ambiguity in Section [2. 11 This Lagrangian 
gives rise to the velocity dependent nucleon and delta propagators 

GN{k) 



k ■ V + ie 



gi^\k) = - — — -, (1.46) 

where P^^'^ is a spin projector which projects out the spin-| components of the Rarita-Schwinger 
field, and is given in d-dimensions (in the heavy baryon formalism) by 

pt^y = yt^v'^ _ _ -^S^" S" , (1.47) 
d-l ^ ' 



where 5^ is the Pauli-Lubanski spin vector [22| EH] . 

The chiral covariant derivative acting on the delta fields is 

(Z)^r)i,fc = a^r,,fe + (vn/'ri,,^ + (v^)/'r,,.,. + {v^k^'T^.y • (i.48) 

We also must include higher order terms which include the coupling of the nucleon and 
deltas to the axial field (these are the operators which give rise to the "pion cloud" of the 
nucleon) , 

C = 2gANS-AN + g^N {t'" A^^N + NA^'T^) + 2gAAT'' S-AT^. (1.49) 

This is the necessary formalism to compute many baryonic observables to LO and NLO in 
the chiral expansion. Before applying this formalism to the baryon masses in Chapter [21 we 
shall first discuss the 3 flavor theory and the partially quenched extensions of these theories. 
Most of the formalism caries over from the 2 flavor theory to the others, and so we shall only 
highlight the differences as compared to the above discussion. 



Octet and Decuplet Baryons 

When considering the octet and decuplet baryons jl22| 112,'^ . the only technical difference from 
the two-flavor theory is the group structure. There are some simpliflcations which occur for 
SU{2) as it has a real representation, which do not hold for SU{3). These will become clear 
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in this section. Another difference is that the inclusion of the strange quark, and the resulting 
mesons and baryons, makes the convergence of the HBxPT more questionable. For two-flavors, 
the expansion parameter is m-Tr/A^ ~ 15% for the physical pion mass, whereas rriK/^x ~ 40%. 
In this sense it is more important to get a precise determination of the SU{'i) HBxPT LECs to 
really test the convergence of the theory. 

The octet baryons transform as an 8 under the SU{2>)v flavor group, 

B^UBU^ (1.50) 
and can be expressed in a traceless 2-dimensional matrix, 

B=\ S- ^A.-ljEO n . (1.51) 

The spin-| decuplet baryons are still constructed with a Rarita-Schwinger field and a totally 
symmetric flavor tensor, T^jf^, and the same normalization of Eq. H1.44|) . with T113 = -^S*"'" 
etc. The decuplet transforms as a 10 under the SU{2>)y flavor group, as in Eq. (|1.43|) . except 
the flavor indices run from 1-3 instead of 1-2 as in the two-flavor theory. The most general LO 
Lagrangian for the octet and decuplet fields is^^ 

C = ti (BivDB) + 2bD tr {B{M+,B}) + Ibp tr (B [M+, B]) +2aMtr(BB) tr(A^+) 
- (T^ [ivD- Ao] T^) + 27M {T" M+ T^) - 2a m (T^T^) tr iM+) . (1.52) 

The chiral covariant derivative on the decuplet field has the same form as in Eq. ()1.48|) . while 
on the octet fields, is defined as 

D,,B = d^B+[V^,B] . (1.53) 

The leading axial coupling of the mesons to the baryons is given by 

£ = 2D tr (BS^" {A^, B}) + 2F tr (BS'' [Af,, B]) + 2n {T'" S ■ A T^) 

+ C (T^ A^,B + BA^ T>') . (1.54) 

Matching onto the two- flavor theory of Eq. p. 49)1 , one learns 

gA = D + F , fifAA = T-i , g/^N = C . (1.55) 

We see that the SU{2>) HBxPT Lagrangian, Eqs. (rr52|l and JTMl), has more operators involving 
the spin-^ baryons than the corresponding Lagrangian in SU{2). These are the differences 
arising from the group structure we have mentioned. Using the Lagrangian, Eq. 11.521 one can 
confirm the Gell-Mann-Okubo (GMO) mass relation amongst the octet baroyns |122| I124j . 
which is violated at the one- loop order in HBxPT from the interactions in Eq. 1)1.54(1 . We now 
move on to describe the partially quenched theories. 

^^We stress here that sometimes we wiU use the same notation for LECs in the three-flavor as in the two-flavor 
theory, however, they are different, and one should not confuse them. This problem could be avoided with 
larger alphabets. 
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1.2.3 Partially Quenched Chiral Perturbation Theory 

Partially quenched chiral perturbation theory (PQxPT) is the low energy EFT of partially 
quenched lattice QCD simulations, as described in Section 11.11 For LQCD simulations in 
which all the quark masses of the sea and valence sectors only perturbatively break chiral 
symmetry, the sicknesses which arise from the lack of unitarity can be systematically studied and 
understood with PQ^PT^^ [T^ [Ml 1123 ESI 1121 11211 11301 II3I1 HSl CSSl II311 • To construct 
the chiral Lagrangians in the previous sections, we made one assumption about the symmetry 
breaking pattern induced by the QCD vacuum, that the chiral symmetry group is broken down 
to the vector subgroup. To construct the partially quenched effective field theories, we will have 
to make a similar assumption and additionally one other assumption. The second assumption is 
that when all the quark masses can be treated perturbatively, the unitarity violation introduced 
by having different masses in the sea and valence sectors of the PQLQCD action is completely 
encoded in the low energy effective field theory, partially quenched chiral perturbation theory 
(PQxPT). This assumption is motivated by the fact that in the limit the sea and valence quark 
masses are degenerate, one recovers QCD. As with the assumption about the chiral symmetry 
breaking induced by the QCD vacuum (and the PQQCD vacuum), this assumption can not be 
proven, and so must be checked by comparing PQxPT expressions to PQQCD lattice correlation 
functions, and verifying wether the predicted behavior is observed. To date, there have been 
no signatures in PQQCD lattice simulations that PQxPT is somehow failing to capture the 
unitarity problems introduced by the sick theory. The real test however, will come in the next 
few years as the pions masses are pushed further into the chiral regime |5H I52j . 

We will now construct the PQ chiral Lagrangian. Where possible, we will be general and 
use A'^ for the number of valence quarks and Ng for the number of sea quarks, and will provide 
specific examples to help clarify certain points. To cancel the contributions from closed valence 
quark loops, we introduce fictitious ghost quarks, which are spin-^, but given Bose-statistics.^^ 
The mass of the ghost quarks is set exactly equal to the corresponding valence quark mass, 
resulting in a complete cancellation of the valence loops (which is manifested as a ghost-quark 
determinant identical in value but inverse to the fermionic determinant of the valence quarks 
in the path integral formulation of the theory). The quark content of the PQQCD Lagrangian 
is 

CPQ= Q' {iIp-mQ)^'' Qk. (1.56) 

As in Eq. 1)1.10(1 . for mq « Aqcd, the PQQCD Lagrangian exhibits an approximate SU{N^ + 
Ns\Ny) (gi SU{N^ + Ns\Ny)ji graded [lH6j chiral symmetry. One then assumes the PQQCD 
vacuum spontaneously breaks this symmetry down to the vector subgroup, giving rise to {2Ny + 

^^In the author's opinion, this is a great example of the power of effective field theory techniques. Using 
EFT, one is not limited to studying rigorously defined quantum field theories, but can also study theories in 
which some of the cornerstone foundations of a quantum field theory are broken, eg. unitarity. One can most 
likely never rigorously prove that partially quenched simulations away from the QCD limit do not produce 
uncontrolled effects, but partially quenched EFT is the only tool we have to study these systems, and to date, 
they have very successfully described partially quenched LQCD simulations. 

^''There is an alternative method to partially quench which is known as the replica method |135| . 
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Ns)'^ — 1 pseudo-Nambu-Goldstone mesons. Then, following the construction of the pions 
in Section II. 2. 11 we construct the partially quenched mesons representing the lowest energy 
excitations about the PQQCD vacuum. 



Mesons 



At leading order, the form of the partially quenched chiral Lagrangian for the mesons has almost 
the exact same form as in the SU{2) or SU (3) theory, with trace operators over flavor indices 
replaced with super-trace operators. In addition, there are operators involving the singlet field, 
$0 oc str($) which contribute to the leading two-point functions for the flavor diagonal meson 
fields. The Lagrangian is 

£ = ^str (a^Sa^St) + ^^^str (mgSt + Sm^ + a$ a^^o d^" $o - ml<^l ■ (1-57) 
The S field here contains many more mesons than in the two or three flavor theory, 



exp 



/2i$ 



M 

X M 



(1.58) 



The matrices M and M contain bosonic mesons while x contain fermionic mesons. 

Specifically, 
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(1.59) 



The upper A'^j, x block of M contains the usual mesons composed of a valence quark and 
anti-quark. The lower Ng x Ng block of M contains the sea quark-antiquark mesons and the off- 
diagonal block elements of M contain bosonic mesons of mixed valence-sea type. Expanding the 
Lagrangian (|1.57j) to quadratic order in the meson fields, one finds that for a meson composed 
of quark and anti-quark of flavor Q and Q' , the LO mass of the mesons is given by 



BoiruQ +mQf) . 



(1.60) 



As in QCD, PQQCD has a strong axial anomaly, and thus the mass of the singlet field, mo 
can be taken to be on the order of the chiral symmetry breaking scale, niQ Ay.. In this limit, 
the two point function between the flavor diagonal mesons, the r] fields, deviate from a simple, 
single-pole form. The momentum space two-point correlator between the ija and r]i, fields is 
found at LO to be [HHl 
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where 

+1 for a = valence or sea quarks 



1 for a = ghost quarks. ^''''^ 

In Eq. (|1.61j) . k runs over the flavor neutral states {(j)jj,4>ii,4'rr) and k' runs over the mass 
eigenstates of the sea sector. Here, we wish to show how the unitarity violations and partial 
quenching can be parameterized. To the author's knowledge, the following discussion was first 
brought up in Ref. (133. 134 as a means to quantify these unitarity violations in the two-nucleon 
system. In Ref. jl37j . this discussion was then expanded to the meson sector for both the two 
quark flavor as well as three quark flavor theory. There is additionally a very similar discussion 
in Ref. [SSI regarding the same issues. We refer the reader to Chapter El for the full details, 
but here highlight the main idea. When computing correlation functions with the partially 
quenched (or mixed action) theory, one can use a basis of fields other than the rjuu, ildd etc. 
which naturally arise from Eq. (|1.57j) . For example, in a theory with two valence quarks, and 
two sea quarks, in the isospin limit, one can instead work with the fields, 

7r° = ^ (r/n - %) , ^ = + ^'^) • (^•^^) 

The usefulness of this basis choice becomes very clear when one constructs the leading propa- 
gators for these fields 

GAp^) = -2 V— -' (1-64) 



1 1 



p2 _ _j_ (•p2 _ j^2 _|_ ^g^2 

T~2 2~~i ■ (1.65) 
[p'^ — + ley 



We see that the propagator for the vr'' has the exact form as in xPT, which is not surprising as in 
the isospin limit, the -k^ is an isospin 1 particle, and thus has no overlap with the vacuum (times 
75). However, the r/-propagator is proportional to the partially quenched breaking parameter 

^PQ = ""^La ~ ""^mience ' (1.66) 

where rusea is the mass of a pion composed of two sea quarks and rriyaience is the mass of 
a pion composed of two valence quarks. The double-pole structure of the ?7-propagator is a 
consequence of the unitarity violations introduced in the partially quenched theory. In the 
theory with two valence and two sea quarks, all unitarity violating terms can be traced back to 
loop graphs involving this ?7-field. We explicitly see here that these effects vanish in the QCD 
limit, 

'^sea — '^valence- similar change of basis can be performed for the two flavor theory 
away from the isospin limit, and also for the theory with three valence and sea quarks, although 
in these cases, the algebra is more involved. However, this parameterization, Eq. (|1.65|) and the 
equivalent parameterizations for the more involved theories provides a means to quantify the 
partially quenched breaking, or unitarity violations which arise in these sick theories. We will 
demonstrate a clear example of the usefulness of this in Chapter EJ 
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Baryons 

In PQxPT the baryons are composed of three quarks, QiQjQk, where i — k can be valence, 
sea or ghost quarks. One decomposes the irreducible representations of SU{Ns + Ny\Ny)v into 
irreducible representations of SU{Ny)yai "X" SU{Ns)sea <^ SU {Ny)ghost ® U{\). The method for 
including the octet and decuplet baryons into PQxPT is to use the interpolating field [138|I130] : 

- {QrQf'Qr - QrQrQt'") ^auci^u. {i.m 

We then require that when the indices z, j, /c, are restricted to the valence sector, Bijk = Bijk, 
where 

Bijk = ^ {eij Nk + Nj) for SU{4\ 2) , 

Bijk = ^ {eiji B^ + eiki S/) for 5C/(6|3), (1.68) 

where Ni is defined in Eq. ()1.34|) . B^ is defined in Eq. H1.51() and e^j and ejjfc are the completely 
anti-symmetric 2 and 3 dimensional tensors respectively. Thus the octet baryons are contained 
as an (8, 1, 1) in the 240 representation of SU{?,)yai ® SU{3)sea ® SU{3)ghost ® ^(1) [13111201, 
while the nucleons are contained as an (2, 1, 1) in the 70 representation of SU {2)vai '^SU (2)sea<^ 
SU (2) ghost jl31j . In addition to the conventional baryons composed of valence quarks, 

Bijk also contains baryon fields composed of sea and ghost quarks. In this paper we only need 
the baryons which contain at most one sea or ghost quark, and these states have been explicitly 
constructed in Ref. ,139, i3Q. for SU{6\3) and in Ref. [HIl for SU{A\2). 
Under the interchange of flavor indices, one finds |138j : 

Bijk = {-f^'^^'^'Bikj and Bijk + {-f^'^^'^^B.ik + {-f+Wmn.+mv.Skji = , (1.69) 

where the commuting or anti-commuting nature of the quark flelds are handled with the grading 
factors which are defined as 

f 1 for = valence or sea quark 

% = W r , 1 . 1 • (1-70) 

[ (J tor k = ghost quark 

Similarly, one can construct the spin-| baryons. The decuplet baryons are embedded as an 
(10, 1, 1) in the 138 representation of SU{3)yai ® SU{3)sea ® SU{3)ghost ^ U{1) ^M^M and 
the delta baryons are embedded as an (4, 1, 1) of the 40 representation of SU{2)yai'S'SU (2)c,ea® 
SU{2)ghost ® U{\) |131j . An interpolating field for the spin-| baryons is 

We require that Tijk = Tijk defined in Eq. H1.44() . when the indices k are restricted to the 
valence sup-space. In addition to the conventional decuplet resonances composed of valence 
quarks, Tijk contains fields with sea and ghost quarks, which have also been constructed in 
Ref. [13311301 for SU{Q\3) and Ref. [Ml for SU{A\2). Under the interchange of fiavor indices, 
one finds that 

T^Jk = {-)'^'''''^T^^k = H'-^'^^'^rik, . (1.72) 



24 



Under SU {6\3)v, both Bijk and Tijk transform as ^ 

Bijk ^ (-)'?'(''^+''-)+('''+^-)('f=+''")c/.'c/^™C/;,"^i^„ . (1.73) 

Leading order partially quenched heavy baryon Lagrangian 

The form of the heavy baryon Lagrangian in PQxPT is independent of whether one is working 
with the 5J7(4|2) or the 5C/(6|3) theory, at least through ©(m^). Therefore, we shah construct 
both Lagrangians simultaneously, but it is important to remember that the values of the LECs 
will differ between the two theories. This is because in the 5C/(4|2) theory, the effects of the 
valence, ghost and sea quarks are implicitly contained in modifications of the S'C/(4|2) LECs. 
One difference between PQ theory and the unquenched theories is the addition of the grading 
factors which keep track of the commuting or anti-commuting nature of a given baryon. In the 
following Lagrangian, we shall use the flavor contractions defined in Ref. jl38j . For a matrix 
Y^j acting in flavor-space, the contractions are 



{BYB) =B'''Y,'Bi,k , {r^YT^) = T^-'^V, ' r^,,,^ , 

(BBY) = (-)('?«+'v.)('?^+''0 b"'' Y,'B^,i , (BYf' %) = b"'' {Y%' T^,ijk • (1-74) 

Using these flavor contractions, the LO partially quenched heavy baryon Lagrangian is given 

by [123 csni inii USE] 

£PQ = (BivDB) +24^^^ {BBM+) +2[i^^j^^ (B M+B) + 2a[^j^^ {BB)stT{M+) 
-{T'-iivD-A] T,) + 27$;^) {T^M^T,) - 2al^«) [T^T,) str(A^+) 
+ 2a(^«) (BS^'BA^,) + 2(5^^^'^ {BS^A^^B) + 2n^^^^ i^''S^'A^%) 

+ y|c(^«) [ {TA,B) + {BAuT-) ] . (1.75) 

One can then determine the relation between these PQ LECs and the LECs in the two and 
three flavor theories by restricting the flavor indices in Eq. (|1.75|) to the valence sector. This 
leads to the relations between the SU{2) Lagrangians, Eqs. (|1.45|1 . (|1.49|1 and the SU{A\2) PQ 
Lagrangian, 

«M = ^a^^^ - , 7^^^ = 7^^ , 

(2) _ (412) 1 (4|2) 2 (4|2) (2) _ (4|2) 
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Similarly, when matching the SU{6\3) PQ Lagrangian to the SU{3) Lagrangians of Eqs. (|1.52j) 
and (|1.54|) . one finds the relations, 

V = ^(5ar + 2;3if') 

^(S) _ (613) , 1 (613) 2 ,613) 
Z)=l(«(6|3)_ 2^(613)^ 
F = _L (^5«(6|3) + 2/3(613)) _ 

With this formalism, we shall now proceed to further develop the Lagrangians as needed, 
and apply these effective field theory techniques to understand various baryonic observables. 
We begin with a detailed look at the masses of the lowest lying spin-i and spin-| baryons in 
HBxPT and the partially quenched analogues. 



(3) 

Tm 



^(3) _ -(6|3) 

c = -£:(6|3) , 



(1.77) 
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Chapter 2 
BARYON MASSES 

Perhaps the simplest question we can ask is what is the spectrum of hadron states which 
emerge from QCD? A quick glance at the Particle Data Book |14()j shows a very rich spec- 
trum. A significant understanding of these hadrons can be made by using the approximate 
5C/(3)Fiavor Symmetry of QCD, the most famous example being the Gell-Mann-Okubo mass 
relations amongst the octet baryons jl241 11411 1122j and the Gell-Mann-Oakes-Renner relation 
for the octet mesons j89l . However, we would like to be able to make a rigorous connection 
between the QCD Lagrangian and the observed particle spectrum, to arbitrary precision. This 
can be accomplished for the lightest hadrons (with a given set of quantum numbers) by utilizing 
lattice QCD in conjunction with chiral perturbation theory, as described in Chapter ^ 

We are interested in extending our understanding of the masses of the lowest lying baryons 
of spin-^ and spin-| to NNLO, or 0{m^)^ in the chiral expansion. This is important for several 
reasons. Firstly, as HB^PT is an asymptotic series, there is no guarantee that it is a convergent 
series, and so it is always good to push the order to which we know the theory to higher orders 
and test the convergence. Moreover, the expansion parameter of HB^PT is not as good as 
xPT with purely mesons, and so the expansion is expected to have larger corrections order 
by order. Thus, to do precision mass calculations (or precise determinations of any baryon 
observable) it is necessary to compute to higher orders in the effective expansion as compared 
to mesonic quantities. Therefore, the focus of this chapter is the computation of the ©(m^) mass 
corrections to the above mentioned baryons, using heavy baryon and partially quenched 
HBxPT. 

To accomplish this task, we must first construct the heavy baryon chiral Lagrangian to the 
next relevant order, beyond the known orders of Eqs. H1.45() and H1.49() . We first do this for 
the two-fiavor chiral Lagrangian (up and down), and then for the theory with three flavors (up, 
down and strange), and finally extend both of these Lagrangians to their partially quenched 
versions. This chapter is based upon the work in Refs. |142| I143| I144j . 

2.1 Nucleon and Delta Masses in SU{2) xPT 
Review of nucleon mass in xPT 

We begin by reviewing the nucleon mass to NLO in the chiral expansion, to set conventions 
and to motivate the need for knowing the NNLO correction to the baryon masses. It is useful 
to parameterize the nucleon mass as follows, 

M^^ = Mo (A) - M« (A, /.) - m(?/') (A, /.) - M^^} (A, /.) + .. . (2.1) 

The leading contribution to the nucleon mass, Mq, is generated by the non-perturbative dy- 
namical interactions between the quarks and gluons and is independent of the quark masses; 



27 



this parameter is often referred to as the nucleon mass in the chiral hmit. The parameter, A, is 
the quark-mass independent sphtting between the delta mass and the nucleon mass, and is phe- 
nomenologically determined to be A ~ 293 MeV.^ This mass splitting arises dynamically from 
QCD and vanishes in the large Nq limit jl2Uj . The terms, are contributions to the mass 

of the i*'' nucleon of the order (m^)". In Eq. (|2.1() . fi is the renormalization scale which arises 
when regulating divergent loop diagrams. Throughout this chapter, we use dimensional regu- 
larization with a modified minimal subtraction (MS) scheme, in which we consistently subtract 
terms proportional to 

7i? + l + log47r, (2.2) 



A-d 

where d is the number of space-time dimensions. 

There are a few points worth noting regarding the chosen form of the mass expansion. First, 
the minus signs in Eq. 1)2. 1|) are purely for notational convenience. Second, the renormalization 
scale dependence appears at each order because we are treating the LECs as polynomial func- 
tions of A. This is motivated by our interest in comparing xPT expressions to lattice QCD 
correlation functions in which the quark masses (but not A) can be varied. Moreover, as we 
will show, because one does not have the ability to vary A, either in nature or with lattice 
QCD, there is some ambiguity in how to include the finite analytic contributions proportional 
to this mass splitting.^ We will now make this explicit by reviewing the chiral contributions to 
the nucleon mass. 

The leading quark mass dependence of the nucleon, comes from the diagram depicted 
in Figure 12.11 It is a simple tree-level diagram giving rise to the mass correction 

M^^^ = 2aMm + 2a M {mu + m^), (2.3) 
where the nucleon dependent contribution is given by 

m. = | . (2.4) 

[ nid, i = n 

The leading pion loop graphs depicted in Fig. I2.ir b) and I2.ir c) give rise to non-analytic 
dependence upon the quark mass. It turns out when using dimensional regularization, the loop 
graph with an internal nucleon is completely finite, giving rise to the self energy correction 

TfV' 

ST^N-ioop = -Stt g\ . (2.5) 
X 

^We note that the experimental mass splitting between the deltas and nucleons is determined to be Aphys ~ 
293 MeV fl401 . but that this value includes quark mass dependence as well. The difference, however, between 
using this value for A in xPT and the true chiral limit value of A is a higher order effect. 

^We note as an aside, that if we instead were to treat this dependence explicitly and expand the nucleon mass 
in powers of q, where A ~ g and ~ q^, the mass (and all baryonic observables) would then take the form 

itf (2) (3) (4) , 

Mn, = mo - my - my - my + . . . , 

with m^' as the renormalization-scale independent mass contribution to the i*** nucleon strictly of the order 
q". The parameter mo is nucleon mass in the chiral limit, but only to leading order in A, i.e. mo = Mo{A — 0). 
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Figure 2.1: Leading order and NLO quark mass dependence of nucleon mass. The LO con- 
tribution, diagram (a), is given by the Lagrangian in Eq. (|1.4U() . while the NLO pion-loop 
contributions, diagrams (b) and (c), arise from the axial coupling of the pions to the nucleons 
and deltas given in Eq. H1.49() . The solid black lines, thick-solid line (red online), and the dotted 
black lines denote nucleons, delta-resonances and pions respectively. 



The loop-graphs with internal delta states are more interesting and give rise to the renormalization- 
scale dependence at this order. The un-renormalized self energy correction arising from dia- 
gram E^c) is given by 
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where we have defined the function T by 
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This is slightly different than the function as defined in Ref. |13nj . Here, we have chosen to 
subtract the chiral limit value of the function, F(0, A,//), as defined in Ref. |13Uj . instead of 
absorbing this contribution into Mq (as defined in Eq. (|2.7|) . .F(0,A,/i) = 0). This has the 
slight advantage that the decoupling of the deltas in the chiral limit 145 is more transparent. 
After applying our renormalization prescription defined in Eq. 1)2. 2() . we observe that we can 
absorb all the analytic dependence upon A into a redefinition of the LECs of the LO chiral 
Lagrangian, Eq. (|1.4U|) . 
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Chiral corrections to the nucleon mass 
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Figure 2.2: Here we plot the quark-mass dependent (chiral) corrections to the nucleon mass 
in the isospin limit. The upper dashing curve (black dashing online) is the LO correction, 
The next two longer dashing curves (longer black and red dashing online) are the NLO 
corrections coming from the nucleon-pion and delta-pion loops respectively. The lowest curve 

(3/2) 

(green on-line) is the total NLO mass correction, Mj^ . The dark solid curve (blue online) is 
the total chiral correction, — Mq, through NLO. 



such that the nucleon mass to NLO is given by 
Mn, = M^{A) -2a^jmi-2 a^iA, //) (m„ + m^) 

X 

To help understand the importance of computing the next order chiral corrections, we will plot 
the contributions to M^r in the above expression, setting /i = 47r/. Here, we are not concerned 
with precision, but rather are interested in the qualitative features of this expansion (of course 
a thorough understanding of the convergence of HBxPT necessitates a precise determination 
of all the LEC's relevant to a given order). Therefore, we estimate our values for a^j and 

from Ref. [IDSj by matching onto SU{2), arriving at the values ~ —1 GeV~^ and 

cr^ -1 

~ —2 GeV . We also use the physical values of the nucleon-pion coupling, qa — 1.25, and 
the nucleon-delta-pion coupling, g^N — 1-5, as the difference between using these values and 
the values at this order are higher order than we work. The resulting contributions are show in 
Fig. 12.21 In this figure, we have plotted the NLO contributions arising from the nucleon-pion 
loops and the delta-pion loops separately as well as the sum. This is simply to demonstrate 
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the importance of including the delta as a dynamical degree of freedom in the effective field 
theory [27,. ,28^ , as can be seen by the comparable size of the nucleon loop as compared to the 
delta loop. The important feature this figure highlights is that for pion masses of m-,^ ~ 250 MeV, 
the NLO contribution begins to compete with the LO contribution, as seen by a flattening of 
the total mass correction to NLO, the solid (blue) curve in Fig. 12.21 This clearly demonstrates 
the need for computing the NNLO mass correction; to have confidence in the fit for values 
of the pion mass m-jr > 250 MeV. The precise value of the pion mass for which one needs 
the NNLO contribution to trust the expansion will be sensitive to the values of and a^j, 
which will in turn depend upon the size of the NNLO contributions. Thus a detailed fitting of 
these parameters, comparing the NLO mass formula to the NNLO mass formula will not only 
provide information about the value of the LECs but also will help determine for which values 
of the pion mass we can trust the convergence of the chiral expansion. The NNLO calculations 
are certainly necessary for the values of pion masses used in current lattice simulations, being 
typically > 300 MeV. 



Higher dimensional operators 

In order to compute the NNLO (©(m^)) correction to the nucleon and delta masses, we must 
first construct the relevant terms of the higher order chiral Lagrangian. It is convenient to 
split these terms into three different categories. The first set of higher-order terms are all 
operators whose coefficients are constrained by reparameterization invariance (RPI) tl46i| . and 
thus known. These terms must be included to insure the Lorentz invariance of HB^PT at 
0{1/M), where M ~ A-^ is the average nucleon mass in the chiral limit. The momentum of 
the nucleon, as given in Eq. ()1.31|) . is not a unique parameterization. In particular, under the 
reparameterization 

+ , k^k-e, (2.10) 

the nucleon momentum, = Mw^ + is unchanged. Reparameterization invariance is the 
requirement that the Lagrangian be invariant under such transformations, which ensures the 
Lorentz invariance of the theory to a given order in 1/M. Moreover, the implementation of 
reparameterization invariance jl46l I102j is a useful tool to generate terms in the higher order 
Lagrangian, as the reparameterization automatically generates the higher terms needed for 
Lorentz invariance. The result of applying these RPI techniques to Eq. ()1.45|) and using the LO 
equations of motion to eliminate terms |147l I148j is the following fixed coefficient Lagrangian 

jr = -N^N + T^^T, + gA(N'^^v-AN-NvA^-^N] 
2Mo 2Mo ^ I Mo Mo i 



+ 9.A|T''!^„.^r„-T%.^^T,.|, (2.11) 



where = — {v ■ D)^. Note that by combining the kinetic term for the nucleon in the 
above Lagrangian and the kinetic term in Eq. (jl.45|) . the nucleon propagator becomes 



(2.12) 
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as promised. 

The next set of operators are the standard terms in the 0{q'^) and 0{q^) Lagrangian jl49j 
with undetermined coefficients and which contribute to the masses at NNLO. Here we use 
different notation for the LECs than in Ref. jl49j . as in that work the delta's were not included 
as dynamical degrees of freedom, and thus the values of the LECs will be different. These 
operators are ^ 

C =-^\bfN MIn + bl^NNti{Ml) + bfN M+NtT{M+) + bl^NN MM+)f 

+ b^ NN tr{A ■ A) + b""^ NN tr{v -Av-A)]. (2.13) 



The LECs bf, b^, and b'"^ are all dimensionless. The choice in numbering the coefficients 
was made to be consistent with the partially quenched version of this Lagrangian jl42j , as will 
become clear in Section [2.31 We have not written down operators of the form, NM^ Ai^ N, 
even though they are formally of this order, as they do not contribute the nucleon masses until 
a higher order. Also in this set of higher-dimensional operators are the equivalent terms for the 
delta-resonances, 

+ tir^T'^tviA^'A,) + r^''(T; • A)/{v • A)/t^^j,, + f/T^T^'tv{v ■ Av ■ A) 
+ tfTf{M+M+V'Ti;^, + tfTf{M+y'{M+)/Ti:^,, + 4'T,T^tT{M+M+) 

+ tf {T^M+T>') tT{M+) + T^r^[tr(7W+)]2|, (2.14) 

and ah of the LECs , tf, t^^, and tf are dimensionless. 

The last set of higher-dimensional operators with undetermined coefficients all involve the 
nucleon-delta mass-splitting parameter A, which is a singlet under chiral transformations. Be- 
cause of this, the inclusion of the spin-| fields in xPT requires the addition of operators involving 
powers of A/A^. We shall not write out such operators explicitly. To account for the effects of 
these operators, all LECs in the calculation must be treated as arbitrary polynomial functions 
of A/Aj(, and expanded out to the required order. For example 



2 



A\ / A A 

^xJ V ^X K 



/ A\ / A A2 . ^ ^ 

QA^ gA\-^j =9A (^1 + gA,i J- + gA,2 jr + ■■■ ) • (^.15) 



■^One could also add operators of the form {^NS ■ AS ■ AN) , but it is straight forward to show that they are 
a linear combination of the operators in Eq. and therefore not distinct. 
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(a) 



(b) 



(d) 



□ ■ ■ □ O 

(e) (f) (g) (li) 

Figure 2.3: Diagrams depicting the various contributions to the NNLO nucleon mass correction. 
The black, thick grey (red onhne) and dashing hues correspond to nucleons, delta-resonances and 
pions respectively. The filled circles and filled squares are insertions of operators from Eq. (|1.45j) 
and Eq. H1.49() respectively. The open squares are RPI fixed operators from Eq. ()2.11|1 and the 
open circles are higher dimensional operators from Eq. (|2.13|) . 



Here we see the ambiguity in the treatment of A. In Eq. 1)2. 8|) we showed how to renormalize 
the LECs of the LO Lagrangian to absorb the A dependence, in particular the LEG aM was 
modified by a known term proportional to A. However, we see from Eq. ()2.15)) that we must 
additionally add a term to ctm linear in A but with an unknown coefficient. And in fact every 
LEG in the HB^PT Lagrangian has this unknown dependence upon A. However, determination 
of the LEGs in Eq. (|2.15)) require the ability to tune the parameter A, which we can not do in 
nature or in lattice QCD. 

To resolve this ambiguity, we advocate the following method. The only A dependence one 
should keep explicitly is that associated with non-analytic dependence upon m-jr, for example 
in the function J^{m.j^, A, fi); this will keep all the important effects of the delta-resonances, 
the large chiral contributions and the cuts in graphs where the internal delta's can go on shell 
(for example in nucleon Compton scattering, Ghapter ^ nucleon pion scattering, and pion 
photo-production off the nucleon). All the finite analytic dependence upon A, whether known 
contributions as in Eq. (|2.8() or undetermined contributions as in Eq. 1)2. 15() . should be absorbed 
in a renormalization of the LEGs. This was consistently done in Ref. jl30) to NLO, and in the 
following section we will demonstrate the consistency of this treatment in the nucleon masses 
through NNLO.^ We now proceed to compute the NNLO nucleon and delta mass corrections. 

2.1.1 Nucleon Masses through 0{rn^) 

The chiral correction to the nucleon mass at 0{m^) is a bit more involved than the LO and NLO 
contributions, receiving corrections from the Lagrangians, Eqs. ()1.45|) . ()1.49|) . ()2.1H) . ()2.13|) and 
combinations of the first two. The various contributions are depicted Figure E31 Before giving 
the result, we first will demonstrate the consistency of the treatment of the nucleon-delta mass 

*As this renormalization prescription is not guided by a symmetry, as is the case for the quark mass dependent 
contributions, it would be nice to check that this prescription holds for other nucleon observables as well, i.e. 
the specific renormalization of the A dependence done for the nucleon masses gives an exact cancellation of 
the analytic A dependence in another nucleon observable like the magnetic moment. 
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splitting parameter, A. To do this we will compare the mass correction of Figure E^fe-f), which 
involves the kinetic correction to the delta-resonances, to the contribution of Figure I2.2r c). 
involving the bare LO delta propagator and resulting in Eq. (|2.6|) . The mass correction from 
this diagram is 
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MoA2 



A-d 



+ 1 + log47r 



-A^ - 2A2m2 
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-J^(m^, A,/i) 



(2.16) 



Note that if we shift the nucleon-delta-pion coupling by a known amount, 

/ 5 A 



(2.17) 



such that the difference between the self energy corrections at NNLO and NLO with the modified 
dAN coupling (to higher orders in dependence of qan) becomes 



5S 



(2) 

A—loop 



5A (3/2) 

2Mo 



^9an 
MoA2 



d 



7£' + 1 + log47r 



8 a4 94 

— A A 

15 15 20 



-mjog 



. (2.18) 



Note that the important contributions from Figure [THte-f) are exactly reproduced (the dimen- 
sional regularization divergences proportional to Eq. (|2.2p and the A dependence associated 
with non-analytic m.,^ dependence as in ^(m7r, A, /i)). There are additional analytic contribu- 
tions, both finite and infinite, as well as a new non-analytic dependence which arise from this 
graph, but under the above stated treatment of the A dependence, it is redundant to keep all 
finite terms in Eq. (|2.16j) . We apply this analysis to all the nucleon mass corrections arising 
from the diagrams in Figure IT^ arriving at the following expression for the 0{mg) nucleon mass 
correction; 
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(2.19) 
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Table 2.1: The coefficients F^', and G^^ in xPT. Coefficients are listed for the nucleons. 







TT 
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2aM{2mu + md) 


auimu + 2md) 


|(7m„ + 2md) 


n 


2aMimu + 2md) 


aM{2mu + md) 


|(2m„ + 7md) 



Here we have used {rrii)'^ which is the square of the tree- level coefficients rrii appearing in 
Eq. (|2.4|) . Above the wavefunction renormalization Zj^. is given by 
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^[j(m^,A,f,)+ml]. (2.20) 



The coefficients in the NNLO contribution, namely -^i^S and are given in Table ITTI 
and depend on whether i = p or i = n. The equations above, Eq. ()2.19|) and Eq. ()2.2U|) . 
also employ abbreviations for non-analytic functions arising from loop contributions. The new 
functions is defined as 

J{m, A, /i) = 2A2 log f — ^ J + m2 log l-j 

- 2A./^^Hog ( ^ + ^^^^^ , (2.21) 

where ^7(0, A, /x) = 0, similarly to .^-"(0, A, /x). 

The expressions we have derived in this section, as well as those throughout this work, are 
functions of the quark masses; e.g. m-,^ above is merely a replacement for the combination of 
quark masses given in Eq. (|1.22)) . These expressions, thus require the lattice practitioner to 
determine the quark masses. In SU{2) xPT away from the isospin limit, there are no equivalent 
expressions in terms of the meson masses, as one cannot independently associate niu and 
to the pion masses, which are degenerate to the order we are working. One can only equate the 
average value, ^{mu + md), with the pion mass. Therefore one cannot plot the baryon masses as 
a function of the physical meson masses, unless one works in the isospin limit, m„ = m^. This 
is unlike the case in the isospin limit of SU{3), where there are only two independent quark 
masses, but three independent meson masses, and one can always convert from a quark mass 
expansion to a meson mass expansion via the Gell-Mann-Okubo relation.^ Lastly we remark 
that if one is using xPT to determine the quark masses from meson masses that are determined 
on the lattice, one must use the one-loop expression in xPT, else one looses contributions to 
baryon masses that are of NNLO. Thus the quark mass expansion is the one to use in SU{2) 
xPT, at least away from the isospin limit. 

^This problem can be avoided in SU(A\2) PQxPT, as there are more independent meson masses than inde- 
pendent quark masses, so one can algebraicaUy convert from the quark mass expansion to the meson mass 
expansion. For example, to leading order, rriu = 2^ (jn\ — m^^ + m^^) . This would require one to know the 
mass of the mesons made of one valence and one sea quark. But if one is interested in the non-isospin 517(2) 
limit of the PQxPT expressions, the problem is unavoidable. 
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We conclude this section by giving the nucleon mass in the isospin hmit, expressed entirely 
in terms of the pion mass, as this is the most relevant and useful formula for lattice QCD 
simulations today. The nucleon mass in this limit is 

M^ = M,{A)-a,, — -2a,,iA,^,) — -3ng^ — -— 



,2 



C^(A, /i) mi + C/^(A, /i) log ( ^ ) + Cf (A, /x) ml J{m^, A, /.) + 0{m['/'^) , (2.22) 



where 



+ 8A3 4S2A^ 2i?2A, 2BIK^ B^A^ + i?oA| + ^^^^ 



^2 



i3oA2 4A3 > MoA2 V 16 ' 2 / ' fioAr^'" ' """"'^ 



~~r1^ ' *^ + ^(^M + 1m - ^<^m) ■ (2-23) 

X 

Each LEG in the above coefficients is an implicit function of A as discussed in the previous 
sections. The LECs of Eq. 1)2. 23() additionally have renormalization scale dependence, which 
we have suppressed, that exactly cancels the scale dependence of the logs in Eq. (|2.22j) . At this 
order, in the isospin limit, there are these three linear combinations of LECs which contribute 
to the nucleon mass at this order, each of them containing LECs which are not well determined 
(if at all). We use naive dimensional analysis (NDA) [HH] to estimate the size of the new LECs, 
arriving at a central value set of ~ afj ~ 5^ ~ 6''^ ~ 1, and \h\j + \hf + ^bf + ~ 1. 
We then do an un-correlated variation of these LECs from — 5 < LEC < 5. The results are 
plotted in Figure [231 The important feature to note is that as can be seen by Figure |23fa), the 
NNLO mass correction, for all values of the LECs used, is smaller than the LO and NLO mass 
contributions, lending confidence to the chiral expansion for this observable. We anticipate that 
soon, all the LECs which contribute to the nucleon mass to this order will be determined from 
lattice QCD simulations, providing the first rigorous determination of HBxPT parameters at 
this order, allowing a chiral extrapolation to the physical pion masses and a prediction of the 
nucleon masses from lattice QCD. 

We now move on to compute the mass of the delta-resonances in the chiral expansion. 



2.1.2 Delta Masses 

Before embarking on the mass computation, it is first necessary to discuss what we mean 
by the mass of a state that can undergo strong decay. The first simple answer is that for 
sufficiently large pion masses, the delta-resonances are actually stable particles. This is because 
the nucleon delta mass splitting, as we have already mentioned, is independent of the quark 
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Nucleon mass contributions through O(m^) 



1 GeV 



Mn -Mq 
1 GeV 



Ara„ (LO) 



Anig (NNLO) 




(a) 



Nucleon mass to 0{m^) 



Amq+Arriq''^ 



. 15 

(GeV) 



(b) 



Figure 2.4: In this figure, we plot (a) the first three chiral corrections to the nucleon mass in the 
isospin limit. Although the NLO correction is not small compared to the NNLO contribution, 
we observe that for the range of LECs used, the NNLO corrections are definitely smaller than 
the LO and NLO corrections for the full range of pion masses plotted, and for most of the range, 
the NNLO contribution is significantly smaller. The grey band is an estimate of the theoretical 
uncertainty at the NNLO order, as there are still undetermined LECs at this order. In (b) we 
plot the total mass of the nucleon through 0{mg), and project the uncertainty in the NNLO 
contribution onto this plot as well. For comparison, we also plot the total mass contribution 
through NLO, which is also plotted in Figure [221 
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masses, and thus for mn > A = Mt — Mn, the delta is kinematicahy forbidden from decaying. 
This occurs for m-^ ~ 300 MeV where the convergence of HBxPT is questionable 52^ and 
thus in general it will be more difficult to extract observable quantities for the deltas than for 
nucleons.^ Nevertheless, the idea is to use the mass expressions derived in this section (or other 
observable quantities) to compare to lattice QCD simulations of delta observables for values 
of the pion (quark) masses such that the deltas are stable particles. One then can determine 
the LECs which describe the properties of these "resonance" particles and then analytically 
continue to the physical quark masses of nature. Using these techniques, for example, one can 
make a prediction with lattice QCD in conjunction with HB^PT of the width of the deltas 
(this is well determined experimentally, so this is really a benchmark test). The real test of 
this theory [271 128j however, will be to determine the LECs governing the delta properties by 
comparing HB^PT calculations to lattice QCD correlation functions, and then make predictions 
of delta observables which play a noticeable role in experimentally measurable quantities. The 
simplest quantity one can determine is the mass of the delta, which plays a large role in the 
majority of nucleon properties, see Figure 1221 for example. Thus in this section, we work out 
the mass of the deltas to 0{m'^). 

The computation of the delta-masses largely parallels that of the nucleon masses presented 
in the previous section. Therefore we will mostly present the results, highlighting the differences 
from the nucleon mass computation. Similar to the nucleons, the delta mass of the r*^ delta 
can be expressed in the chiral expansion as 

Mt, = Mo (A) + a + M^J^ (A, fi) + Mg^^^ (A, ^^) + M^^ (A, /i) + . . . (2.24) 

Here, Mq (A) is the renormalized nucleon mass in the chiral limit from Eq. (|2.1j) . and A is 
the renormalized nucleon-delta mass splitting in the chiral limit. Both of these quantities are 
independent of niq and also of the Tj. M^"'* is the contribution to the i^^ delta baryon of the 

order rri~^\ and ^ is the renormalization scale. ^ 

The diagrams relevant to calculate the delta masses to NNLO are depicted in Figure 12.51 
The leading-order contributions to the delta masses are 



M^T^ = - 7M mf - 2a M tr(mg), (2.25) 



where the tree- level coefficients mj are given in Table EH for the deltas T. 

The NLO contributions to the delta masses are more interesting. It is at this order that we 
first see the imaginary part of the delta mass (for m^r < A), and see that it decays. The mass 
corrections are depicted in Figure ES^b) and l2.5f c^. resulting in 

^(3/2) _ 25gi^ 3 _ 2g\^ _ 
~ 4327r/2"'" 3(47r/)2-^^"'"' ^'^>' 

^There are other tricks one can use in lattice QCD to stabilize the deltas for < A, but we will not discuss 
them here. 

^As with the nucleon masses, the renormalization scale appears in each term contributing to the delta masses 
because we implicitly treat the LECs as polynomial functions of the mass parameter A. If we expand out the 
LECs, the /i-dependence disappears at each order. 
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Table 2.2: The tree-level coefficients in xPT. The coefficients mf , {m^)J , and {mm')J are listed 
for the deltas T. 
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(a) (b) (c) 

Figure 2.5: Leading order and NLO quark mass dependence of delta-resonance mass. The LO 
contribution, diagram (a), is given by the Lagrangian in Eq. (|1.45|) . while the NLO pion-loop 
contributions, diagrams (b) and (c), arise from the axial coupling of the pions to the nucleons 
and deltas given in Eq. (|1.49|) . The diagram (c), with an internal nucleon, is the diagram 
which gives rise to the imaginary contribution to the delta mass at this order, for m-,^ < A. 
The solid black lines, thick-solid line (red online), and the dotted black lines denote nucleons, 
delta-resonances and pions respectively. 



where the function .F(m,rj A, /i) is defined in Eq. (|2.7j) . Note that in Eq. (|2.26j) . it is .F(m^, — A, /i) 
which appears, arising from Figure I2.5f c) with an internal nucleon loop. For m-j^ > A, this 
function is real valued,* but at m.,^ = A, this function develops a cut, and thus has a real and 
imaginary component for < A, and thus the deltas become unstable and decays.^ As with 
the nucleon masses, the LEG, Wm, becomes scale dependent to cancel the scale dependence 
introduced in Eq. (|2.26|) . The derivation of the scale dependence follows exactly from Eq. (|2.8j) . 



**For m-jr > 0, jr(m,r, A,/i) is always real valued. 

^In lattice QCD simulations, the geometry of the box will also play a role in the decay of the delta (or other 
resonance states) 150 . Because the lattice box only allows for quantized momentum modes, one can be in a 
region of pion masses where the delta would decay in infinite volume, but the geometry of the box (the box 
size) does not allow the necessary kinematics for the decay, and thus the deltas are stable. Another way to see 
this is that in finite volume, the imaginary piece of the function .7-"(m,r, — A,/i), comes from a series of poles 
determined by the kinematics allowed by the box, instead of a cut as in infinite volume. 
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After applying the renormalization prescription of Eq. (|2.2|1 . the LO delta LECs become 



Mo 
Ao 

7M 



Mo^(A) = Mo + 



9 



A2 ' 

X 



A = Ar 



9 



7m 



7M 



<^ M = M 



A2 

X 



+ 



A2 

X 



log 



(2.27) 



Here we see that the parameter Mq is the same as in Eq. ()2.8p . We restate that this was a 
choice to make this equal to the nucleon mass in the chiral limit. A different choice would 
simply result in a different value for A (but then also the mass splitting parameter in the J- 
function would be Mt — Mb instead of A). We also reiterate that the only important part of 
this renormalization prescription is the ^ dependence, as there is arbitrary dependence upon A 
in all these LECs just as with the nucleon LECs. 

The NNLO (©(m^)) corrections are given by diagrams similar to those in Figure but 
with an interchange of the nucleon and delta resonance lines for both the internal and external 
states. The mass correction is then 
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(2.28) 



Here we have used {m?)f , {mm')f to label the tree-level coefficients that appear in Table 
Above the wavefunction renormalization Zy. is given by 
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The coefficients in the above expression, namely Cj', and are given in Table ESI and 
depend on the delta state T. Given these delta mass expressions, Eqs. (|2.25j) . (|2.26|) and (|2.28j) . 
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Table 2.3: The coefficients Cj', FJ% and G^" in xPT, for the various delta states, T. 
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one can predict the width of the deltas with the knowledge of a small number of LECs, f,BQ, 
A, aM, ctm, IM, and most of which can be determined independently from the delta 
masses. The expression for the delta width is given to this order by 
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(2.30) 

Additionally, one can compute the quark mass dependent nucleon delta mass splitting, which is 
the quantity which should be used in place of A in the functions, J- and J ^ as for a given quark 
mass, this is the parameter which governs the stability of the deltas. Thus, the real component 
of the nucleon delta mass splitting is given through 0{m'^) by the expression, 
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Strong isospin breaking of the nucleon and delta masses in lattice QCD 

Having derived the nucleon and delta masses to NNLO in the chiral expansion, we now focus 
on the mass splittings between these states. To begin, we consider the nucleon mass splitting, 
which to our knowledge was first theoretically addressed in The degeneracy between 

the proton and neutron is broken by leading-order effects in the chiral theory, see Eq. ()2.3|1 . 
Beyond this order pion loops contribute, but to the order we are working, all the pions are 
degenerate, even away from the isospin limit. Thus the NLO contributions to the neutron and 
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proton masses are the same, see Eq. (|2.9|) . and the mass sphtting, M„ — Mp, is given to NLO 
accuracy, entirely by the difference of the LO mass contribution in Eq. (|2.3|) . and is hnear in 
rrid — mu- Any deviation from this hnear mass sphtting seen in lattice simulations of the nucleon 
masses should be a signature of the NNLO mass contributions and certain LECs that arise at 
this order. Additionally, the nucleon mass splitting can be enhanced from that in nature by 
increasing the quark mass splitting on the lattice, — rriu- This enhancement, combined with 
the vanishing of the NLO contribution to the mass splitting, provides us with a means of cleanly 
determining the NNLO nucleon mass contributions and isolating certain LECs arising at this 
order. These effects are normally obscured by the NLO contributions. 
We find the nucleon mass splitting is given to NNLO by-*^"^ 
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Bo l + T] 
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7] ml 



1 + 71 Bo 



m^ 
A2 

X 



8 ( C/A"Af + gAN 



y 

2 
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J{m^,A,fi) 2 
H — 09 AN 
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(2.32) 



where = An identical situation arises for the deltas, as their degeneracy is broken at 
leading order in the chiral expansion, while the next contribution to their splittings occurs at 
NNLO. Their mass splittings are given by 5Mt, which stands for M^- — M^o , M^o — 
and Ma+ - Ma++. We find, 
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27Af 1 
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J{m^,-A,fi) 2g\j^ 
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(2.33) 



where 5t^^ is given in Table ITU Considering these splittings, there are a few things to note. 
We know from experiment that the neutron is more massive than the proton, and we expect 
that this is mostly due to > m„. Similarly, we expect the deltas to follow a similar pattern 



^'^We note that during the completion of this thesis, the first determination of the strong isospin breaking in 
the nucleon masses from a fully dynamical lattice QCD simulation was completed 60 , although unfortunately 
there were not enough different pion masses used to make use of this full mass splitting expression. 
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Table 2.4: The coefficient 5t , which encodes the violation of the delta equal spacing rule at 
0{ml). 



6Mt 


6t^ 


Ma- - Mao 
Mao - Ma+ 
Ma+ - Ma++ 


2md 
m„ + md 
2m„ 



with A~ being more massive than A , and so on, although this is presently undetermined 
experimentally. Prom these expectations and the expressions above, it is expected that aM is 
negative and 7a/ is positive. Also, the delta masses, to a good approximation are expected to 
follow an equal spacing rule. In HB^PT, this rule is first violated at NNLO, and even then, 
only by the operator in Eq. ()2.14|) with coefficient . This coefficient can be isolated by taking 
successive differences. 

Lastly we should comment that in nature isospin violation in the baryon masses has another 
source of the same size as the NNLO chiral effects, namely electromagnetic contributions. In 
lattice QCD calculations one can turn off the electric charges of the quarks. This is the sce- 
nario for which our calculations are applicable. There have been a few lattice computations of 
electromagnetic contributions to hadronic masses |1521 11531 1154j . Por a recent discussion of the 
electromagnetic effects in hadrons, see Ref. |155j . and for a comprehensive phenomenological 
review, |156j . 



2.2 Octet and Decuplet Baryon Masses in SU{3) xPT 



In this section and the following section, Sec. 12.31 we wish to extend the baryon mass compu- 
tations to the SU{3) theory as well as the partially quenched theories, S'C/(4|2) and SU{6\3). 
There are many more baryons in each of these theories than the two nucleons and 4 deltas in 
SU{2). To highlight the differences between these theories we will only explicitly give the masses 
of the nucleons and deltas in the following sections and refer the reader to Refs. |142[ I143[ I144j 
for a complete listing of the formulae for all the baryon masses. 

There are no conceptual differences in the computation of the baryon masses in SU{3) vs. 
SU{2). The only differences are the addition of more operators in the chiral Lagrangian, due to 
the larger symmetry group of SU{3), which leads to kaon and eta loops, in addition to the pion 
loops of the previous section. We have already shown the Lagrangian necessary to compute the 
octet and decuplet baryon masses to NLO in Eqs. (|1.52|) and (|1.54j) . To compute the NNLO 
mass corrections, we must first construct the analogue of Eqs. ()2.11|) . (|2.13|) and (|2.14|) in SU{3). 
The treatment of the octet-decuplet mass splitting parameter is identical to the discussion in 
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Sec. 12.11 The reparameterization fixed Lagrangian is given by 




2Mq J \ 2Mo 
{D + F)tT (bvA^-^b] +{D + F)tT (b'-^^vAb] 

y Mo J y Mo J 

(D-F)tr ( B^^^BvA^ +{D-F)tv ^'J^BvA^ 

The other terms in the Lagrangian with undetermined LECs are given by 

C= -^tT(BA-AB) +^tr (BAuBA^) + -^tT(BBA-A) + ^tr (BB) tr M • ^) 

^t>j4 h"^ 
+ ^tiiBv Av- AB) + ^iT{BvABvA) + -r-tr CBBvAv-A) 
Ax ^ ^ A^ ^ ^ A^ ^ ' 

+ -^tr {BB) ti (v Av ■ A) + ^tv {B M+ M+ B) + ^tv{BM+BM+) 
Ax A^ 

— b^ — fo*^ — 

+ ^tv{BBM+M+) + ^tT{BM+B) tv {M+) + ^ti {B B M+) tv{M+) 
Ax A;,. A;^ 

5*^ _ b^^ — 

+ ^tr{B B) tr {M+ M+) + tr B) tr tr {M+) 

Ax ^x 
b^^ — 

+ tr M+) tr {B M+) . (2.35) 
Ax 

There are a few points worth making about the above Lagrangian. Firstly, the operators with 
coefficients b^_^ form an over-complete set. One can make use of Cayley-Hamilton identities 
for SU{2>) matrices to reduce these operators by one |157j . However, there is no consensus in 
the literature about which 7 of the 8 operators to keep, and so we hst them all for to form 
the complete basis. This over-complete basis also has a more obvious extension to the partially 
quenched operators we will determine in Sec. 12.31 Also, as in Sec. 12.11 one can show that the 
operators of the form ii {B S ■ AS ■ AB) are not linearly independent from the operators with 
coefficients b^_^ and ^1^4. 

The operators of the form {Bv ■ Av • AB) have not been kept explicit in the literature 
so far for calculations of the octet baryon masses to 0(m2) |2Z1IIM1IISZI- These operators 
have identical flavor structure to the corresponding tr{BA-AB) operators. However, these 
two sets of operators have different Lorentz structure, which gives rise to different finite 
contributions to the octet baryon self energy calculations. One can choose a renormalization 
scheme such that the contributions to the octet baryon masses from these different operators can 
not be distinguished. Therefore, with a suitable redefinition of the bf^_^ and b^j coefficients, 
the operators with coefficients b^^^ can be neglected in the baryon mass calculations, as their 
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contributions to the masses can be absorbed by the other operators to this order in the chiral 
expansion. 

However, the operators with coefficients bf_^ and b^^^ can be distinguished in vrA^ — > vrA^ 
scattering at tree level, for example. It is therefore useful to keep both types of operators explic- 
itly in the Lagrangian, allowing them to be distinguished in the octet baryon mass calculation. 
This also provides a consistent means to determine these LECs, as they are the same coefficients 
which appear in the Lagrangian used for the calculation of other observables like vrA^ ttN 
scattering. 

In principle, additional 1/Mq operators with the same chiral symmetry properties as those 
contained in Eq. 1)2. 35(1 can be generated. However, these I/Mq operators do not have their 
coefficients constrained by RI, therefore they can be absorbed by a re-definition of the Jj^''"^'^^ 
coefficients. For example, the Rarita-Schwinger field used to describe the decuplet baryons 
contains un-physical spin-^ degrees of freedom which can propagate when the decuplet baryons 
are off mass-shell. However, these off-shell degrees of freedom are suppressed by I/Mq and are 
implicitly included in the operators with coefficients b^_^^. 

There are a similar set of operators for the decuplet fields. The greater flavor symmetry of 
the decuplet fields, however, reduces the number of operators as compared to the octet baryons 
but the Lorentz structure allows for more types of operators. These operators are given by 



£ = ^ (A . A)/ + 'f- T^W ■ AfTi;, + 'f- {T'T,) tr (A • A) 

-i-vA . . j-vA , . . 4.vA 

+ 'j^T'^\v.Av. A)/ + 'j- T^v -A^^v. AfT,;, + 'f- T,) trivAvA) 

X XX 

+ f- {A^ . A)/ Ti;, + f- T^'A^^ . Al^ Tr,, + f- (r n) tr {A, A^ 
^^x ^^x ^^x 

^^x ^^x ^^x 

fM _ fM _ 

+ f - {T,M+T^^) tr {M+) + (T^T'^) tr {M+) tr {M+) . (2.36) 



All of the above LECs as well as those in Eq. ()2.35() are dimensionless. We now have the 
complete set of operators needed to compute the masses of the octet and decuplet baryons 
through 0{mg). As stated above, we shall only give the masses of the nucleons and deltas, and 
we shall restrict ourselves to the isospin limit, m„ = = rh. The quark mass expansion of the 
nucleon mass is given exactly as in Eq. 1)2. For simplicity, we treat rriK ~ rriTt in the power 
counting, although phenomenologically this may not be the ideal counting. The nucleon mass 
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is then 
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(2.37) 
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where the leading order masses of the octet and decuplet baryons are 



M^^^ = 2m {oM + Pm + 2crA/) + 2ms cfm , 



5 2^ , \ 

-"M + -PM + ^CFM I + 



m. -c 



-OtM + -PM + 2crM 

au + 213m + ^cfm) + ms{aM + '^ctm) , 

M^^^ = 2fh{-jM - 2a m) - 2ms^M , 

2 _ _ 
My.- = -{2m + ms){-fM - Scjm) , 



Mi" = 

f(l) 



(2.38) 



and the linear combination of LECs is 



i^A,vA _ '^iA,vA . rA.vA rA.vA . oi,A,vA 

"tt ~ 2 2 ~ ' 

j^A,vA _ 1 T^A,vA 1 ,A,vA , 1 r A,i;A , r A,«;A 

On - + + gf^s "T 04 ' 

^AvA^^^AM^^AM^^^AM^ 
^fnfn ^ ^ + 2hf + 46f , 

^msms ^ ^ ^ ^ _ (2.39) 



In Eq. (|2.37p . the first two lines are the LO and NLO nucleon mass corrections, and the rest 
of the expression is the NNLO mass contribution. Comparing to the nucleon mass in SU{2), 
Eq. (|2.22|) . we can see that Eq. 1)2. 37() is significantly more complicated. In the SU{3)v limit, in 
which the quark masses are all degenerate, this expression must reduce to a form more similar 
to Eq. (|2.22l) . but the group structure will still give rise to a splitting amongst the octet. The 
full listing of the octet baryon masses to 0{m^) can be found in Ref. |142j . 

In a similar fashion, we can compute the mass of the delta baryons in the isospin limit of 
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SU{3). The mass of the deltas are 

Ma = Mo + a + 27M fh - TaMi^rh + mg) 
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where the Unear combination of LECs is 



j.A,vA,A _ 3 A,i)A,A , ^ ,A,vA,A . nj.A,vA,A 
'■tt — 2 1 2 
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,A,vA,A _ ,A,vA,A , .,A,vA,A 
% ~ ^1 + ' 

^mm ^ ^ 2tf + Itf + 4tf , 

3 

I — ^3 + 2''4 + ^'^S ' 

.Af , .M 



tt+tf. (2.41) 



2.3 Baryon Masses in SU{4\2) and SU{6\3) PQxPT 

In Sec. 11.2.31 we showed the construction of the partially quenched heavy baryon Lagrangian 
to the order necessary to compute the baryon masses through NLO. In this section we will 
construct the next order partially quenched heavy baryon Lagrangian necessary to compute 
the masses to NNLO. Because of the symmetries of the partially quenched theories, at this 
order, the form of the PQ Lagrangian is independent of whether one is working with S'C/(4|2) 
or SU {6\3). Therefore, when we construct the operators, we will do it simultaneously for both 
theories and show how to match these Lagrangian operators onto the SU{2) and SU{3) theories 
respectively. 

The PQxPT Lagrangians in general have more linearly independent operators than their 
respective xPT Lagrangians (13(H 11311 1132j . Therefore, to determine the LECs of xPT from 
partially quenched LQCD simulations one must also match the operators order by order, 
which can be done by restricting the flavor indices of the partially quenched operators to the 
vector sub-space. The only other computational difference from the nucleon mass computation 
outlined in Sec. 12.11 is the addition of the hairpin interactions of the flavor neutral mesons 
which give rise to the double pole pieces of the propagators of these mesons, Eq. (|1.61|) . There 
are various ways to handle these double pole terms and for these mass computations, we use 
the method in Refs. |13n| I131j . To be specific, the two-point function of the flavor neutral 
propagators in momentum space, for three flavors of sea-quarks, is given by 

g fp2. ^ i^ah^a _ i_ (p^ - m|^. ) (p^ - m^, ) ^ 

— m^^ +ie 3 (p^ — 7712^ + ze)(p2 — rri^^ + ie){p^ — m\ + ie) ' 

where the mass m\ = ^{m'jj + 2m^^) and the masses m'jj are the masses of mesons composed 
of quark-antiquark flavor /. We can rewrite this propagator as 



^VaVb — ^a. Sab Pa + Hab [Pa, Pb^Px] , 



(2.43) 
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where the Hat operator is defined as 
Hab [A,B,C] 



A ; „ , , ID 
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(2.44) 



This then allows one to simply evaluate the integrals arising in the mass computation (and other 
observables) , as the double pole features are captured by the now mass dependent coefficients 
of the propagators in the Hab operator. A similar procedure can be used in the SC/(4|2) 
theory [HT] . 



2.3.1 Partially Quenched Heavy Baryon Lagrangian 

In Sec. 11.2.31 we showed the construction of the LO PQHBxPT Lagrangian. We now will 
construct the Lagrangian at the next order. The various operator structures that enter will 
involve the and A fields. We shall use a labeling convention very similar to that in 
Sees . 12 . l1 and 12 . 2l We shall then match the PQ Lagrangian onto the xPT Lagrangians to make 
the labeling conventions clear. The fixed coefficient PQ Lagrangian is 
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(2.45) 



where the flavor contractions are defined in Eq. 1)1.74(1 . The operators relevant to the mass of 
the spin-^ baryons with two insertions of the axial- vector pion fields are 
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and the operators with two insertions of the mass spurion are 



,M(PQ) ,M{PQ) 
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+ B^'^Bijk str {M+M+) + -^-r B'''\M+)i''Bnjk str {M+) 
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(2.47) 



We can then match these operators onto the xPT operators which gives the following relations 
amongst the LECs, 



2 1 2 2 33 12 ^ ^ 
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and for the matching of SU{Q\3) to SU{3) we find, 
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One must similarly match the operators relevant for the spin-| baryon masses. The relevant 
PQ Lagrangian is 
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Matching this Lagrangian onto the SU (2) Lagrangian, one can show that all the coefficients 
involving the two mass spurions have the same numerical values as the QCD coefficients, 
^M(4|2) _ ^*^(2)^ ^ g_ there is a one-to-one correspondence between these partially quenched 
operators and the SU{2) operators. This also holds for the matching of the SU{&\'i) opera- 
tors onto the SU{2>) operators. The partially quenched operators involving the axial-vector 
pion fields also match identically from SU{<6\'i) to 5C/(3), ^^'^^'^(^l^) _ ^^,-"-4,^(3)^ However, 
because of the symmetries of the SU (2) theory, there is one less operator, Eq. 12.141 such that 
the matching is 

A,vA,A{2) _ A,-uA,A(4|2) 
'■2 ~ ''2 

^Ai,A,i(2) ^ l^Ai,A,i(4|2) ^ ^A,vAMm _ (2.51) 



We now proceed to compute the mass of the nucleon in partially quenched HB^PT. 
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2.3.2 Nucleon Mass in S'?7(4|2) 

We find that the mass of the nucleon, in the isospin limit of the sea and valence sectors, with 
the valence quark mass, m„ and the sea quark mass, rrij, is given by 



M^l^^ = Mo- (2aM + Pm) rriu - 4aM mj 
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(2.52) 



where in the above expression we have used various replacements, 



L[m^,m^/,fj] = H 



m\ In 



m. 



rrij 



4 , / Y 
m V In — ^ 



(2.53) 



and we have used known linear combinations of the LECs, b^'^^, 6"", etc. |144j . If we extend 
this computation to SU {6\3), we will find a similar expression. The difference is that now there 
are 3 more quark flavors in the theory, the strange s-quark, the ghostly s, the s-quark, and the 
sea strange quark, the r-quark. Therefore this expression will simply be longer, as there will be 
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more mesons which can propagate in the loops, but the general structure of the mass will be 
identical to Eq. (|2.52|) . The full expression for the nucleon mass in SU (6|3) as well as the other 
members of the octet-baryons can be found in Ref. jl42j . Similarly, the mass of the deltas can 
be found for 5C/(4|2) in Ref. dS] and the fuh decuplet in Ref. 

We would like to conclude this chapter by mentioning that while this thesis was being writ- 
ten, the NPLQCD collaboration has computed both the deviation from the Gell-Mann-Okubo 
baryon mass relations [5^1 as well as the strong isospin breaking contributions to the nucleon- 
proton mass splitting j^Ol- Therefore, we anticipate that within one year, these expressions for 
the baryon masses to 0{ml) win be necessary for a more precise determination of these quanti- 
ties as well as the masses themselves. This is about 5 years earlier than the author anticipated 
when first computing these mass expressions. 
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Chapter 3 

INCLUDING BARYONS IN TWISTED MASS CHIRAL 
PERTURBATION THEORY 

This chapter is based upon the work in Ref. |158j . 
3.1 Introduction to Twisted Mass Lattice QCD 

Twisted mass lattice QCD (tmLQCD) |31j is an alternative fermion discretization technique for 
lattice QCD that has recently received considerable attention.^ It has the potential to match 
the attractive features of improved staggered fermions (efficient simulations |161j . absence of 
"exceptional configurations" 162 , 0{a) improvement at maximal twist |16,'-{j . operator mixing 
as in the continuum |311 11641 il65j ) while not sharing the disadvantage of needing to take roots 
of the determinant to remove unwanted degrees of freedom. Thus tmLQCD offers a promising 
and interesting new way to probe the properties and interactions of hadrons non-perturbatively 
from first principles. 

Due to the limitations in computational capabilities, the quark masses, niq, used in current 
simulations are still unphysically large. Thus extrapolations in ruq are necessary if physical 
predications are to be made from lattice calculations. This can be done in a systematic and 
model independent way through the use of chiral perturbation theory (xPT). Since xPT is 
derived in the continuum, it can be employed only after the continuum limit has been taken, 
where the lattice spacing, a, is taken to zero. However, when close to the continuum, it can be 
extended to lattice QCD at non-zero a, where discretization errors arising from the finite lattice 
spacing are systematically included in a joint expansion in a and rUq |166j . For tmLQCD with 
mass-degenerate quarks, the resulting "twisted mass chiral perturbation theory" (tm^PT) has 
been formulated previously |167[ 11681 1169] . building on earlier work for the untwisted Wilson 
theory [MH 11701 [T7T] . 

So far, tmxPT has only been applied to the mesonic (pionic) sector. There have been studies 
on pion masses and decay constants for niq ^ aAQ(-,j-, |167j . the phase structure of tmLQCD for 
niq ~ a 

^^QCD [miMlIinHlIinnillZll, and quantities involving pions that do not involve final 
state interactions up to next-to-leading order (NLO) in the power counting scheme where rriq ~ 
oAqqj-, jl75j . However, as pointed out in Ref. |175j . many of the pionic quantities considered are 
difficult to calculate in numerical simulations because they involve quark-disconnected diagrams. 
This motivated us to extend tm^PT to the baryon sector, enabling us to analytically study 
baryonic quantities that do not involve quark-disconnected diagrams. Numerical studies of the 
baryons in tmLQCD are already underway, and the first results from quenched simulations 
studying the nucleon and delta spectra have been obtained recently in Ref. |176j . 



^For recent reviews see Refs. |159llltiLl1 . 
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In the baryon sector, the extension of xPT to the lattice at finite lattice spacing to 0{a) |177j . 
and to 0{a^) has b een done for a theory with untwisted Wilson fermions. We extend that 

work here to include the effects of "twisting" , i.e. our starting underlying lattice theory is now 
tmLQCD. Specifically, we study the parity and flavor breaking effects due to twisting in the 
masses and mass splittings of nucleons and deltas in an SU{2) chiral effective theory. The mass 
splittings are of particular interest to us as they allow one to quantify the size of the parity-flavor 
breaking effects in tmLQCD; furthermore, they present less difficulties to numerical simulations 
than their counterparts in the mesonic sector, which involve quark-disconnected diagrams. 

We consider here tmLQCD with mass non-degenerate quarks |179j . which includes an addi- 
tional parameter, the mass splitting, eg. This allows us to consider the theory both in and away 
from the strong isospin limit. With simulations in the near future most likely able to access the 
region where rriq ~ oAqqj-,, the power counting scheme we will adopt is 

1 > ^ aAqcD ~ ~ (3.1) 

Aqcd Aqcd 

with denoting the small dimensionless expansion parameters. In the following, we will work 
to 0{e'^) in this power counting. 

The remainder of this chapter is organized as follows. In Sec. 13.2.11 and 13.2.21 we briefly 
review the definition of tmLQCD with mass non-degenerate quarks, and we show how the mass 
splitting can be included in the Symanzik Lagrangian and the 0{a) meson chiral Lagrangian. 
Higher order corrections from the meson Lagrangian are not needed for the baryon observables 
to the order we work. In Sec. 13.2.31 we extend the heavy baryon (HBxPT) to include the 
twisting effects to C(e^). In Sec. I3.3l we present the nucleon and delta masses in tm^PT in the 
strong isospin limit, including lattice discretization errors and the flavor and parity breaking 
induced by the twisted mass term. In Sec. 13.41 we extend the calculation to include isospin 
breaking effects and discuss the subtleties that arise. We conclude in Sec. 13.51 



3.2 Mass Non-Degenerate Twisted Mass Chiral Perturbation Theory 

In this section, we work out the extension of the baryon chiral Lagrangian to 0{a'^) given in 
Ref. |178j in tmLQCD. We start by briefly outlining the construction of the Symanzik La- 
grangian in the mass non-degenerate case, which follows the same procedures as those in the 
mass-degenerate theory with minimal modifications. 



3.2.1 The effective continuum quark level Lagrangian 

The fermionic part of the Euclidean lattice action of tmLQCD with two mass non-degenerate 
quarks is 

Sf= J^V5K^)[^2^7/.(V;^ + V^)-^J^V;;V^ + mo + i75nMo-r3eo]V'K^), (3-2) 

where we have written the action given in Ref. |179j for a general twist angle (not necessarily 
maximal), and in the so-called "twisted basis" |lti3j . The quark (flavor) doublets ipi and ipi 
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are the dimensionless bare lattice fields (with "/" standing for lattice and not indicating left- 
handed), and and V* are the usual covariant forward and backward dimensionless lattice 
derivatives, respectively. The matrices Tj are the usual Pauli matrices acting in the flavor space, 
with r3 the diagonal matrix. The bare normal mass, mo, the bare twisted mass, fiQ, and the 
bare mass splitting eo, are all dimensionless parameters; an implicit identity matrix in flavor 
space multiplies the bare mass parameter mo- The notation here is that both mo and eo are 
positive such that the upper component of the quark field is the lighter member of the flavor 
doublet with a positive bare mass. 

Note that in the mass-degenerate case, twisting can be done using any of the Tj, the choice 
of r3 is merely for convenience. Given the identity 

vr vr 

ex.p{-i-Tk) Ta exp{i-Tk) = ekabU + Sak Ta , (3.3) 

one can always rotate from a basis where the twist is implemented by Ta, a = 1, 2, to a basis 
where it is implemented by using the vector transformation 

- - vr 
ip ^ ip ex.p{i9Tk) , tp ^ ex.p{—i6Tk)ip , fe = l, 2, 3, 9 = ±—, (3.4) 

where the appropriate sign for 6 is determined by the index o. However, with used here to 
split the quark doublet so that the mass term is real and flavor-diagonal, it can not be used 
again for twisting if the fermionic determinant is to remain real.^ 

Following the program of Symanzik (180if il81) , and the same enumeration procedure detailed 
in Ref. jl69j . one can obtain the effective continuum Lagrangian at the quark level for mass 
non-degenerate quarks that describes the long distance physics of the underlying lattice theory. 
Its form is constrained by the symmetries of the lattice theory. To O^e"^) in our power counting, 
in which we treat oAqqj-, ~ m^ ~ Cq, we find that the Pauli term is again the only dimension 
five symmetry breaking operator just as in the mass-degenerate case jl69| (the details of this 
argument are provided in Appendix I A. 

Ces = J~-g + iPiP+ m + nsn/i - egTa)^ + bia-ip ia^^F^u tp + 0{a^) , (3.5) 

where Cg is the continuum gluon Lagrangian, m is the physical quark mass, defined in the usual 
way by 

m = Zmimo - mc)/a , (3.6) 

H is the physical twisted mass 

fi = Zf,fj.o/a = ZpVo/a , (3-7) 

^One way to see this is to note that the mass terms mo + i'fsrsfio — Tgeo can be written as {xo — 
rsyo) exp(iQ75) exp(i/375r3), where xo/yo = tan/3/tana. Thus, the twisted mass term can be transformed 
away leaving just the normal mass term and the mass sphtting term. However, since this involves an (7(1) 
axial transformation which is anomalous, an iaFF term is introduced into the action which we see now is 
complex (because of the factor of i). Thus, since the gauge action is real, this means that the fermionic action 
(before the transformation) must be complex, and so the fermionic determinant obtained from it must also be 
complex. This also implies that a theory, where both the twist and the mass splitting are implemented by ra, 
is a-dependent. 



59 



and Cq is the physical mass sphtting 

eg = Z^eo/a = Z~'^eo/a. (3.8) 

The factors Zp and Zs are matching factors for the non-singlet pseudoscalar and scalar densities 
respectively. Note that the lattice symmetries forbid additive renormalization to both and 
eo jl79j . The quantity rhc is the critical mass, aside from an 0{a) shift (see Ref. jl741 11751 11821 
I183j and discussion below). 

Anticipating the fact that the mesons contribute to the baryon masses only through loops, 
and so will be of 0{e^) or higher, we only need to have a meson chiral Lagrangian to 0{a) for the 
order we work; £eff as given in Eq. (|3.5j) is sufficient for its construction. To build the effective 
chiral Lagrangian for baryons to 0{e^) on the other hand, terms of 0{a'^) in Eq. ()3.5() are of the 
appropriate size to be included. However, except for the operator which breaks 0(4) rotation 
symmetry, a^'(/57^-D^-D^Z)^'i/') the 0(0^) operators do not break the continuum symmetries in a 
manner different than the terms explicitly shown in Eq. (|3.5j) . and thus their explicit form is 
not needed. The 0(4) breaking term will lead to operators in the baryon chiral Lagrangian at 
the order we work. However, it is invariant under twisting and thus contributes as those in the 
untwisted theory |178j . 

3.2.2 The SU{2) Meson Sector 

The low energy dynamics of the theory are described by a generalized chiral Lagrangian found 
by matching from the continuum effective Lagrangian 1)3. 5|) . As usual, the chiral Lagrangian is 
built from the SU{2) matrix-valued field S, which transforms under the chiral group SU{2)l x 
SU{2)ji as Eq. H1.17() The vacuum expectation value, Sq = (S), breaks the chiral symmetry 
spontaneously down to an SU{2) subgroup. The fluctuations around Sq correspond to the 
pseudo-Goldstone bosons (pions). 

From a standard spurion analysis, the chiral Lagrangian at O(e^) is (in Euclidean space^) 

= ^Tr(9^Sa^St) - ^Tr(xtS + St^) - ^Tr(its + sti) , (3.9) 

where / is the decay constant (normalized so that = 132 MeV). The quantities x and A are 
spurious for the quark masses and discretization errors respectively. At the end of the analysis 
they are set to the constant values 

X — > 2BQ{m + ifiTi - CqT^) =m + ifiTi - igT^ , A — > 2Woa = a , (3.10) 

where Bq ~ O(Aqcd) and Wq 0(Aqqj-,) are unknown dimensionful constants, and we have 
defined the quantities m, fi and a. 

As explained in Ref. |175j . since the Pauli term transforms exactly as the quark mass term, 
they can be combined by using the shifted spurion 

X=X + A, (3.11) 



^We will work in Euclidean space throughout this Chapter. 
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leaving the 0{e'^) chiral Lagrangian unchanged from its continuum form. This corresponds at 
the quark level to a redefinition of the untwisted component of the quark mass from m to 

m' = m + aWo/Bo. (3.12) 

This shift corresponds to an 0{a) correction to the critical mass, so that it becomes 

rric = Zmrhc/a — aWo/ Bq . (3.13) 

Since the O(e^) Lagrangian takes the continuum form, and the mass splitting term does not 
contribute at this order, the vacuum expectation value of S at this order is that which cancels 
out the twist in the shifted mass matrix, exactly as in the mass-degenerate case: 

(0|S|0)z.o = So = ^ + ^ eMi^on) , (3.14) 

where 



M' = V(^?r+a)2+^. (3.15) 

Note that M' is the leading order result for the pion mass-squared, i.e. = M' at 0{e'^). If 
we define the physical quark mass by 



nin 



^/m/^TJ^, (3.16) 



then it follows from (|3.14)) that 

coscjo = m' /rriq , sinc^o = n/mg . (3-17) 

Details of the non-perturbative determination of the twist angle and the critical mass can be 
found in HZi CZSl IMl CHa CHHl, and wih not be repeated here. 

At 0{e'^), the mass non-degenerate chiral Lagrangian for the pions retains the same form 
as that in the mass degenerate case |169( 1175) , because the mass splitting does not induce any 
additional symmetry breaking operators in C^g. The 0{e^) pion Lagrangian contains the usual 
Gasser-Leutwyler operators of 0{m^, mp^, p^) Eq. (|1.23)) . where m is a generic mass parameter 
that can be m, /i, or e^, as well as terms of 0(om, ap^, a^) associated with the discretization 
errors. Now as we stated earlier, since the pions will enter only through loops in typical 
calculations of baryon observables, keeping the pion masses to 0{e^) will lead to corrections of 
0{e^), which is beyond the order we work. As our concern is not in the meson sector, the 0{e'^) 
pion Lagrangian (|3.9|) is thus sufficient for our purpose in this work. 

3.2.3 The SU{2) Baryon Sector 

With the effective continuum theory and the relevant part of the eff'ective chiral theory de- 
scribing the pions in hand, we now include the nucleon and delta fields into tm^PT by using 
HBxPT llHlIlZlIinni, which we wih refer as the twisted mass HBxPT (tmHBxPT). The con- 
struction of the tmHBxPT Lagrangian parallels that in Sections 11.2.21 and 12.11 Recall that the 
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nucleoli is a doublet and the deltas a quartet under the SU{2)v transformations, with normal- 
izations given in Eqs. (|l.H4j) and (|1.44|1 . The free Lagrangian for the nucleons and deltas to 
0{e'^) consistent with the symmetries of the lattice theory is (in Euclidean space) 

=Niv DN -2aM NM^N - 2 aM NN tv{M+)-2aw NN tr(W+) 
+ {T^iv ■DT^) + A (T^T^) + 2 7m (T^M'^T^) 

- 2aM {T^T^) tr(7W*+-) - 2 aw {T^T^) tr(>V+) . (3.18) 
The "twisted mass" spurion field is defined by 



I- 



2 



^t„i*-e^ ± i{m'^)^i] , m*Q- = ^ , (3.19) 

with rn^Q being the "twisted" mass spurion for the baryons. The "Wilson" (discretization) 
spurion field is defined by 

W± = ^ {£,^wq£} ± ew^^e) , WQ = , (3.20) 

with wq being the Wilson spurion for the baryons. Note that we have made simplifications using 
the properties of SU{2) matrices when writing down Eq. 1)3. 18(1 . When setting the spurious to 
their constant values, W+ is proportional to the identity matrix in flavor space. Thus the 
operators NW+N and (T^W+T^), although allowed under the symmetries of tmLQCD, are 
not independent operators with respect to A^A^tr(yV+) and (T^T^) tr(yV+) respectively. This 
is also true of the nucleon and delta operators involving in the isospin limit (but not away 
from it). The independent operators we choose to write down are those with the simplest flavor 
contractions, and this will be the case henceforth whenever we make simplifications using the 
properties of SU{2). 

In Eq. 1)3. 18(1 . the four- vector, f^, is the Euclidean heavy baryon four- velocity, and our 
conventional here is that in Euclidean space, v ■ v = —1. The parameter. A, is the mass 
splitting between the nucleons and deltas which is independent of the quark masses and we 
treat A ~ m^r ~ following |281 1271 llL)5j . The dimensionless low energy constants (LECs), 
aM, CTM, 7Af, and aM have the same numerical values as in the usual untwisted two-flavor 
HBxPT. As was noted in Ref. jl75j . the shifting from x to x' = X + which corresponds 
to the shift of the physical mass m to m! at the quark level does not, in general, remove the 
discretization [A) term, and this is seen explicitly here with the presence of the discretization 
terms. 

At this order, the Lagrangian describing the interactions of the nucleons and deltas with the 
pions is still given by Eq. (|1.49j) . Note that the 0{e^) free Lagrangian (|3.18|) and interaction 
Lagrangian (|1.49() are the same as those given in Ref. jl78j when the twist is removed, i.e. 
when /i = 0. With non-vanishing twist, the mass operators carry a twisted component and 
the vacuum is "twisted" from the identity to point in the direction of the twist (the flavor 
Ti-direction here) jl75j . 

Following Ref. jl75j . we expand S about its lattice vacuum expectation value, deflning the 
physical lattice pion fields and the physical lattice ^ fields by 

S = TT,phT , ^ = TiphU{iph) , T = exp(zu;on/2) , 

Sp/, = exp(i^/2 TT • r//) , T,UeSU{2), (3.21) 



62 



If we make the following chiral transformation for which the fields transform as in Eqs. (|1.17|1 . 
(fT^Kl) . (fT3?71) . (fnpi and fTI^ . and use the particular SU{2) matrices L = i?t = Tt, we have 
in the transformed effective chiral Lagrangian 



ph, 



(3.22) 



and 



tw 




Note that since L = R^i = T^i e SU{2), and ^ = L^^phU = WS,phR, 
S = L%hR = e = (L^^phU) ■ {U^iphR) = L^^phR = 



Kph 



-'ph 



(3.23) 



(3.24) 



We see that the ^ field is now ^p/j, the field associated with the physical pions, and the twist 
is transferred from the twisted mass (A1?^) term to the "twisted Wilson" (W^) term, making 
the mass term in the HB^PT now the same as that in the untwisted theory. The new mass 
spurion, mg, and the "twisted Wilson" spurion, u^g", now take constant values 



..tw 



a exp(— icJoTi) 



(3.25) 



We will call this the "physical pion basis" since this is the basis where the pions are physical as 
defined by the twisted lattice action jl75j . and we will work in this basis from now on, unless 
otherwise specified.^ A technical point we note here is that, in the isospin limit where twisting 
can be implemented by any of the three Pauli matrices, say Tk, the physical pion basis can be 
found following the same recipe detailed above but with ri in T replaced by Tk throughout. 

Rotating to the physical pion basis where the ^ field is now the physical ^p/i field in all field 
quantities, the form of the interaction Lagrangian remains unchanged as given in (|1.49j) . while 
the O(e^) free heavy baryon chiral Lagrangian 1)3. 18(1 changes to 



NivDN -2aM NM+N - 2aM NNtv{M+) 
+ {T^iv ■DT^) + A (T^T^) + 2 74/ (T^M+T^) 



2 aw NNtr{W: 



tw\ 



2Wm {T^T^) ti{M+) - 2 aw {T^T^) tr(W. 



tw \ 



(3.26) 



Note that is also proportional to the identity matrix in flavor space when set to its constant 
value. Thus if we build the free chiral Lagrangian directly in the physical pion basis, the same 



''As detailed in Ref. |175) . the twist angle that one determines non-perturbatively in practice, call it u), will 
differ from ujo by 0(a). This will give rise to a relative 0(a) contribution to the pion terms. But since the 
pions come into baryon calculations only through loops, the correction will be of higher order than we work. 
Thus to the accuracy we work, we may use either u) or a;o- 
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simplifications due to SU{2) we used in writing down Eq. (|3.18|) apply. Note that at this point, 
one can not yet tell whether the nucleon (A^) and the delta (T^) fields are physical. This has to 
be determined by the theory itself. We will return to this point when calculating the nucleon 
and delta masses below. 



At ©(e^), there are contributions from 0{avn) and 0{a^) operators. The enumeration of 
the operators is similar to that set out in Ref. jl78j . except now the Wilson spurion field carries 
a twisted component. The operators appearing in the 0{e^) chiral La grangian will involve 
two insertions off the following: A4±, , and the axial current A^. Note that since parity 
combined with flavor is conserved in tmLQCD, any one insertion of Ai- or must be 

accompanied by another insertion of or W^. Now operators with two insertions Mj^ or 
A^, which contribute to baryon masses at tree and one-loop level respectively, have the same 
form as those in the untwisted theory (and so give the same contribution). These have been 
written down in ,144^ and will not be repeated here. Operators with an insertion of either a 
combination v ■ A and (which have the same form as in the untwisted theory), or a 
combination of v ■ A and Wi"", will also not contribute to the baryon masses at 0{e^). 

At (D{am), there are two independent operators contributing to the masses in the nucleon 
sector 



and two independent operators contributing to the masses in the delta sector 



C -J- 



tf^+(r^A^+T^)tr(W: 



tw 



{T^T^)ii{Wf)iT{M 



(3.27) 



(3.28) 



where we remind the reader A-^ = 47r/. Note that there are no operators involving the commu- 
tator, [A^+,yV^], because it is identically zero. There are also operators involving ^Wi"" 
at 0{am), but these again do not contribute to the baryon masses at the order we work. 

At 0{o?), there are operators that do not break the chiral symmetry arising from the 
bilinear operators and four-quark operators (see e.g. Ref. |171j for a complete listing) in the 
mathcalO{a'^) part of £eff- These give rise to the tmHBxPT operators 



bNN + t{T^T^ 



(3.29) 



There are also chiral symmetry preserving but 0(4) rotation symmetry breaking operators 
which arise from the bilinear operator of the form a^tpj^D^Dfj^D^tp in the 0{a'^) part of £efr- 
These give rise to the tmHBxPT operators 



(3.30) 



Note that these chiral symmetry preserving operators are clearly not affected by twisting (the 
0(4) symmetry breaking operator at the quark level from which they arise involve only deriva- 
tives with no flavor structure, and {7^, 75} = 0), and so they have the same form and contribute 
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to the baryon masses in the same way as in the untwisted theory. The chiral symmetry breaking 
operators at 0{a'^) are those with two insertions of the Wilson spurion fields. For the nucleons, 
there are two such independent operators 



1 



^1 NNtr{Wf) tr(Wf ) + n^" NN ti{W^_^W^_!") 
and for the deltas, there are three such independent operators 



(3.31) 



C -J- 



tf+ (T^T^) tr(W7) tr(Wf ) +tr' {T^T^) triWL^W] 



ituii A jtw\ 



(3.32) 



In the isospin limit where the mass splitting vanishes (e^ 0), more simplifications occur 
in the 0{e'^) chiral Lag rangian. The nucleon operators in Eq. (|3.27j) with coefficients ^ 
are the same up to a numerical factor, and the same holds for delta operators in 



and n. 



Eq. (|3.28j) with coefficients 1^^^ and ^2*^*^^, and for operators in Eq. (|3.32j) with coefficient 



,WM. 



ij^ and 



Note that in the untwisted limit, the 0{e^) chiral Lagrangian reduces to that given in 
Ref. |178j . In particular, with the twist set to zero, operators with two insertions of Wi"" will 
not contribute to the nucleon or the delta mass until 0{a'^m) ~ O(e^), but for non- vanishing 
twist, they contribute at 0{a?). 



3.3 Nucleon and Delta Masses in the Isospin Limit 

In this work, we are concerned with corrections to the masses of the nucleons and the deltas 
due to the effect of the twisted mass parameter. We will therefore only give expressions for 
the mass corrections arising from the effects of lattice discretization and twisting in tmLQCD. 
A calculation of the nucleon and delta masses in the continuum in infinite volume to 0{m'^) 
can be found in Ref. |144j . The mass corrections due to finite lattice spacing to O(a^) in the 
untwisted theory with Wilson quarks can be found in Ref. [178j . and the leading finite volume 
modifications to the nucleon mass can be found in Ref. j52j . 

In this section, we present the results of nucleon and delta masses calculated in tmHBxPT, 
in the isospin limit, where the quark doublet is mass-degenerate, and the twist is implemented 
by T^. As we discussed in Sec. 13.2.11 in the isospin limit, the content of tmLQCD is the same 
regardless of which Pauli matrix is used to implement the twist - the action for one choice is 
related to another by a flavor-vector rotation. This must also hold true of the effective chiral 
theory that arises from tmLQCD. Indeed, the heavy baryon Lagrangian constructed in Sec. 13.2.31 
with Ti-twisting can be rotated into that with rs-twisting by making a vector transformation, 
which is given hy L = R = V = exp(i^r2). 

3.3.1 Nucleon Masses in the Isospin Limit 

In the continuum, the mass of the nucleons in infinite volume HB^PT with two flavor-degenerate 
quarks are organized as an expansion in powers of the quark mass, which can be written as 
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Figure 3.1: Diagrams depicting mass contributions to the nucleons through next-to- next-to- 
leading order (NNLO) in tmHBxPT in the physical pion basis. The solid, thick solid (red) 
and dashed lines denote nucleons, deltas and pions respectively. The solid triangle denotes 
an insertion of the twisted Wilson operator as given in (|H.26|) . The solid squares denote the 
coupling of the baryons to the axial current whose form is given in (|1.49|) . The clear triangle 
denotes a tree level insertion of the operators given in Eqs. (|3.27j) and (|3.29j) - (|3.31|) . 



Eq. 1)2. Here, we are interested in the corrections to this formula due to the effects of lattice 
discretization and twisting arising from tmLQCD. We denote these lattice corrections to the 
nucleon mass at 0(e^") ~ 0{m'^) ~ 0(a") (factors of Aqcd needed to make the dimensions 

correct are implicit here) as 5M^\ and the nucleon mass in tmHBxPT is now written as 

M*,7 = Mo - (M« + <5mJ;) ) + . . . (3.33) 

The leading correction in tm^PT comes in at tree level, arising from the twisted Wilson 
nucleon operator in the free heavy baryon Lagrangian (|3.26|) and is depicted in Fig. 13.11 

6M^^]{to) = 4:awak\cos{Lo) , (3.34) 

where to the accuracy we work, lo can either be or the twist angle non-perturbatively 
determined. Note that this correction is the same for both the proton and the neutron. At 
leading order, the nucleon mass is automatically 0{a) improved, as 5M^-* vanishes at maximal 
twist, uj = -k/2. At zero twist, a; = 0, it reduces to that in the untwisted theory jl77l 1178] , 

The next contribution to the nucleon mass comes from the leading pion loop diagrams shown 
in Fig. 12.11 However, at the order we work, the form of the ©(m^/^) nucleon mass contribution 
is unchanged from the continuum, and given in Eq. 1)2. 9|) . 

(2) 

The corrections to Mj^ come from both the tree level and the one-loop diagrams, as shown 
in Fig. 13.11 The twisted Wilson operator in the free Lagrangian ()3.26)) gives rise to a tadpole 
diagram, which produces a contribution of 0{am). The leading Wilson spurious also contribute 
to 0{am) when inserted inside the pion-nucleon loops, and are partly cancelled by wavefunc- 
tion corrections. The tree level contributions come from the operators given in Eqs. 1)3.27)) 
and ()3.29)) - ()3.31)) . Just as in the untwisted continuum theory, these act both as the higher 
dimensional operators and as counter terms that renormalize divergences from the lower order 
loop contributions. For instance, coefficients n^^^ and n^^^ are renormalized to absorb di- 
vergences from the tadpole and one-loop contributions mentioned above. These coefficients are 
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taken to be the renormalized coefficients (finite) in the mass calculations, and to contain the 

(2) 

counter terms needed in our renormalization scheme. The corrections to Mj^ read 

-5M^^\uj) = 12 aw amllog (^^^ cos(w) 

+ 165AAf i'^w - <yw)a [j{m^,A,^j.) + ml] cos{uj) 

WM+ I „ WM+\ . .X 2a3 , I, \ 

ii + 2 rig I aA^ mq cos^w) — a (o + b^) 

+ a^Kl (2nf- sm^{uj) - 4nf+ cos2(t^)) . (3.35) 

Note that the 0{£^) corrections are again the same for both the proton and the neutron. At 
maximal twist, the 0{e'^) corrections are given by 

-5M^^]{uj = 7r/2) = a^A^ (^2nf- - 6 - 6^) , (3.36) 

while at zero twist, these reduces to the corrections given in Ref. jl78j . We see that the nucleon 
masses are also automatically 0{a) improved at 0{e'^). 

To the order we work, the expressions for the nucleon mass corrections in tmHBxPT given 
in Eqs. (|3.34|) and (|3.35|) . together with the untwisted continuum HB^PT expressions for the 
nucleon masses, provide the functional form for the dependence of the nucleon masses on the 
twist and angle, a;, and the quark mass, mq, which can be used to fit the lattice data. 

3.3.2 Delta Masses in the Isospin Limit 

We are again interested in the corrections to the delta masses arising from the twisted mass 
parameters, such that the delta mass expansion is written 

= Mo + A + {M^^ + 6Mff ) + ... (3.37) 

The leading mass correction arises at tree level from the twisted Wilson delta operator given in 
Eq. (|3.26j) , which we depict in Fig. 13.21 

5M^^ (uj) = -4 aw a K\ cos(u;) . (3.38) 

Just as for the nucleons, this does not split the delta masses and vanishes at maximal twist. 

The 0{e^) delta mass contributions are similarly given as for the nucleons. They do not 
cause any splitting between the deltas, and receive no discretization corrections, and are shown 
in Fig. Vn\ 

At ©(e^), contributions due to the effects of twisting arise from similar diagrams as in the 
nucleon case, and are shown in Fig. 13.21 A splitting in the delta masses first arises at this order, 
which comes from the operator with coefficient ~ given in (|3.32|) . The mass corrections read 

5M^^{uj) = 12 aw amllog cos(w) +4^^^ (aw - aw) a J{mT,,-/S., ^) cos{u:) 

+ 2 (tf ^+ + 2 *'^+) aA^ mq cos(w) + a^A^ [t + t^) 

+ o^A^ (At^+ cos2(cj) -2tY~ sin2(u;) + t^- 5t, sm^{u)) , (3.39) 
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Figure 3.2: Diagrams depicting mass contributions to the deltas at NNLO in tmHBxPT in 
the physical pion basis. The solid, thick solid (red) and dashed lines denote nucleons, deltas 
and pions respectively. The solid triangle denotes an insertion of the discretization operator in 
Eq. (|3.26j) . The solid squares denote the coupling of the baryons to the axial current whose 
form is given in Eq. (|1.49|) . The clear triangle denotes a tree level insertion of the operators 
given in Eqs. (IT^ - (FOn|l and (jSSl- 



(b) 



where 



-1 forri = A++,A- 



3 



(2) 

Note the appearance of the mass splitting, (Jt^, in 6Mj./ . We see from above that starting at 
O(e^), the delta multiplet is split into two mass-degenerate pairs, with one pair containing 
and A~, and the other. A"*" and A*^. At maximal twist, 6M^^ becomes 

5M^^ (to = tt/2)= a^Al + tf" + "^t,) , (3.41) 

while at zero twist, it reduces to that given in Ref. |178j . Just as in the nucleon case, the delta 
masses to 0{e^) are also automatically 0{a) improved. 

As is the case with the nucleons, to the order we work, the expressions for the delta mass 
corrections in tmHBxPT given in Eqs. (|3.38|) and (|3.39|) . together with the untwisted continuum 
HB^PT expressions for the delta masses, provide the functional form for the dependence of the 
delta masses on the twist angle, uj, and the quark mass, rriq, which can be used to fit the lattice 
data. 



3.3.3 Mass Splittings 

Having derived the expressions for the nucleon and delta masses in tmHBxPT to order 0{£^) 
in the isospin limit, we now focus on the mass splittings between the nucleons and between the 

(n) (n) 

deltas. The mass contributions in the continuum, Mj^/ and Mj._ , clearly do not give rise to 
mass splittings for the nucleons and deltas, since they are calculated with degenerate quarks. 
Therefore, any mass splitting can only come from the mass corrections arising from tmLQCD. 

From the results of Sec. 13. 3. H and Sec. 13.3.31 we find that to 0{e'^), the protons and neutrons 
remain degenerate, while the delta multiplet splits into two degenerate pairs, with A^^ and 
A~ in one pair, and A"^ and A" in the other. The splitting between the degenerate pairs in the 
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delta multiplet is given by 

- M^++.- = ^ t^'a^Al sm\co) = ^ tf'a^A^ 4 • (3-42) 
S ^ S ^ 

We reiterate here that the 0{a) uncertainty inherent in the definition of the twist angle results 
in a correction to M^+,o — M^++,- of O(a^) ~ 0{e^), which is of higher order than we work. 
Hence to the accuracy we work, we may use ujq or any other non-perturbatively determined 
twist angle for lo above. 

Just as the case of the pion mass splitting worked out in Ref. jl75j , this delta splitting must 
vanish quadratically in afj, = orrig sin(aj)^ on general grounds, since the masses do not violate 
parity. One would therefore expect, naively, the splitting to be O(a^mg) ~ 0{e^). But as our 
results show, there is, in fact, a mass dependence in the denominator such that the effect is 
0{e'^). Suppose we take = 2 GeV, then we would find a mass splitting 

Ma+,0 - Ma++,- ~ 1.1 tf" GeV. (3.43) 

Using naive dimensional analysis, we expect ~ 0(1), giving a mass splitting of the delta 
pairs on the order of 1 GeV. This is rather large, and in fact unexpected. In fact, a recent study 
of the quenched tmLQCD spectrum found the splitting on the order of 50 to 100 MeV 176 . 
If however, we were to use Aqcd as the dimensionful parameter in Eqs. 1)3. 20() and (|3.32|) . and 
guess some reasonable size of Aqcd ~ 500 MeV, then we would have found the mass splitting 
to be 0.04*2 GeV, and then naive dimensional analysis estimates for would place this 
LEG as 0(1), in much better agreement with the twisted mass lattice data |158j . We can not 
make a direct comparison with the quenched data, as our computation is for full tmLQCD. 
Nevertheless, the quenched spectrum is expected to be within 20% of the full QGD result, and 
so it would be nice to determine this LEG, t^' ■ 

Now, the degeneracies we found for the nucleons and the delta multiplet above hold not only 
at O(e^), but in fact they hold to all orders in tm^PT. This can be understood by considering 
the lattice Wilson-Dirac operator associated with the action of tmLQGD given in Eq. 1)3. 2|) in 
the isospin limit with rs-twisting 

DwD = ^ 5Z + ^m) - ^ K^f^ + "^0 + ^TsrsMo , (3.44) 

M At 

which has the self-adjointness property |179j 

n75 DwD 75n = D\y ^ . (3.45) 

It follows then that the propagator for the upper and lower component of the quark doublet, 
tpi{x), call them Su{x,y) and Sd{x,y) respectively, satisfy the relations 

l5Su{x,y)j5 = Sl{y,x) , -y^ Sd{x,y) j5 = Sl{y,x) . (3.46) 



^Note: here n is the twisted mass parameter, not to be confused with the renormahzation scale. 
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This means that any baryon two-point correlator which is invariant under the interchange of 
the quark states in the quark doublet combined with hermitian conjugation, leading to the 
degeneracies mentioned above. An argument of this type has been given in Ref. |176j . 

The same can also be shown in a chiral Lagrangian treatment, as must be the case. Now one 
of the symmetries of tmLQCD with two flavor-degenerate quarks and ra-twisting is the pseudo- 
parity transformation, Vp, where ordinary parity is combined with a flavor exchange jlB.'Sj . 

''Uo{x) Uo{xp) , xp = (-x,t) 
Uk{x)^ul{xp), fc = l, 2, 3 ^^^^^ 

iplix) ^ iTl'Jolplixp) 

ipiix) -ii>i{xp)loTi 



where are the lattice gauge fields. At the level of HB^PT, this is manifested as the invariance 
of the chiral Lagrangian under the transformations 

< — >n{xp) , A++'-(x)< — ^A-'++(xp) , A+'°(x) < — > A°'+(xp) , 

Ofc(x) ®fliTiOfc(xp)ri , (3.48) 

where for an operator in the chiral Lagrangian, Ok is any operator matrix that contracts with 
the flavor indices of the the nucleon (N) or delta (T^) flelds in the operator. If the or the 
fields contained in an operator have a total of 2F flavor indices, ^^^^Ok is the tensor product 
of F operator matrices which contract with the F distinct pairs of these flavor indices. The 
degeneracies in the nucleons and the delta multiplets discussed above would then follow if all 
the operators in the chiral Lagrangian that contribute to the baryon masses have a structure 
that satisfies the condition 

^LiOkix) = ^k=inOk{xp)n , (3.49) 

Consider first the case for the nucleons. Since the nucleon fields are vectors in fiavor space, 
we can take F = 1 without loss of generality (the nucleon fields can only couple to one operator 
matrix). Since the chiral Lagrangian is built with just M.±, , Af_i, and V^, the operator 
matrix Ok can only be constructed from combinations of these fields. We need not consider 
combinations involving just and W^, since the fiavor structure of both is trivial, i.e. pro- 
portional to the identity. We also need not consider mass contributions arising from pion loops, 
because they must have the same flavor structure as the tree level local counterterms used to 
cancel the divergences in these loops. Therefore, we do not have to consider operators involving 

and V/i, which give rise to mass contributions only through pion- nucleon interactions. This 
leaves us with only combinations involving and as possible candidates to break the 
degeneracy in the nucleons. As was discussed in Sec. 13.2.31 because of the parity-flavor sym- 
metry of tmLQCD, Ok can not contain just a single A4„ or , but must always have an 
even number from the set {A4_,Wi"'}. Now any such combination would indeed have a pure 
tree level part, however, it is also proportional to the identity in flavor space. Thus there is no 
operator matrix. Ok, that one can construct which violates the condition Ok{x) = TiOk{xp)Ti. 

The arguments for the case of the deltas runs similar to that for the nucleons. For the same 
reason given in the nucleon case, we need not consider operator structures that involve 
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Afj,, and V^. We need only consider operator structures involving an even number from 
the set {M-,W^J^}. For the deltas, F can be three since each delta field has three flavor indices. 
But since two of the in (^l^j^O^ must come from the set {Ai^,yV^} to satisfy the parity- 
flavor symmetry of tmLQCD, we can take F to be at most two without loss of generality. Now 
each of M.- and W!!" has a tree level part that is proportional to rs, thus, under Vp, Oi ® O2 
where Ok can be either or , satisfies the symmetry condition, Eq. (|H.49|) . Therefore, 
one can not construct operators for the deltas that break the degeneracy between the pairs in 
the delta multiplet. 

3.4 Nucleon and Delta Masses Away from the Isospin Limit 

In this section, we present results for mass corrections due to twisting away from the isospin 
limit, where the quarks are now mass non-degenerate. To the order we work, the corrections 
due to the mass splitting come in only at tree level. For clarity, we will only point out the 
change arising from the quark mass splitting; we will not repeat the discussion on the nucleon 
and delta masses that are the same both in and away from the isospin limit. 

3.4-1 The Flavor- Diagonal Basis for the Mass Matrix at C'(e^) 

The natural choice for splitting the quark doublet is to use the real and flavor-diagonal Pauli 
matrix, rs, since the quark states one uses on the lattice correspond to the quarks in QCD in 
the continuum limit. But as was discussed in Sec. 13.2.11 above, twisting can not be implemented 
with r3 in this case (the fermionic determinant would be complex otherwise), and so ti is used 
instead. 

Since the twist is implemented by a flavor nondiagonal Pauli matrix away from the isospin 
limit, flavor mixings are induced for non-zero twist: the quark states in tmLQCD are now linear 
combinations of the physical quarks of continuum QCD. At the level of the chiral effective theory, 
this manifests itself in that the hadronic states described by the tm^PT Lagrangian are linear 
combinations of the continuum QCD hadronic states we observe, viz. the pions, nucleons, 
deltas, etc. 

If the effects from twisting are perturbative as compared to the isospin breaking effects, 
the hadronic states described by tmxPT will be "perturbatively close" to their corresponding 
continuum QCD states, i.e. the difference between them is small compared to the scales in the 
theory (see Add endix I A . 21 for an explicit demonstration). In this case, we can still extract QCD 
observables directly from tm^PT, as the corrections will be perturbative in the small expansion 
parameter. However, if the twisting effects are on the same order as the isospin breaking effects 
so that the flavor mixings are large, these corrections will not be perturbative.^ Nevertheless, 

qualitative guide to the size of the flavor mixings can be found in the ratios of two-point correlation 
functions. Deflne the ratio of QCD delta states by 

^ (A' A^) + (A^' A') 
" ~ (A« A«) + (Aj Aj) ■ 

Flavor mixing should be small if the off diagonal elements of Rij are small. To determine the size of the flavor 
mixings quantitatively, one has to look at the splitting in the delta multiplet. We will discuss further in the 
text below. 
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one can still extract information for the QCD observables: One can still measure the masses 
of these tmLQCD hadronic states in lattice simulations, and one can fit these to the analytic 
expressions for these masses calculated in tm^PT to extract the values of the LECs. The 
LECs associated with the continuum xPT contributions have the same numerical values as in 
tmxPT. Therefore, if one determines these from tmLQCD simulations, one knows the masses 
of the QCD hadronic states. 

At the order we work, flavor mixings are manifested in the appearance of flavor non- 
conserving pion-baryon vertices in the Feynman rules of tmHBxPT, and in that the baryon 
mass matrix is not flavor-diagonal. Since we work in the physical pion basis where the twist 
is carried by the Wilson spurion (now flavor non-diagonal) instead of the mass spurion (now 
flavor-diagonal) , flavor mixings can only arise from operators with one or more insertions of the 
Wilson spurion field. Because of this, the flavor non-conserving pion-baryon vertices and the 
non-diagonal terms in the mass matrix must be proportional to a, the lattice spacing, and so 
must vanish in the continuum limit where the effects of the twist are fake and can be removed 
by a suitable chiral change of variables |3 H I162j .^ 

For the nucleons, flavor mixings induce only unphysical flavor non-conserving pion-nucleon 
vertices which vanish in the continuum limit; the nucleon matrix is still flavor diagonal at the 
order we work. In fact, this is true to all orders in tmHBxPT. The reason is the same as that 
given in Sec. 13.3.31 We need only consider the tree level part of the possible operator structures 
that one can construct from the spurion fields in tmUB^PT. Now the only spurion field that 
has a tree level part with non-diagonal flavor structure is W^, and as we discussed above, it 
must be paired either with another or with Ai-, which renders the flavor structure of the 
tree level part of the combination trivial. Thus we may take the basis of nucleons used in the 
tmHBxPT Lagrangian as the physical nucleon basis. 

For the deltas, not only are there flavor non-conserving pion-delta vertices, at the order 
we work, the delta mass matrix is already flavor nondiagonal at tree level. This happens for 
the deltas because the tensor nature of the field allows more freedom in the way the flavor 
structure of the delta operator can be constructed. Thus, in order to have only physical tree 
level mass terms for the deltas, we must change to a basis where the delta mass matrix is 
diagonal, which can now only be done order by order. 

When diagonalizing the delta mass matrix, we need, in principle, to diagonalize the mass 
matrix that contains all the mass contributions from both tree and loop level to the order that 
one works. But we find the difference between diagonalizing the delta mass matrix including 
both tree and loop level contributions at the order we work, and diagonalizing that with only the 
tree level mass contributions, give rise to corrections only to the loop level mass contributions, 
which are higher order than we work. Thus, we will diagonalize the delta mass matrix containing 
just the tree level mass terms in our calculation for the delta masses. 

To the order we work, if the tree level mass is given by 

v^Mava, va = {A++ AO A+ A-) , WA = (A++ A" A+ A")^ , (3.50) 



^This shows again the convenience of the pion physical basis, where all the effects of symmetry breaking 
in the lattice theory are parametrized and contained in the Wilson spurion fields, which then vanish as the 
symmetries are restored in the continuum limit. 
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where and are vectors of the delta basis states used in the tniHBxPT Lagrangian, and 
Ma is the tree level mass matrix, the physical delta basis is given by 

Va = -VA, v'^ = v^-S, (3.51) 

where S is the matrix of eigenvectors of Ma such that 

S-Ma-S-^=V, (3.52) 

with D the corresponding diagonal eigenvalue matrix. This implies that 

v^Mava = (v'^ ■ S-^) -{S-V- S-^) ■ [S ■ v'a) = v'-^Vv'a ■ (3.53) 

The full details of the diagonalization are provided in Appendix lA. 21 In the following sections, 
we will work in this basis for calculating the delta masses. 

The Feynman rules in the new basis are obtained from the same tmHE^PT Lagrangian 
given above in Sec. 13.2.31 but with each of the delta flavor states now rewritten in terms of the 
new delta flavor states given by the defining relations Eq. (|3.51|) . Note that changing to the 
new delta basis induces new unphysical flavor non-conserving vertices in the delta interaction 
terms given in (|1.49j) . because in terms of the new basis states, flavors are mixed. However, 
these flavor mixing components are proportional to the off-diagonal elements of 5, which are 
proportional to the lattice spacing as well as the twist angle (see Appendix IA.2|) . Thus they 
vanish in the limit of vanishing twist or lattice spacing, and so the unphysical vertices arising 
from them also vanish in these limits. 

We note and reiterate here that in the isospin limit, this order by order mass matrix diago- 
nalization is unnecessary as one can always rotate to a basis where the twist is flavor-diagonal 
from the outset, and issues of flavor nonconserving vertices and non-diagonal mass matrices due 
to flavor mixings do not arise. ^ 



3.4.2 The Nucleon Masses 

Away from the isospin limit, the first change caused by the mass splitting occurs in the contin- 
uum mass contribution since the quark masses 

ruu = niq - eg , nid = nig + eg , niq ,eq> 0, (3.54) 

are no longer equal. 

At the order we work, the only other change due to the mass splitting appears at 0{e'^) in 

the contribution from the 0{am) nucleon operator with coefficient n}^^+ given in Eq. (|3.27jl . 

(2) 

In the isospin limit, its contribution to 6M^ (uj) is proportional to rUg, but away from the 
isospin limit, it becomes 

2nJ^*^"'"aA^mgCos(a') — > 2 nj^^^a A^ mj cos(a') , (3.55) 

where rrii is defined in Eq. (|2.4j) . The corrections to the nucleon masses from the effects of lattice 
discretization and twisting are otherwise the same as those given in Eqs. (|3.34|) and 1)3. 35(1 . Note 
that the nucleon masses are automatically 0{a) improved, just as in the isospin limit. 

*In fact, as is shown in Appendix IA.2I if one insists on remaining in the basis where the twist in flavor non- 
diagonal, one would find that the unphysical terms arising from flavor mixings do not vanish in the continuum 
limit. 
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3.4-3 The Delta Masses 

Away from the isospin limit (e^ 7^ 0), we calculate the delta mass and mass corrections in the 
basis where the delta mass matrix is diagonal to the order we work. This diagonalization is 
worked out in Appendix IA.21 where we have obtained general expressions for the new delta 
basis that are valid in the range from = to eg ~ oAqqj-,. Here, we present the case where 
Eg > and Eg ~ ruq ~ aAq^j-, ^ o^Aqqj-,, which is a regime that simulations in the near future 
can probe. To the order we work, we may take the new delta basis states in this regime to be 

Ti ^ |Ai) =Ci 
^ IA3) = C3 
T2 ^ IA2) = C2 
^ IA4) = C4 

where T/ = Aj denote the deltas in the new basis, Ci are normalization factors, and 

A = 2eg (7Af + ^+aA^ cos(u;)) , B = tf^a^A^ sin^iuj) . (3.57) 

Note that A ~ 0{e'^) and B ~ ©(e^) in our power counting, so B /A ^ C'(e^) and the effects of 
the flavor mixings are perturbative. 

The masses of these states are comprised of the continuum expressions given in Chapter |21 
Ref . jl44j , and corrections due to the effects of discretization and twisting. Note the continuum 
expressions for the delta masses here are necessarily changed from that in the isospin limit 
because 7^ m^. The mass corrections due to the effects of lattice discretization and twisting, 
come in at tree level; the loop contributions remain unchanged from that in the isospin limit. 
The tree level mass contributions to the order we work have been worked out in Eq. ()A.11|1 . in 
the process of diagonalization. We list here the full delta mass corrections to ©(e^), which we 
denote by SMt^, to the mass of the delta state denoted by T': 

— ^ J cos(u;) 
+ ^9an i^w - crw) a J {niT, , -A, /i) cos{uj) 

+ 2tf*'^+aA^^cos(w) +4tf^+aA^ mqcos{oj) + kUt + t^) 
3 

+ a^Al (4tf+ cos^iuj) - 2tf- sm^{uj) + tf" sm^{uj)) , (3.58) 
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(3.56) 



74 



for T = 1 ... 4 and 



Sruu 

ruu + 2md 
3md 



for i = 1 
for i = 2 
for 2 = 3' 
for 2 = 4 



5^, 



for z 
for i 



1,4 
2, 3 



(3.59) 



Note that (JM-p. (w) as given in Eq. I|3.58() . is the same as the sum of 6m!^^ and 6M!j^,' as given 
in Eqs. (|3.38|) and (|3.39|) respectively, but with the changes 
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The full expressions for the delta masses can be obtained when the continuum contributions 
are included. To the order we work, one can obtain the complete mass expression for delta 
denoted by T' to 0{£^) in tmHE^PT by adding the mass corrections, (JM-p. (a;), to the continuum 
mass of the delta denoted by T, whose expression can be found in Section f2. 1.21 

We stress here that one can not take the isospin limit from any of the expressions give above 
in this subsection. They have been derived for 7^ and with the assumption that the twisting 
effects are much smaller than the isospin breaking effects. One must use the general formulae 
given in Eqs. ()A.lfl|l and ()A.11|) when considering cases where these conditions are not true. 

Observe that away from the isospin limit, the delta masses are also automatically 0{a) 
improved at maximal twist [u = 7r/2), as all terms proportional to a in 5M are proportional to 
cos(a;) as well. Hence, to the order we work, the contributions due to the isospin breaking are 
the same as that in the continuum at maximal twist. 



3.5 Summary 

In this chapter we have studied the mass spectrum of the nucleons and the deltas in tmLQCD 
with mass non-degenerate quarks using effective field theory methods. We have extended heavy 
baryon chiral perturbation theory for SU (2) to include the effects of the twisted mass, and we 
have done so to 0{e'^) in our power counting, which includes operators of 0{am'^ ^ap^ ,0^). Using 
the resulting tmHE^PT, we have calculated the nucleon and the delta masses to 0(6^), and we 
found them to be automatically 0{a) improved as expected from the properties of tmLQCD. 

Because of the twisting, the vacuum is no longer aligned with the identity in flavor space, 
which has non-trivial effects on the physical excitations (pions) of the theory. Also, depending 
on whether the quarks are mass degenerate or not, the way twisting is implemented determines 
what the physical baryon states are in the theory. We have highlighted these subtleties when 
doing calculations in tmHE^PT. 

In order for the pions in the theory to be physical, we have to make a particular (non- 
anomalous) chiral change of variables to undo the twisting effects. This requires the knowledge 
of the twisting angle, but once that is determined, the physical pion basis can be determined d 
priori. However, whether or not the nucleons and deltas are physical must still be determined 
from the theory. In the isospin limit, both the nucleon and the delta mass matrices are diagonal, 
and so the nucleon and delta states contained in the A*" and fields are physical. However, 
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away from the isospin limit, only the nucleon mass matrix remains diagonal. Thus, the N field 
can still be regarded as physical, but the physical deltas are now linear combinations of the 
flavor states contained in the field. This can be understood from the fact the at the quark 
level, the physical QCD states, the u and d quarks, are eigenstates of T3 but not of ri. So 
only in the isospin limit, where the twist can always be implemented by the flavor-diagonal 
Pauli matrix, T3, are the states contained in the quark doublet physical quarks. Away from 
the isospin limit, the twist can not be implemented by rs anymore, and the eigenstates of the 
Hamiltonian of the theory are composed of linear combinations of the u and d quarks. 

The physical states in tmxPT are in general a mixture of those in the untwisted xPT. The 
size of the mixture is determined by the relative sizes of the discretization effects, which are 
0(0^), and the isospin splitting effects, which are C(eq). In this chapter, we have given general 
expressions for the nucleon and delta masses with respect to this mixing of states that are valid 
in the range from = to ~ oAqcd- 

The quantities which provided the motivation for this work and turned out to be most 
interesting are the mass splittings between the nucleons and between the deltas. We found 
that in the isospin limit, the nucleon masses do not split to any order in tm^PT, while the 
delta multiplet splits into two degenerate pairs. This can be understood from the symmetries 
of tmLQCD at the quark level, and as we have shown, also at the level of tm^PT. The mass 
splitting between the multiplets, M^+,o — M^++,-, first arises from a tree level contribution at 
0{a^), and it gives an indication of the size of the flavor breaking in tmLQCD. This splitting 
in the delta multiplet will be easier to calculate in lattice simulations than the corresponding 
quantity m^^ — m^^ ^ the meson sector jl75j . since it involves no quark disconnected diagrams. 

Twisted mass HB^PT can also be extended to partially quenched theories (such extension 
of tmxPT for pions has recently been done jl85j ). This will be useful as unquenched twisted 
mass lattice QCD simulations are underway jl86j . 
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Chapter 4 

ELECTROMAGNETIC POLARIZABILITIES OF THE NUCLEON FOR 

LATTICE QCD 

The work in this chapter is based upon Ref. jl87j 
4.1 Introduction 

Compton scattering at low energies is an invaluable tool with which to study the electromagnetic 
structure of hadrons. At very low photon energies, the Compton amplitude is dominated by 
point-like photon scattering from the total charge and magnetic moment of the target hadrons. 
As the frequency increases, contributions beyond point-like scattering enter and one begins to 
resolve the hadronic response to an applied electromagnetic field. For unpolarized scattering 
on spin one-half objects, the first structure dependent contributions in this energy expansion 
of the amplitude are the electric polarisability, a, and the magnetic polar isability, /?. These 
quantities reflect the ability of the hadron's components to align or anti-align themselves in 
response to an applied electric or magnetic field. For the proton and neutron, the positivity 
of the accepted experimental values of these polarizabilities {op = 12.0 ± 0.6, Pp = 1.9 =F 0.6, 
an = 12.5 lb 1.7 and betan = 2.7^ 1.8 in units of 10~^ fm^ |188j ^ indicates that both nucle- 
ons are diamagnetic objects. Recent experimental advances |189| I188j have also allowed the 
extraction of certain combinations of target polarisation-dependent observables in Compton 
scattering. These involve the spin polarizabilities [HI], conventionally labeled 71-74, and they 
have consequently been investigated in numerous theoretical and further experimental studies. 
Although the classical interpretation of spin-dependent Compton scattering is less clear, the 
spin polarizabilities encode additional fundamental properties of the nucleon. Compton scat- 
tering observables, however, are not limited to these six parameters. Higher order quasi-static 
properties of the nucleon appear from further terms in the energy expansion of the amplitude. 
These higher-order polarizabilities jSH]) as well as generalized polarizabilities _35^ (which arise 
in the singly (doubly) virtual Compton scattering process, j*X — > j^*^X) allow for an even 
finer resolution of the electromagnetic structure of hadrons at low energies. 

While experimentally one is hoping to open further windows through which to view hadronic 
electromagnetic structure, theoretically one ultimately hopes to understand how hadronic po- 
larizabilities arise from the basic electromagnetic interaction of the photon with quarks that are 
bound to form the hadrons. The electric and magnetic polarizabilities should naively scale with 
the volume of the hadron. However, this expectation overestimates the observed polarizabilities 
by four orders of magnitude, indicating that the nucleon's constituents are strongly coupled. 
Lattice techniques provide a method to investigate the non-perturbative structure of hadrons 
directly from QCD. In particular, the various hadron polarizabilities can be computed. Com- 
parison of these results with experimental determinations would provide stringent tests of the 
lattice method's ability to reproduce the structure of physical hadronic states; for the individual 
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spin polarizabilities that have not been measured, the lattice approach may be the only way 
to determine them. On the lattice, direct calculations of the required hadronic current-current 
correlators are difficult and so far have not been attempted. However significant progress has 
been made |19()| [TDl, 192 in extracting the electric and magnetic polarizabilities by performing 
quenched lattice calculations in constant background electric and magnetic fields respectively 
and studying the quadratic shift in the hadron mass that is induced (essentially an application 
of the Feynman-Hellman theorem). These studies have investigated the electric polarizabilities 
of various neutral hadrons (in particular, the uncharged vector mesons and uncharged octet and 
decuplet baryons), and the magnetic polarizabilities of the baryon octet and decuplet, as well as 
those of the non-singlet pseudo-scalar and vector mesons. As we shall discuss below, generaliza- 
tions of these methods using non-constant fields allow the extraction of the spin polarizabilities 
from spin-dependent correlation functions and also allow the electric polarizabilities to be de- 
termined for charged hadrons. More generally, higher-order polarizabilities and generalized 
polarizabilities are accessible using this technique. 

As with all current lattice results, these calculations have a number of limitations and so are 
not physical predictions that can be directly compared to experiment. For the foreseeable future, 
lattice QCD calculations will necessarily use quark masses that are larger than those in nature 
because of limitations in the available algorithms and computational power. Additionally, the 
volumes and lattice-spacings used in these calculations will always be finite and non-vanishing, 
respectively. For sufficiently small masses and large volumes, the effects of these approximations 
can be investigated systematically using the effective field theory of the low energy dynamics 
of QCD, chiral perturbation theory (xPT) 11001 EOl EH- ^ In this chapt er we shall perform an 
analysis of the nucleon electromagnetic and spin polarizabilities at next-to-leading order (NLO) 
in the chiral expansion. We do so to discuss the infrared effects of the quark masses and finite 
volume in two-flavour QCD and its quenched and partially-quenched analogues (QQCD and 
PQQCD). The polarizabilities of the hadrons are particularly interesting in this regard since 
they are very sensitive to infrared physics and their quark mass and finite volume dependence is 
considerably stronger than that expected for hadron masses and magnetic moments. This should 
be physically evident given that the polarizabilities scale with the volume. In essence, chiral 
perturbation theory provides a model independent analysis of the modification of the nucleon's 
pion cloud in a finite volume. When the charged pion cloud is infiuenced by to the periodic 
boundary conditions imposed on the lattice, the nucleon's response to external electromagnetic 
fields is altered compared to that at infinite volume, and in most cases the effects are dramatic. 
A particularly striking oddity that we find in this analysis is a modification of the Thomson cross 
section at finite volume. This can be explained through the physics of chiral loop corrections 
to point-like hadron structure. 

If future lattice QCD simulations are to provide physical predictions for the electromagnetic 
and spin polarizabilities, careful attention must be paid to both the chiral and infinite volume 
extrapolations. To illustrate this point, we present our results at representative values of the 
quark mass, finding significant effects. We also use our quenched chiral perturbation theory 
results to assess the volume dependence of the quenched data at the lightest pion masses used in 

^The effects of the lattice discretization are short distance in nature, and while some of them can be analyzed 
in an extension of the effective field theory described here |166l 117(71 11711 11771 11781 1193| . Here we will assume 
that a continuum extrapolation has been performed. 
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Refs. |19H 1192] . While the quenched theory contains unphysical low energy constants (LECs) 
and the convergence of the chiral expansion is questionable at these pion masses, we can still 
provide an estimate of the volume dependence of quenched data for the nucleon polarizabilities 
using our results. Such an estimate is achievable because the corresponding polarizabilities in 
the unquenched theory do not depend on phenomenologically undetermined LECs at the order 
of the chiral expansion to which we work. At the lightest quark masses used in the existing 
quenched lattice simulations, m.,^ ~ 0.5 GeV, we find strong sensitivity to the lattice volume 
(as large as 10%). The effects will only increase as the pion mass is brought closer to that in 
nature. Clearly careful chiral and volume extrapolations of polarizabilities are mandated to 
connect lattice calculations to real world QCD. 

To begin our investigation of nucleon polarizabilities in lattice QCD, we discuss in Sec. 14.21 
the kinematics of Compton scattering and define the electromagnetic and spin polarizabilities 
that are the primary focus of this work. In Sec. 14.31 we perform a general analysis of the 
external field method pertaining to all electromagnetic and spin polarizabilities. We discuss 
how suitable background fields can be used in lattice QCD simulations to determine the spin 
polarizabilities and, more generally, generalized polarizabilities (though we limit our discussion 
of these in the present chapter). Following this we introduce the low-energy effective theories of 
QCD (xPT), quenched QCD (Q^PT) and partially-quenched QCD (PQxPT). These effective 
theories provide the model independent input necessary for calculating the quark mass and lat- 
tice volume dependence of polarizabilities. We focus primarily on PQ^PT in Sec. 14.41 discussing 
the relation to xPT where relevant, and relegating the peculiarities of QxPT to Appendix lB.il 
Our results for the dependence of the nucleon polarizabilities on quark masses and the lattice 
volume are presented in Sec. 14.51 We provide detailed plots relevant for full QCD simulations of 
polarizabilities showing the dependence on quark masses and lattice volumes. We also estimate 
the quenched QCD volume dependences of the polarizabilities at a pion mass typical of existing 
quenched lattice data. A glossary of finite volume functions required to evaluate the polariz- 
abilities in a periodic box appears in Appendix IB. 21 Lastly, Sec. 14.61 consists of a concluding 
discussion of our results. 

4.2 Nucleon Compton Scattering and Electromagnetic Polarizabilities 

The real Compton scattering amplitude describing the elastic scattering of a photon on a spin- 
half target such as the proton or neutron can be parameterized as 



T^N = Ai{uj,9)e' ■ e + A2{uj,9)e' ■ ke- k' + iA:i{uj,e)a ■ {e' x e) 

+ i Aiiio, 9) a-{k' X k) e -e + i ^5(w, 9) a ■ \{e x A:) e • k' -{ex k') e ■ k 





where we have chosen to work in the Breit frame of the system^ and the incoming and outgoing 
photons have momenta k = {uj,k = u)k) and k' = {uj,k' = cvk'), and polarisation vectors e 



^The Breit frame is not actually a reference frame, but a convenient choice of coordinates in momentum 
space for a fixed momentum transfer, such that the initial hadron momentum is given by pi = p — | and the 
final momentum by p/ — p + ^. Here, p is the average of the initial and final momentum of the hadron and 
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and e', respectively. The Ai{uj,9), i = 1 . . . 6, are scalar functions of the photon energy and 
scattering angle, cos 6 = k ■ k' . It is convenient to work in Coulomb gauge throughout where 
^0 = ^0 = (the physical amplitudes are gauge invariant). 

The functions, Ai, determining the Compton amplitude can be separated into a number of 
pieces. The Born terms, or tree level graphs, describe the interaction of the photon with a point- 
like target with mass, M^r, charge, eZ (where e > 0), and magnetic moment, /i. These terms 
reproduce the Thomson-limit (the zero frequency limit) and quadratic frequency pieces |194j of 
unpolarized scattering as well as the Low-Gell-Mann-Goldberger (LGMG) low energy theorems 
[1951 1196j for spin-dependent scattering.^ The remaining parts of the amplitude describe the 
structural response of the target. Expanding the amplitude for small photon energies relative 
to the target mass and keeping terms to 0{u}^) one can write 



Ai{u;,e 
A2{u,e 
A^{uj,e 
^4(0;, 6* 



+ 



Mn 4M3 



+ cos e) - Z^) (1 - cos 9) 47r(a ^ P cos 6)u! + 0{u! ) , 
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(/i^ - Z^)uj^ cos 9 - AnPio'' + 0{uj'' 



2Ml 



(Z(2/i - Z)- ^? cos 9) + 47ra;^(7i - (72 + 274) cos 9) + 0\ 



e uo 
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/i^ + 47ra;^74 + 0{iJ' 



'2MI, 



(4.2) 



describing the target structure in terms of the electric, magnetic and four spin polarizabilities, a, 
/?, and 7i-4, respectively. In the conventions above, the spin polarizabilities receive contributions 
from the anomalous decay vr'^ — > 77 (shown in Fig. 14.21 below) . This contribution varies rapidly 
with energy and is omitted from the polarizabilities in some conventions. Higher order terms in 
the energy expansion can be parameterized in terms of higher-order polarizabilities |36j . The 
more general process of virtual (and doubly-virtual) Compton scattering at low energies can 
similarly be described in terms of generalized polarizabilities |35| . We will focus in the six 
polarizabilities defined above. 



p — m " which highlights the fact that this is not a reference frame. The total momentum transfer to 
the hadron is q. 

''The LGMG theorems relate for example, the two-photon-nucleon interactions which arise from two successive 
photon insertions to the square of the single photon insertion. These are the one-nucleon reducible graphs, 
or the one-particle reducible graphs in quantum field theory. In Figure 14.11 the third diagram is related to 
the square of the second diagram, while the fourth diagram is a new contribution, which gives rise to the 
nucleon polarizabilities. To be completely explicit, at the second order in the nucleon-photon interactions, one 
contribution comes from the square of the nucleon magnetic dipole moment, which is related by the LGMG 
theorems to the square of the single magnetic dipole interaction, while the spin-dependent polarizabilities are a 
qualitatively new feature of the nucleon electromagnetic structure, and given by the interaction of the photons 
with the pion-cloud of the nucleon. 
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The goal of this chapter is to determine the quark mass and volume dependence of the 
polarizabilities defined above to allow accurate extraction of their physical values from lattice 
calculations. Before we do this we shall discuss how these lattice calculations may be imple- 
mented. 

4.3 Compton Scattering and Polarizabilities on the Lattice 

Lattice QCD provides a way to study the polarizabilities of hadrons from first principles. There 
are two ways to do this. The method most reminiscent of the experimental situation is to study 
the (Euclidean space) four point Green function defining the Compton scattering tensor directly 
(the photon fields are amputated). By measuring the large Euclidean time behaviour of this 
correlator, the hadron matrix elements of the two vector currents can be extracted. In princi- 
ple, by calculating particular Lorentz components of the Compton tensor with various different 
source and sink spin states, all six electromagnetic and spin polarizabilities and their higher 
order and generalized analogues can be extracted. However, this is a complicated task, requir- 
ing the evaluation a large number of quark propagator contractions resulting from quark-line 
disconnected diagrams which are statistically difficult to determine. At present this approach 
is too demanding for the available computational resources and has not been attempted. 

The second method is based on measuring the response of hadronic states to fixed external 
fields. A number of exploratory quenched QCD studies have been performed in this approach. 
The pioneering calculations of Refs. |197|ll98|ll99|[78l[79ll2n0j attempted to measure the nucleon 
axial couplings, magnetic moments and electric dipole moments by measuring the linear shift 
in the hadron energy as a function of an applied external weak or electromagnetic field. As 
discussed in the Introduction, various groups |19U1 11911 1192j have also used this approach to 
extract electric and magnetic polarizabilities in quenched QCD by measuring a quadratic shift 
in the hadron energy in external electric and magnetic fields. The method is not limited to 
electroweak external fields and can be used to extract many matrix elements such as those that 
determine the moments of parton distributions and the total quark contribution to the spin of 
the proton j2Ulj . Here we focus on the electromagnetic case. 

The Euclidean space (x4 = r) effective action describing the gauge and parity invariant 
interactions of a non-relativistic spin-half hadron of mass M and charge q with a classical U(l) 
gauge field, A'^{x,t), is 




(4.3) 



for the Lagrangian 



27ri (-7Sii?i<7 ■ExE + jm^m^^ -H xH + -fM,E2(^' E'^ W + jE.Ah^^'H'^ E^) *(^, t)+. . . , 



(4.4) 



81 



where E = —-^A{x,t) — V^4(x, r) and H = V x A(x,t) are the corresponding electric and 
magnetic fields, X = -^X denotes the Euclidean time derivative, X^^ = ^{d^X^ + d^X^), 
and the ellipsis denotes terms involving higher dimensional operators. By calculating one- 
and two-photon processes with this effective Lagrangian, it is clear that the constants that 
appear in Eq. (|4.4|) are indeed the relevant magnetic moment and electromagnetic and multipole 
polarizabilities (202] [these are simply related to the polarizabilities defined in the previous 
section as: 'JEiEi = -(71+73), 7A/iAfi = 74, lEiAh = 73 and-fMiE2 = 72 + 74]- The Schrodinger 
equation corresponding to Eq. (|4.4j) determines the energy of the particle in an external U(l) 
field in terms of the charge, magnetic moment, and polarizabilities. Higher order terms in 
Eq. (|4.4|) (which contain in part the higher order polarizabilities [301) t)e neglected for 
sufficiently weak external fields. For a magnetic field, the minimally coupled terms generate 
towers of Landau levels and for a constant electric field the same terms accelerate charged 
particles. 

Lattice calculations of the energy of a hadron in an external U(l) field are straight-forward. 
One measures the behaviour of the usual two-point correlator on an ensemble of gauge configu- 
rations generated in the presence of the external field. This changes the Boltzmann weight used 
in selecting the field configurations from det m] exp [—Sg] to det lp+ i + m exp [— 

where Ip is the SU(3) gauge covariant derivative, Q is the quark electromagnetic charge operator, 
and Sg is the usual SU(3) gauge action. Since calculations are required at a number of different 
values of the field strength in order to correctly identify shifts in energy from the external field, 
this is a relatively demanding computational task (although it is at least conceptually simpler 
than studying the four-point function). In general one must worry about the positivity of the 
fermionic determinant calculated in the presence of a background field, however for weak fields, 
positivity is preserved. The exploratory studies of Refs. [WllTIMITIMITIMITRllTmiMTllTimiTn^ 
used quenched QCD in which the gluon configurations do not feel the presence of the U(l) field 
as the quark determinant is absent. In this case, the external field can be applied after the 
gauge configurations had been generated and is simply implemented by multiplying the SU(3) 
gauge links of each configuration by link variables corresponding to the fixed external field: 
{Ua{x)} — > {[/a (x) exp[i e a A'^]}, where a is the lattice spacing. These studies are interesting 
in that they provide a proof of the method, however the values of the polarizabilities extracted 
have no connection to those measured in experiment. 

It is clear from Eq. ()4.4() that all six polarizabilities can be extracted using suitable space 
and time varying background fields if the shift of the hadron energy at second order in the 
strength of the field can be determined. One can also see this because the Compton tensor 
appears explicitly as the second-order connected term in the expansion of hadronic two-point 
correlation function in weak background fields |2Ulj . Previous studies |19Ul 11911 I192j have 
employed constant electric and magnetic fields to determine the corresponding polarizabilities 
in quenched QCD. Here we perform a more general analysis to show how the spin polarizabilities 
and the electric polarizabilities of charged particles can be obtained. 

In order to determine the polarizabilities, we consider lattice calculations of the two-point 
correlation function 

CsAp,r;A) = J d^xe'P^-^{0\xs{x,T)xU0,0)\0)A, (4.5) 
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where Xs{x,t) is an interpolating field with the quantum numbers of the hadron under con- 
sideration (we will focus on the nucleons) with z component of spin, s, and the correlator is 
evaluated on the ensemble of gauge configurations generated with the external field, A^. 

For uncharged hadrons at rest in constant electric and magnetic fields, it is simple to show 
that this correlator falls off exponentially at large times with an energy given by the appropriate 
terms in Eq. (|4.4j) owing to the constancy of the effective Hamiltonian. However for space-time 
varying fields, charged particles or states of non-zero p, a more general analysis is needed. This 
is most easily formulated using the effective field theory (EFT) defined by Eq. (|4.4|) . For weak 
external fields (such that higher order terms in Eq. 1)4. 4() can be safely neglected), the small 
p and large r dependence of this QCD correlation function is reproduced by the equivalent 
correlator calculated in the effective theory corresponding to the Lagrangian, Eq. (|4.4jl . That 
is 

CsAp,t;A) = jd^xe'P-'^^^^jv¥Vi'^s{x,T)^l,{0,0)expi-S,s[A]), (4.6) 

where Z^slA] = f V'^^'V^ exp {-Ses[A]). Since the right-hand side of Eq. (|4.6() is completely 
determined in terms of the charge, magnetic moment and polarizabilities that we seek to ex- 
tract, fitting lattice calculations of Cssi{p,t; A) in a given external field to the effective field 
theory expression will enable us to determine the appropriate polarizabilities. In the above 
equation we have assumed that the ground state hadron dominates the correlator at the rele- 
vant times. For weak fields this will be the case. However one can consider additional terms in 
the effective Lagrangian that describe the low excitations of the hadron spectrum that have the 
same quantum numbers as the hadron under study. This will lead to additive terms in Eq. (|4.6|) 
that depend on the mass, magnetic moment and polarizabilities of the excited hadron instead 
of those of the ground state. With precise lattice data, the properties of these excited states 
can also be determined. 

In many simple cases such as constant or plane-wave external fields, the EFT version of 
Css'ip, t; A) can be determined analytically in the infinite volume, continuum limit j2()Hj . How- 
ever in finite lattice spacing and at finite volume, calculating Css'{p,t; A) in the EFT becomes 
more complicated. In order to determine the EFT correlator, we must invert the matrix /C 
defined by 

Si-^tt[A] = Y.^Y1 ^^)^ss'[x, T^, y, Ty; A]^s{y, Ty) , (4.7) 

y,Ty s,s' 

where Siattl^] is a discretization of the EFT action in which derivatives are replaced by finite 
differences (the time derivative is given by a forward difference as we can ignore anti-particles). 
/C has dimension ANf where Ni is the number of lattice sites. For the most general space-time 
varying external field, this must be inverted numerically; given a set of lattice results for the 
correlator, Eq. (|4.6|) is repeatedly evaluated for varying values of the polarizabilities until a 
good description of the lattice data is obtained. 

For weak fields such that |^^(x, r)p <C ^qc^ for all x and r, a perturbative expansion of 
/C~^ in powers of the field can be used. This corresponds to the series of diagrams in Fig. 14. 11 To 
extract all six polarizabilities using such an analysis, we need to consider a number of different 
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fields; lattice calculations of the correlators in Eq. (|4.5|) using 
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for a number of different choices for the strength parameters, aj and bi (with |aj|, <C Aq^^j-,), 
are sufficient to determine the full set of polarizabilities.^ By measuring correlators for different 
spin configurations (including those that flip spin), we can reduce the number of fields required 
to extract the polarizabilities. 

As an example, the behaviour of the correlator in the field A'^^-^ (x) (which corresponds to a 
constant electric field in the xi direction) is given by 



Ss,s' exp {-^ [ai {q'r' + UMna) - Siqrp,] } e-^^^e" ^irbl' + 0(a?) 



Ip1->o 



6s,s' exp 



6M 

-{M + 2iraal)T 



1^3 



6M 



+ 0(4). 



(4.9) 



In this case, the perturbative series has been resummed exactly in the continuum, infinite volume 
limit and the higher order corrections come from terms omitted in Eq. (|4.4jl . For electrically 
neutral particles, the exponential fall-off of this correlator determines the polarisability a once 
the mass M has been measured in the zero-field case. When a charged particle is placed in such 
a field it undergoes continuous acceleration in the xi direction (this is described by the term 
in the exponent). However at times small compared to the correlator essentially falls off 

exponentially. Matching the behaviour of Eq. 1)4.9(1 to lattice data for a charged hadron will 
again enable us to determine the electric polarisability, a. 



*These fields correspond to real E and H fields in Minkowski-space for real- valued and bi. Since periodic 
spatial boundary conditions are envisaged for the link variables, there are quantization conditions that must be 
satisfied by the Ui 198 199 . For example, 5^02 = for each of the quark charges qi. The more complicated 
fields in Eq. 114.811 require two parameters to satisfy these conditions. 
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As a second analytic example, we consider one of the multipole polarizabilities. In the 
presence of the field A^gj(x), which corresponds to a more complicated electric field E^q-^{x) = 

{^T,i^,0), we find that 



7^ — 7^ A — T = exp 

<-u(p,^;^(6)) 



— 0.6 be 1 El El T 
a 



+ ..., (4.10) 



independent of p and the ellipsis denotes terms cubic in the field that have been neglected in 
Eq. H4.4() . Whilst the individual correlators, C-^^ and C^, have relatively complicated time- 
momentum behaviour involving q and a as well as "fEiEiy tliis becomes very simple in the ratio 
and "fEiEi can be determined cleanly. 

Analogous results can be derived for the other fields in Eq. (|4.8j) . however to take into 
account the finite lattice spacing and periodic finite volume nature of the underlying lattice 
simulations to which the EFT description is matched, the correlator is most easily calculated 
by inverting the matrix /C numerically. This also allows for more general choices of fields. If we 
seek to extract higher order polarizabilities, the Lagrangian in Eq. (|4.4j) must be extended to 
include higher dimension operators j36j . At this order, relativistic corrections and three-photon 
couplings also need to be included. Correlation functions similar to those in Eq. 1)4. 5|) involving 
two different external momenta will allow us to also extract the generalized polarizabilities [^^ . 

4.4 Heavy Baryon Lagrangian with Electromagnetic Interactions 

To calculate the quark mass and volume dependence of the nucleon polarizabilities, we use heavy 
baryon chiral perturbation theory (HBxPT) as was first constructed in Refs. [271 1281 [1081 l2U4j . 
In current lattice calculations, valence and sea quarks are often treated differently, with sea 
quarks either absent (quenched QCD) or having different masses than the valence quarks 
(partially-quenched QCD).^ The extensions of HB^PT to quenched HBxPT |138| 11391 and par- 
tially quenched HB^PT |130l I131j to accommodate these modifications are also well established 
and have been used to calculate many baryon properties. In this section, we will primarily focus 
on the two flavour partially-quenched theory and briefly introduce the relevant details following 
the conventions set out in Ref. |131j . Since QCD is a special limit of the partially-quenched 
theory, our discussion also encompasses two flavour xPT. Additional complications in quenched 
xPT are relegated to Appendix IB. II 

4-4- 1 Pseudo-Goldstone mesons 

We consider a partially-quenched theory of valence (n, d), sea {j, I) and ghost (■u, d) quarks 
with masses contained in the matrix 

mq = diag{mu,md,mj,mi,mu,m^) , (4.11) 

^At finite lattice spacing, different actions can even be used for the different quark sectors (e.g., staggered sea 
quarks and domain wall valence quarks). As was shown in Refs. |177ll2Ui3lll93| . the lattice spacing corrections 
to baryon electromagnetic properties are expected to be small, as they can not enter at tree level, and for 
current simulations with aAq^j:, ~ m,, they generally enter at leading loop order through valence-sea meson 
masses. In our work we assume a continuum extrapolation has been performed. 



85 



where m- j = rriu^d such that the path-integral determinants arising from the valence and ghost 
quark sectors exactly cancel. The corresponding low-energy meson dynamics are described by 
the PQxPT Lagrangian, Eq. (|1.57|) . with the addition of minimally coupled electromagnetism 
(the U(l) gauge field is again denoted by A'^ and its field strength tensor F'^'^ = d^A'^ — d'^A'^) 
to the theory through the chiral, and U(l) gauge covariant derivative 

Vf" = ^^' + [V^, ] , (4.12) 

with the vector current 

' ^ C{d^- ieQA'') + ^Qf^ _ ieQA^") i 1 , (4.13) 



= - ^ 
2 



depending on the quark charge matrix, Q. In coupling electromagnetism to this theory, we 
must specify how the quark charges are extended to the partially-quenched theory. We choose: 

Q = diag{qu,qd,qj,qi,qu,qd) , (4.14) 

though other arrangements are possible. However, one must set qj + qi ^ in order to retain 
sensitivity to the full set of LECs that appear in two flavour xPT j2U61 12U7] . In addition to 
the Lagrangian, Eq. (|1.57|) . the anomalous couplings of the Wess-Zumino-Witten Lagrangian 
|2()8| f2fl9) will also contribute to the spin polarizabilities. These terms are described below. 



4.4- 2 Baryons 

The baryon Lagrangian is given in Eq. (|1.75|) . again with the addition of minimally coupled 
electromagnetism, 



A" 



(4.15) 



As with the mesons, at leading order the photon is minimally coupled to the baryons with 
fixed coefficients. At the next order in the expansion there are a number of new electromagnetic 
gauge invariant operators which contribute to the Compton amplitude and the polarizabilities. 
Here, we display the relevant terms at this order, 



C 



(1) 

B 



2Mn 



+^^str[Q5+] (i3[5^,5n B) 



(4.16) 



where ^«,/3,7 are magnetic moment coefficients |1391 113 J] , /ir is the coefficient of the Ml tran- 
sition 70-44 operator HTUl HOSl and 



(4.17) 
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The partially quenched magnetic moment coefficients are related to the isoscalar and isovector 
magnetic coefficients, /io and fii, in standard two flavour xPT as 

^0 = ^ (/"a + + , /Ul = ^ (2/Ua - , (4.18) 

where the xPT Lagrangian describing the magnetic moments of the nucleons [the proton and 
neutron magnetic moments are fip^n = ^(/^o i fii)] is given by 

>C = ^F^.(MoiV[5^5^7V + /illV[5^5^r|+iV), (4.19) 

for = \ (eV'^e±^^"C"^)- 

There are other operators formally at this order which do not contribute to the polariz- 
abilities at the order to which we work. There are kinetic operators and higher dimensional 
couplings of the baryons to the axial current whose coefficients are exactly fixed by the repa- 
rameterisation invariance of the baryon four-momentum |146|I144] . These operators give the Z 
dependent pieces of the Compton amplitudes in Eq. ()4.2() . There are also additional operators 
with unconstrained coefficients such as (B A ■ AB^ that contribute to the Compton amplitude 
at higher order. In two-flavor there are two such operators, and in the SU(4|2) case there 
are ten jl44j . 

The leading operators which contribute to the electromagnetic polarizabilities at tree level 
occur at 0{Q'^) and are given by the general form, 

X 

(where the F^^ are spin structures) while the leading tree-level contributions to the spin polar- 
izabilities occur at 0{Q^). The complete set of such operators in the case of two- flavour j(PT is 
given in Ref. j211j . Again there are significantly more such operators in Q^PT and PQxPT. We 
do not explicitly show these operators, as they do not contribute at the order we are working 
and will not modify volume dependence until 0{Q^). 

4.5 Nucleon Compton Scattering in Finite Volume 

Using the Lagrangian of the preceding section, we can calculate the amplitudes defined in 
Eq. 1)4. 1() for Compton scattering from a nucleon (extensions to full octet and decuplet of 
baryons are straight-forward although the convergence of HBxPT with three-flavours is not 
clear). We work with a power counting such that 

g~e~M~^~^ (4.20) 
A^ A^ 

(it is also convenient to count A/A^ as the same as Q as it is numerically similar at the masses 
relevant for current lattice calculations).^ Below, we will also restrict ourselves to the low 

^Loop and pole |107| contributions with 44-plet intermediate states must be included since A is a small-scale. 
Any A dependent terms analytic in arising from the loop diagrams, and additional operators proportional 
to powers of can be resummed into the appropriate LECs of A independent operators (the LECs then 

depend on A) |212l I144| . Keeping these contributions explicit is redundant as A can not be varied in a 
controlled manner. 
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Figure 4.2: Anomalous contributions to the polarizabilities. The crossed circle corresponds 
to the insertion of an operator from the Wess-Zumino-Witten Lagrangian, Eq. (|4.25|) . and the 
crossed meson line corresponds to a hairpin interaction |213j . 




Figure 4.3: Born diagrams involving internal 44-plet states that give contribution to the 
polarizabilities. 



frequency limit uj <C m.j^ in order to extract the polarizabilities from the Compton scattering 
amplitudes defined in Eqs. (|4.H) and (|4.2|) . For larger energies, the concept of polarizabilities 
breaks down and the target essentially becomes a dispersive medium. Working to order 
in the chiral expansion, Compton scattering requires the calculation of the diagrams shown in 
Figs. 14.21 14.3l and l4.4l (and a corresponding set involving internal 44-plet baryons). By definition, 
tree level contributions from nucleon pole diagrams do not contribute to the polarizabilities; 
their contribution to the amplitudes are given explicitly in Eq. (|4.2|) . For each polarisability 
X = a, (3, 7i-4, it is convenient to separate the different contributions as 

X = X™^^y + X^ + Xi°°P , (4.21) 

corresponding to the contributions from Figs 14.21 14.31 and 14. 4| respectively. We discuss these 
contributions in the following subsections. At order Q^, all contributions are expressible in 
terms of a small set of LECs that contribute in many other processes and are thus reasonably 
well determined (at least in the xPT case). The total 0{Q^) loop contribution is finite, but 
loop-contributions at higher orders are divergent; as discussed in the preceding section, the 
counter-terms specific to Compton scattering that absorb these divergences and the associated 
scale dependence enter at 0{Q'^) for the electric and magnetic polarizabilities and 0{Q^) for 
the spin polarizabilities. 
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(a) 




(c) 



(d) 



(e) 




(f) 






(h) (i) 

Figure 4.4: Diagrams contributing to nucleon polarizabilities at order Q^. The solid and dashed 
Unes denote nucleons, and pions respectively. Additionally, our results include a similar set of 
diagrams in which the internal 70-plet propagator is replaced by a 44-plet resonance. 



4.5.1 Volume independent contributions to polarizabilities 



Anomalous contribution to jN ^yN: vr^ 
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The anomalous decay of flavour neutral mesons to two photons [2141 1215j has important con- 
sequences in Compton scattering in non- forward directions. These contributions arise from the 
meson pole diagram shown in Fig. 14.21 Anomalous decays are well understood in xPT, entering 
through the Wess-Zumino-Witten (WZW) Lagrangian [2U81 l2Ut?] . However, these effects have 
not been investigated in the quenched and partially quenched theories and some interesting 
subtleties arise. 

In SU(2) xPT, the one-pion, two-photon piece of the WZW Lagrangian, is given by 



where 



3e^ 



167r2/ 



vr 



tr 



Q 



nvpa TP p 
c -L /iiy-t per. 



V2/ 



(4.22) 



(4.23) 



This Lagrangian is completely determined as its coefficient can be fixed by directly matching 
to the perturbative QCD calculation of the relevant triangle diagram (the one loop calculation 
is exact |21Hj . in accordance with Witten's geometric quantization condition |2()9j ). At higher 
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orders, additional anomalous operators appear |217j but they do not contribute to Compton 
scattering until 0{Q^). 

It is well known that quenched and partially-quenched chiral perturbation theories generally 
have more complicated operator structure than in the case of QCD {e.g., one can not use Cayley- 
Hamilton identities jl32j ). Thus, in order to generalize Eq. ()4.22|) to the partially-quenched 
cases, we might imagine the extended vr'^ — > 77 Lagrangian to be of the form 



ci str [$Q2j _^ ^2 str [$Q] str [Q] 



C3 str [^>] str [Qf + C4 str [^>] str [Q^] 1 , (4.24) 



(in the quenched case only the first operator is non- vanishing, but a similar discussion applies). 
With the condition that in the QCD limit where the sea-quark and ghost-quark masses and 
charges are set equal to those of the valence-quarks, matrix elements of Eq. (|4.24|) reproduce the 
matrix elements of Eq. 1)4.22(1 . As discussed in Section f4. 41 the singlet field, $0 acquires a large 
mass from the strong U(1)a anomaly j218j and can be integrated out of the partially-quenched 
theory; consequently, the operators proportional to 03 and 04 can be ignored. Additionally, 
from the multiple super-trace structure, one can deduce that the operators 



str str [Q] , str [$] str [Qf , str [$] str [Q^] , 



have at least two closed loops at the quark level. Following the arguments in Refs. [2141 12151 
l21Hj . one can show that these operators do not correspond to anomalous quark level processes. 
Moreover, the leading dependence of the underlying quark-level diagrams is proportional to the 
quark mass, and thus the coefficient of these operators must scale as, 02,3,4 ~ "T-g/^x" Although 
these operators contribute to r]a — > 77, they are not anomalous, and only contribute at higher 
orders in the chiral expansion. We can thus conclude that the only operator in the Lagrangian, 
Eq. (|4.24() . which contributes to the anomalous decay of the neutral mesons at leading order 
is str [^>Q^]. The coefficient is easily determined by matching to either perturbative partially- 
quenched QCD or to the x^T expression in the QCD limit.'' The appropriate Lagrangian is 
therefore 

<t = -T|4^tr[$Q2]e'^^^-F,.F,.. (4.25) 



From the above Lagrangian, it is apparent that all of the flavour diagonal fields in Eq. ()1.59() . 
have anomalous couplings to two photons. Calculating the diagrams in Fig 14.21 leads to the 
following anomalous contribution to Compton scattering on the proton^ in partially-quenched 



^We can thereby bypass the need to extend Witten's global quantization condition |2U9| to non-compact 
graded Lie groups. 

*The anomalous contribution to neutron- Compton scattering is given by Eq. 14.261 with the interchange of 
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T-iPQ, anomaly 



-I e 



k"" kf^r-S 



29A 



2* 



+ 51 



+ 



q] - if 



mt 



~ 2 I 2 z 

+ ~ 



dd 



It 



(r^ — m\) \ {r 



+ 



rat 



(r2 



m 



dd' 



(4.26) 



In the above expression, r = q' — q, is the momentum transfer to the nucleon and Afj = mfi — injj 
is a measure of the isospin breaking in the sea sector. In the sea isospin hmit (mi — > mj), the 
double pole structure of the amplitude vanishes, and in the QCD limit all dependence on gi 
vanishes. 

Expanding Eq. 1)4. 26(1 in frequency and comparing with Eq. (|4.1[) leads to the following 
anomalous contributions to the polarizabilities: 



^anomaly 


= 0, 


(4.27) 


^anomaly 


= 0, 


(4.28) 


anomaly 


>Jt; >-^anom 


(4.29) 


^anomaly 


= 0, 


(4.30) 


anomaly 


"-^anom 


(4.31) 


^anomaly 


■Jc; "-Tanom 


(4.32) 



where the coefficients, Ganom; are given in Table ITT] fat the end of this chapter) for the different 
theories under consideration. These contributions vanish in the iso-scalar combination of proton 
and neutron targets in the QCD limit. 



A resonance contributions 



The contributions to the amplitude from the Born-terms involving the 44-plet resonance (which 
contains the A-isobar), Fig. 14.31 are identical in xPT, PQ^PT and QxPT as they are purely 
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valence quark processes. They are given by 



a 



7f 
1^ 



it 
it 



0. 



367r(2Miv)2A ' 



2 e^iqu-qdf 



727r(2MAr)2A2 



0, 

2 e^{qu-qdf 



727r(2M7v)2A2 ' 

where is the magnetic dipole transition couphng of Eq. (|4.16|1 . 



(4.33) 
(4.34) 
(4.35) 
(4.36) 
(4.37) 
(4.38) 



4-. 5. 2 Infinite volume 

The loop contributions to infinite volume chiral expansion of the polarizabilities in \PT are 
well known at order [m [TOni HH CIl HSl HTZl and at 0{Q^) 1121 1^ HH 
I225( I226| I115| I118j (at this order, the A-resonances have not been included as dynamical de- 
grees of freedom, restricting the range of applicability to ^ A). Since the photon only 
couples to charged mesons, the results in the quenched and partially-quenched theories are 
similar to those in xPT. In particular, no quenched or partially-quenched sicknesses (double 
pole contributions from neutral meson propagators) enter expression for the loop diagrams. 
In general, the quenched power counting presents differences for electromagnetic observables 
|139[ 1205} I227[ l2T?1] . however no new contributions appear at the order we work. 

Using the eff'ective couplings Gb-, G'^, Gt and G'rp given in Table we find that the loop 
contributions to the polarizabilities are 



^loop 
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Here we have used dimensional regularization, however the results are finite and hence indepen- 
dent of the regulator without the addition of counterterms. These loop contributions vanish at 
zero photon frequency, preserving the Thomson limit. They are identical for both proton and 
neutron targets, though isospin breaking effects from loops enter at 0{Q^) in the expansion. In 
the xPT case, these results reproduce those of Refs. |113( 1114] . 

4-5.3 Finite volume 

In momentum space, the finite volume of a lattice simulation restricts the available momentum 
modes and consequently the results differ from their infinite volume values. These long-distance 
effects can be accounted for in the low-energy effective theory. Here we shall consider a hyper- 
cubic box of dimensions x T with T ^ L. Imposing periodic boundary conditions on 
mesonic fields leads to quantized momenta k = (/cq, k), k = = ^(ji, ^'2, J3) with G Z, but 
ko treated as continuous. On such a finite volume, spatial momentum integrals are replaced by 
sums over the available momentum modes. This leads to modifications of the infinite volume 
results presented in the previous section; the various functions arising from loop integrals are 
replaced by their finite volume (FV) counterparts. In a system where 777,^^ ^ 1, the power 
counting of the infinite volume low-energy effective theory remains valid and finite volume 
effects are predominantly from Goldstone mesons propagating to large distances where they 
are sensitive to boundary conditions and can even "wrap around the world". Since the lowest 
momentum mode of the Goldstone propagator is ~ exp(— 777^L) in position space, finite volume 
effects will behave as a polynomial in 1/L times this exponential if no multi-particle thresholds 
are reached in the loop (as is the case in these calculations provided the photon energy is small 
enough, to < tr^t). 

Repeating the calculation of the loop diagrams using finite volume sums rather than integrals 
leads to the following expressions for the loop contributions to the polarizabilities: 
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(4.54) 



(4.55) 
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(4.57) 
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(4.59) 



and the finite volume sums Ifs{m), . . . , C/^{m) are defined in Appendix IB. 21 These expressions 
reduce to the results of Eqs. ()4.39I) - H4.44p above in the infinite volume limit. 

To illustrate these effects, Figs. 14.51 and 14.61 show the volume dependence of the various 
polarizabilities in the proton and the neutron, respectively. Here we have specialized to QCD, 
setting qu = 2/3, qa = -1/3, gA = 1.25, Iqna] = 1-5, ^lT = 10.9, / = 0.132 GeV, Mn = 
0.938 GeV and A = 0.3 GeV.^ In each plot we show results for the ratio 

X{L)-X{oo) 



AX{L) 



X{oo) 



(4.60) 



^The value of fir is chosen to correspond to that found in analysis of Ref. |114j {^t = 2V2bi of that reference). 
In principle this LEG can be determined from an analysis of lattice polarizabilities or A'^-A transition matrix 



elements. 
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Figure 4.5: Volume dependence of the proton polarizabilities. Here we show the ratio of the 
difference of the finite and infinite volume results to the infinite volume results for three values 
of the pion mass using the parameters described in the text. The curves terminate at rriT^ L = 3. 



for the six polarizabilities at three different pion masses, = 0.25, 0.35, 0.50 GeV. The overall 
magnitude of these shifts varies considerably; generally volume effects are at the level of 5-10% 
for = 0.25 GeV and smaller for larger masses. Larger effects are seen in a number of the spin 
polarizabilities but there are as yet no lattice calculations of these quantities. The magnetic 
polarisability has a particularly small volume dependence which can be understood from the 
large decuplet resonance contribution that is independent of the volume. 
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Figure 4.6: Volume dependence of the neutron polarizabilities. The various curves are as in 
Fig. 1131 



The above expressions also allow us to calculate the finite volume effects in the quenched 
data on the various polarizabilities calculated in Refs. |19H 1192] . The quenched expressions 
involve a number of undetermined LECs (quenched qa-, gi, On A and /^t are unrelated to their 
PQxPT/xPT values), so we can only estimate the volume effects. To do so, we choose Qu = 2/3, 

= -1/3, gA = 1.25, gi = 1, \gNA\ = 1-5, fiT = 5.85, / = 0.132 GeV, Mn = 0.938 GeV 
and A = 0.3 GeV and take the pion masses corresponding to the lightest used in these lattice 
calculations, ttItt ~ 0.5 GeV (we ignore the issue of the convergence of xPT at such masses). The 
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Figure 4.7: Volume dependence of the proton polarizabilities in quenched QCD at the hghtest 
quark mass used in the lattice calculations of Refs. jl911 11912] . The central curves and shaded 
region correspond to the parameters quoted in the text. 



results for the volume dependence of the various polarizabilities of the proton and neutron are 
shown in Figs. E77I and ITHl In each plot, the shaded region corresponds to reasonable variation 
of the unknown couplings, \ < qa < 1.75, —\<g\< 1, 0.8 < I^ataI < 2 and 2.8 < \[it\ < 8.5. 
From these figures, we see that the calculations on a (2.4 fm)^ lattice with rn-jr =0.5 GeV may 
differ from their infinite volume values by 5-10% in the case of the electric polarisability and a 
few percent for the magnetic and spin polarizabilities. 
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Figure 4.8: As in Fig. 14. 7( but for the neutron. 



An interesting effect that arises at finite volume is that the Thomson-limit and other Born 
terms in the frequency expansion of the scattering amplitude (terms in Eq. ()4.2() that are not 
polarizabilities) receive finite volume contributions from the loop diagrams in Fig. 1131 that 
vanish exponentially as the volume is increased. As an example, the amplitude for Thomson- 
limit (zero frequency) scattering on the neutron (which is identically zero at infinite volume) is 
shown in Fig 14.91 This result is somewhat counter-intuitive, but arises from the effects of the 
periodic boundary conditions on the long range charge distribution of the hadron. It does not 
imply the non-conservation of charge. One way to understand this feature is that in a finite 
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Figure 4.9: Volume dependence of the Thomson limit of photon neutron scattering. 



box with periodic boundary conditions, only quantized momentum are allowed, including the 
momentum of the photon field. The Thomson limit is defined as the limit of a zero-frequency 
photon scattering off a charged particle, and thus can only resolve the total charge of an object. 
However, to take the limit of zero-frequency, one must first take the box size to infinity, such that 
one can smoothly vary the photon energy to zero. But at infinite volume, as seen in Fig. 14.91 
the deviation of the zero-frequency limit scattering from the naieve infinite volume Thomson 
limit goes to zero. Thus, this is another example where the order of limits is important, and in 
this case one must take the infinite volume limit before the zero-frequency limit. 

The results presented here all assume that the higher order terms in the Q expansion provide 
small contributions to the volume dependence of the polarizabilities. This may or may not be the 
case as diagrams that are formally of higher order in the infinite volume power-counting can 
have volume effects that are enhanced over those at lower infinite volume order (see Ref. j2U7j 
for a detailed discussion) . Such issues may be particularly relevant for the polarizabilities where 
the convergence of the chiral expansion is tenuous. In this regard, studying the FV behaviour 
of the lattice results may in fact be a useful diagnostic tool with which to determine if or why 
the convergence is poor. 
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4.6 Conclusion 

We have investigated Compton scattering from spin-half targets from the point of view of lattice 
QCD. We first discussed how external field methods can be used to probe all six polarizabilities 
of real Compton scattering for both charged and uncharged targets. Such calculations will tell 
us a lot about the low energy QCD structure of hadrons and will be of much use in phenomeno- 
logical studies requiring the full set of polarizabilities as only certain linear combinations are 
available from current experiments jl88l I189j . The techniques discussed here also allow us to 
extract other electric properties of charged particles using external fields including the electric 
dipole moment of the proton and the quadrupole moment of the deuteron. 

Our second major focus was on the effects of the finite volume used in lattice calculations 
on the polarizabilities. Since polarizabilities are infrared-sensitive observables (they scale as 
inverse powers of the pion mass near the chiral limit), the are expected to have strong volume 
dependence. This is indeed borne out in the explicit calculations presented here. In QCD, 
we generically find that the polarizabilities experience volume shifts of 5-10% from the infinite 
volume values for lattice volumes ~(2.4 fm)^ and pions of mass 0.25 GeV. The electric and first 
spin polarizabilities are particularly sensitive. In the case of quenched QCD (relevant to the 
only existing lattice data), we find significant effects even at pion masses ~ 0.5 GeV. Future 
lattice studies of the polarizabilities should take these effects into account in order to present 
physically relevant results. 

As extensions of this work, one can also consider the generalized polarizabilities, higher-order 
polarizabilities and parity violating polarizabilities (see Ref. |228j ) all of which can be extracted 
from appropriate lattice calculations similar to those detailed in Section f4. 31 Such information 
would lead to a further-improved understanding of the low-energy structure of the hadrons and 
prove very useful in directing the next generation of precision Compton scattering experiments. 
The lattice provides a novel opportunity to study the neutron polarizabilities directly instead 
of from nuclear targets and extending the lattice methods of Section [4.31 to the deuteron (along 
similar lines to those discussed in Ref. |229j ) will also prove useful for comparison to experiment. 



Table 4.1: Effective couplings for the various contributions to the polarizabilities. 
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Chapter 5 

TWO-HADRON INTERACTIONS: 
TTTT SCATTERING IN A FINITE VOLUME 

The work in this chapter is based upon Refs. |23Uj 

5.1 Introduction 

Except in the case of infinitely heavy baryons, where an adiabatic potential can be defined, 
the interaction between two hadrons is studied with lattice QCD by numerically calculating 
energy levels of the system in a finite box. This is because in the infinite volume limit and 
away from kinematical thresholds, the two-hadron Euclidean correlator gives no information 
about the Minkowski space amplitude [10] • The alternative is to consider the system in a 
universe with finite spatial extent (a finite box), as is the case with numerical calculations 
anyway. The energy levels of a system composed of two hadrons are not simply the sum of 
the energies carried by each hadron, but there is an additional (usually small) shift that arises 
due to the interaction between them. The smaller the box, the larger the shift in energy 
levels. This volume dependence is inversely proportional to the volume and furthermore, there 
is a relation between the energy level shifts and the scattering phase shifts [371 1391 138j . This 
relation, valid for energies below the first inelastic threshold is a consequence of unitarity and 
is thus model independent.^ In addition to this power law shift in the energy levels, there are 
exponentially suppressed corrections which are not model independent and are the analogue 
of the exponentially suppressed corrections to the mass, decay constants, etc., in the single- 
hadron sector |23H I232| I233| I234| I235j . These exponential volume effects arise because the 
off-shell propagation of intermediate states is altered by the presence of the finite box, which 
allows them, for instance, to "wrap around" the lattice. As such, these effects are dominated 
by the lightest particle, the pion in QCD, and are proportional to e"™""^ with m^r the pion 
mass and L the linear dimension of the box. For simulations done with small enough quark 
masses such that the pions are within the chiral regime, these soft pion effects can be computed 
using the chiral perturbation theory (xPT) |701 171j . The vrvr scattering phase shifts have been 
computed using lattice QCD following the universal finite volume method mentioned above |236[ 
IMl USa HSni 12101 Hm HH HH HH mSl mi mi • As the chlral Umlt is approached jUl and 
more precise calculations appear, these exponentially suppressed corrections will need to be 
understood. 

Our goal in this paper is first to show the modification of the universal scattering formula for 
a hadron-hadron system in a box due to the exponentially suppressed finite volume corrections. 
Second, we compute the dominant exponential volume dependence explicitly for the two-pion 



^By model independent relation we mean a relation valid whether one is considering QCD or some other 
theory, as long as this theory obeys unitarity, locality, etc. 
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(d) (e) 

Figure 5.1: The one-loop diagrams which contributing to the vrvr scattering amplitude. Only 
the s-channel diagram, (a), contributes to the power-law volume dependence. Diagrams (b) and 
(c) are the t-, and n-channel diagrams, respectively, while diagram (e) represents wavefunction 
renormalization. All these diagrams contribute to the exponential volume dependence. 



system in / = 2 channel near threshold by use of the leading loop-order two-flavor xPT. 
5.2 Finite Volume vrvr Scattering 

5.2.1 Power law and exponential volume dependence 

As discussed above there are two types of volume dependence of the energy levels of two hadrons 
in a box: power law (proportional to ^/L^) and exponential (proportional to e~™""^). The 
first is exploited by the finite volume method to extract information about scattering parame- 
ters |371l39ll^ . The second, usually numerically smaller, appears as a correction to the relation 
between energy levels in a box and scattering parameters. In order to compute the exponentially 
suppressed terms we need to separate them from the larger power law contribution. 

Figure 15.11 shows all vrvr scattering diagrams which contribute at one- loop order. As we 
will discuss in more detail in the next section, the power law corrections arise only from s- 
channel diagrams as shown in Fig. 15.11 (a). where the intermediate particles can be on-shell, and 
thus propagate far and "feel" the finiteness of the box. In all other diagrams the intermediate 
particles are very off-shell, cannot propagate farther than a distance of order l/rriT^ and therefore 
have only small, exponentially suppressed sensitivity to the size of the box. 

Let us now discuss the general form of finite volume corrections. Consider first the pion 
propagator at finite volume which is a function of the spatial momentum k = 27rn/L and the 
energy E. It will have poles for values of E corresponding to the values of the energy of a pion in 
the box. In particular, for A; = the pole will be at mTj{L) (the "finite volume mass"), differing 
from the (infinite volume) mass of the pion m^r by an exponentially small quantity proportional 
to (m2/(47r/)2)e-™-^/L^/^ f^l^SSlllSSl- The extra suppression factor rn^ / [An f-j^)'^ is due 
to the fact that only loop diagrams contribute to the finite volume corrections. 
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The volume corrections for systems with more than one hadron are more subtle. The reason 
is that there are two kinds of volume corrections to the energy levels: a power law one described 
by the Liischer formula and the exponentially suppressed ones. To understand how to separate 
them let us first look at the infinite- volume, S-wave scattering amplitude, T{s), with on-shell 
external pions. It is given at one loop by^: 

r(.) ^ T(')(^s) + 4^{s)+T^l{s) + ^T^){s), (5.1) 

and unitarity allows us to re-sum this amplitude as a bubble series, 

- r(o)W -ri>) -</(.)■ 

where T^") is n-th loop contribution, the s-channel contribution is separated into its real part, 
rj^(s), and imaginary part, Tjy'(s), and all other contributions at one- loop including t- and 

M-channels are denoted by T^]j{s). The imaginary part, T^^j{s), comes from picking in the 
loop integration, the particle poles in both pion propagators, such that both the loop pions are 
on-shell. The loop integral is then proportional to the phase space volume and is given by 

^i' <•> = ^^V^. (5.3) 

The fact that the imaginary part is determined by the tree level amplitude is a consequence of 
the optical theorem. 

It is useful to define the X-matrix |247j . which at one-loop is given by: 



K{s) ^ T('\s) + Ti'^{s) + T^iis) ^ ^'^Z^'^^' . (5.4) 

Since the scattering amplitude can be written in terms of the phase shift, S{s), as 

^ 32vrVi 1 ^ 32vr^ e^^W - 1 

- 4ml cot 6{s) - i - 4ml ' 

the relation between the K-matiix and the phase shift is then given by 



1 1 Is 4mJ.^^^^^^^_ ^^g^ 



K{s) 327r V s 

Now we look at the finite volume amplitude T{s)^. It is computed in the same way as the 
infinite volume amplitude, except the loop integrals are substituted by sums over the momenta 

■^We are considering the s-wave projected amplitude and disregarding the mixing with higher partial waves 
induced by the breaking of rotational symmetry. 

''By finite volume scattering amplitude we mean the amputated four-point correlator since, of course, there 
is no scattering at finite volume. 
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allowed in a finite box. The important point to keep in mind is that sums where the summand is 
regular are, at large enough L, well approximated by the analogous integral, up to exponentially 
small terms. If the summand, however, contains a singularity, power law dependence on the 
volume arises. As mentioned before, only the kinematics of the s-channel diagram allows for 
both of the intermediate pions to be on-shell simultaneously. This implies that the summand 
in the sum over the loop momentum contains a singularity and leads to power law volume 
corrections. For the remaining diagrams no singularities are present and only exponentially 
suppressed corrections can arise. As it will be shown explicitly below, the finite volume ampli- 
tude then has the following form: the tree term remains the same, the t-, u- and the real part 
of the s-channels pick only exponential corrections but the imaginary part turns into the term 
with power law L-dependence (and is real at finite L): 



T[s) ^ T(0)(.) + 7;(i)(.)+Ta(.) + A7ii)(.) + g{^ 



{T^){s)f 



(5.7) 



where s = 4{k'^ + m^), ATexp(s) is the finite volume correction to T^]}{s) + T^^^{s), and 5 is a 
universal (independent of the interaction) function of s [421 143j : 



r3 



L3 ^^^^ J (2^s 



= lim ^^7:5- — 47r A„ . (5.8) 



|n|<An 1^ 



Combining Eqs. (|5?2|) . (fOI) and (|5?T|) . we have 
T{s) ~ ^ 



_l AtW(^) _ 1 ^ ( {s-iml)L^ 

K{s) (T(0)(s))2 Wtt^L^"^ \ 

16'^■^/s 

k cot Sis) - l6.V-s0f^ - izS (^^^" 



(5.9) 



The energy of the states in the box are determined by the location of the poles of the finite 
volume amplitude determined by the solution of 

fccot 5{s) - 167r^— ^j^^— - = —S —2 . (5.10) 



(r(o)(s))2 vrL V 16^^ 
We recognize in Eq. ()5.10|) the familiar form of the Liischer relation modified by the finite volume 

at'"*"' is') 

correction: the quantity — IGvr (-r(oT(s))2 sought-after (exponentially small) correction 

to k cot 5{s). In this work, we will focus on the 2-pion correlator near threshold in the center 
of mass frame, for which the infinite volume energy is given by ^/s = 2m^. The solution of 
Eq. (|5.10j) . s* , will be away from threshold by an amount given by — 2m^ f« y/s* — 2mT^{L) ~ 
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l/f^L^. Therefore, for s ^ s* the correction term, A{k cot S{s)) = —16'^■^/sJ;^^^^^, can be 

approximated by —32Trm.^-^^j^^^^, the difference being suppressed by 

It is customary to expand Eq. ()5.10() in powers of k"^ ~ Up to the first three orders of 

this expansion (near threshold), kcot6 can be approximated by the inverse scattering length, 
1/a, resulting in 

^/;^-2m. = ^(^l + c,- + c,(-) (5.11) 

where ci_2 are known numerical factors. A generahzation of this formula including the expo- 
nentially suppressed corrections, however, is not useful. The error in using Eq. ()5.11|) instead 
of Eq. (|5.10|) . that is, the error in the extrapolation from s* to 4m^, is of order 1/L^, which 
is parametrically larger than the exponential corrections we are interested in. Numerically, it 
may be the case that, for a set of simulation parameters, the exponential term is larger than the 
terms. But the analogue of Eq. (|5.11|) one would obtain by formally counting, for instance 
^-m^L ^ 1/L^, would involve the effective range in addition to the scattering length. In any 
case, it is unclear a priori that a set of simulation parameters exist where this kind of expansion 
is useful and we will not pursue this line of thought in this paper. We now will compute the 
exponential corrections to k cot 6{s). 



5.2.2 The vrvr scattering amplitude 



The vrvr finite volume correlator in the 1 = 2 isospin channel, for arbitrary momentum (in the 
chiral regime) is given by 



2 ({3s-2m^-Et=iPl 
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32 



(5.12) 



In the above expression, m is the (volume independent) tree level pion mass, / « 132 MeV is the 
(volume independent) tree level decay constant and the £j's are the Gasser-Leutwyler coefficients 
of counter terms appearing in the chiral Lagrangian at next-to-leading order (NLO) j7Uj . The 
loop integrals/sums are given by T, J^{P), J'^{P), and J^u[P) and will be defined below shortly. 
The external momenta pi, i = 1 • • • 4 are described in Fig. 15.11 and the Mandelstam variables 
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are employed: s = with ps = pi +P2, t = Pf with pt = pi —ps, and u = p'^ with pu = pi —pA- 
In the above equation, the first term is the leading-order (LO) tree level contribution, and the 
remaining terms come from the one-loop diagrams shown in Fig. 15.11 and tree level diagrams 
with vertices of the 0{p'^) Lagrangian j7()|l71j. In the NLO contributions, we have approximated 
pf = m^, as the corrections to this at finite volume are beyond the order we are working. In 
the one- loop terms, we have expressed all contributions in terms of the bare pion mass as the 
difference is higher order in the chiral expansion. One can, of course, choose to express the 
scattering amplitude in terms of the "lattice quantities" such as m.,^{L) and fniL), which are 
measured directly from Euclidean correlation functions by lattice simulations. Converting the 
bare quantities into lattice quantities involves additional tadpole loops, which will affect the 
form of AT^lxp , the finite volume corrections to the two-pion amplitude. The exponential volume 
dependence of course doesn't depend upon whether one expresses the amplitude in terms of 
either the bare or physical parameters, and so it is useful to use the form which is simplest, 
that in terms of bare parameters. In what follows, we will be interested in the 2-pion correlator 
near threshold, for which the external pion momentum are given by pi ~ {^^/s,0). 

5.2.3 Loop integrals / sums at one-loop 

The loop integrals/sums appearing in Eq. 1)5. 12() are defined by 

(5.13) 



L 
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1 


1 


1 


q2 _ j^2 
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q2 _ j^2 


(P -t- qY — rn? 




1 



JAP) - I ll ^3 Y^^^ qTT^2 (p + qY _ ^2 ■ (5-16) 



Note that an integral is taken along the 0**^ component whereas sums over discrete momenta 
are taken with cubic symmetry. Finite volume effects in the loop integrals/sums in Eqs. (|5.13j) - 
1)5. 16() can be computed by first evaluating the go contour integral and then using the Poisson 
resummation formula. 



The difference between the finite volume and infinite volume loop integrals/sums is given by 
the second term in the right-hand side of Eq. (|5.17|) . and is always ultraviolet finite. If the 
function f{q) is regular, this difference is exponentially suppressed in the large L limit. Power 
law dependence on L can however appear if f{q) has a singularity, i.e., the case when P = pg. 
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Let us now evaluate the difference between the finite and infinite volume integrals/sums 
given in Eqs. (|5.1I-{|) - (|5.16)1 . We shall define the difference between the finite and infinite volume 
integrals/sums as 



1=^ 



_fq_ 

(2^)3 



(5.18) 



where it is implicit that we regulate the sum and integral in the same manner such that the 
UV divergences cancel. The tadpole integral in Eq. (|5.13|) . which contributes to and /,r and 
the loop diagrams of Fig. 15.11 has the following volume correction; 
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(5.19) 



where ujq = \/ q'^ + rn^- The mass and decay constant measured in lattice simulations are thus 
given to NLO by dSDHS]^ 
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(5.20) 
(5.21) 



Using the asymptotic form of the Bessel function, one can see that for large mL, the volume 
shift of the pion mass is exponential 



Ami _ ^^2: 



1 



772 2^ 



ml /I 25/2vr3/2L/2 ^ n^l^ 

n=\n\^0 



\ 3 1 

8 nm-j^L 



(5.22) 



where the ellipses denote more terms in the asymptotic expansion of the Bessel function and c{n) 
is the multiplicity factor counting the number of times n = |n| appears in the 3-dimensional 
sum. Note, this sum is not over integers, but rather over the square-roots of integers. In 
Table 123 we list the first few values of the multiplicity factors. 



Power law L-dependence can only occur through the integrals/sums in Eqs. (|5.14|) - (|5.16j) 
when p2 > 0. For the center-of-mass scattering kinematics we are considering here this can 



*This relation is known up to two loops |235l 12481 - 
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Table 5.1: Here we list the first few multiplicity factors which arise when converting the 
three-dimensional sum to a scalar sum. 



n 


1 ^2 


VS V4 


V5 


V6 


V7 




V9 


Vw 


c(n) 


6 12 


8 6 


24 


24 





12 


30 


24 



only occur for P = Ps, since = s > 0. As argued above we will only need the amplitude 
at threshold, i.e., Ps = (2/71-7^,0) and pt = Pu = 0, except for the terms with power law L- 
dependence. Consequently, we will need only the values of AJ'{P = 0), AJq{P = 0), and 
^Joo{P = 0) for t- and ■u-channels as well as AJ'{P = Ps), AJq{P = Ps), and Aj7oo(-P = Ps) 
for s-channels. The JT" integrals/sums at P = can be shown to be related to 1, giving the 
volume difference: 
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(5.24) 
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^AJoo(O) = z / ^ 
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(5.25) 



The power law volume dependence appears in the remaining integrals/sums. In those we 
keep s away from the threshold value and take ^/s = 2\/A?^~+~m^. After performing the qo 
integral, we separate the singular piece of the summand from the rest as 



iJ{Ps) = -773 V — ^2~~ 
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1 1 J_ V 

„ r,3 ^2 _ 1,2 + AT?, ■ 



1 UJa 
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(5.26) 



The first term contains a singularity when the internal momentum coincides with the external 
momentum, while the second term is regular. The difference AJ'{ps) is then 



iAJ{p,) = - 



g 2 - A;2 4a;fca;g 



(5.27) 



The first piece above is the promised universal function containing the power law volume depen- 
dence. The summand in the second term contains only exponential finite volume corrections. 
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This term, contributing to AT^xp, can be computed at the s = Am?, k = threshold point, 
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, (5.28) 



where Ly is the Struve function. To get the second line of Ea. H5.28() . we have neglected terms 
which are suppressed by 0{k'^ /m?) relative to the first. For the two-pion system, this is ap- 
proximately given by ^ ~ trfiiji ^- '^^^ asymptotic expansion of iA^exp(Ps) is given by 



iAJ'exp(Ps 
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(5.29) 



Again, one can see that these volume corrections to the integral are exponentially suppressed. 

The finite volume dependence of the other s-channel loop integral functions, iAj7o(ps) and 
iAj7oo(Ps) become simpler to evaluate by first observing that the summands can be separated 
into the following pieces: 
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One then obtains, 
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(5.30) 



lAJooips) = iAI 
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(5.31) 



Having these tensor integrals/sums written in terms of the scalar integrals/sums iAJ'{ps) 
and iAI and using the scattering amplitude in Eqs. 1)5.2(1 and l|5.12() . one can now verify that 
the coefficient of 5(^^) in the amplitude is what was promised in Eq. ()5.2|) . 



5.2.4 Exponential volume correction to the / = 2 vrvr correlator 



Collecting the results for the sums/integrals in Eqs. 



and using the amplitude in 



Eq. (j5.12j) we can now compute the correction term in Eq. (j5.1()|) for 1 = 2 two-pion system 
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Figure 5.2: Ratio of the magnitude of the exponential correction term A(/ccot(5) 
— 327rm^ (r"(oT(4m?))2 to A; cot 5 for different values of the pion mass. 



near threshold. We find 
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Equation 1)5.32(1 is our main result (the first line being the exact one-loop answer and the 
second line the asymptotic expansion in mj^L). In this expression, one can use either the bare 
parameters, the physical parameters or the finite volume parameters as the difference is higher 
order than we work in either the chiral expansion, or in the exponential dependence. It is most 
convenient to use the values oimn{L) and fniL) directly measured in a given lattice simulation. 
In Fig. I5.2l we plot the ratio of A(k cot S{s)) to the one-loop value of k cot 6{s) using Eq. (|5.32() 
as a function of L for some reasonable values of mT^. We find the finite volume corrections to be 
relatively small, a few times smaller than the statistical and systematic errors quoted in recent 
simulations. An error of about 10% was quoted in reference 0ZI for the determination of the 
scattering length for a pion mass of rriT^ ~ 290 MeV and a box size of L ~ 2.5 fm.^ The finite 



^In Ref. \47l, Beane et. al. determined the 1 — 2 tttv scattering length for various pion masses using a mixed 
action simulation with Domain- Wall valence quarks and staggered sea quarks |249| . Because of the mixed- 
action, the mesons composed of sea quarks and the mesons composed of valence quarks receive different mass 
shifts from the finite lattice spacing. This means that even when the sea and valence quark masses are tuned 
equal, there are still partial quenching effects in the simulation. In Ref. |137| . the partial quenching and lattice 
spacing corrections to the I = 2 nn scattering length were worked out for this mixed action theory. It was 



Ill 



volume correction from Eq. (|5.32|) for these parameters is approximately 1%. These corrections 
however grow fast with the approach to the chiral limit, and they become non-negligible as 
smaller pion masses are used and statistical errors are reduced in simulations. 

5.3 Discussion 

We have described the leading exponential volume dependence explicitly for the scattering 
parameter of a / = 2 two-pion system near threshold in a box, by extending the one-loop 
xPT calculation of pion scattering [JHI to include the volume dependence. The exponentially 
suppressed volume corrections can distort the universal relation between the infinite volume 
scattering parameters and the power-law volume dependence of the two-particle system, espe- 
cially as the chiral limit is approached. An important point we want to stress is that the useful 
way to add the exponential volume dependence to the relation between infinite volume scat- 
tering parameters and the energy of the two-particle system in a box, is via Eq. 1)5. 1U() . which 
allows an understanding of the leading exponential volume dependence to k cot S{s). This is 
contrast to the notion of studying the exponential volume dependence of the scattering length, 
the effective range etc., separately. 

It is important to stress the limits of validity of the present calculation. On one hand, 
the pion masses should be small enough so that the chiral expansion is converging. From the 
experience acquired in the three flavor case, where kaon loops are a borderline case for the 
convergence of the expansion, one expects chiral perturbation theory to be useful for m-,^ < 
500 MeV (of course the exponential volume dependence for a 500 MeV pion, or kaon will be 
negligible). Also, the box size has to be large enough so the usual power counting used here (the 
so-called "p-counting" ) is valid. When L is much smaller than the inverse pion mass, another 
power counting is required such as the e- |250| or the e'-regime [251 j . Additionally, we have 
neglected corrections which occur from higher loops, all of which are suppressed by additional 
factors of {m.,^ / Att f-jr)'^ and some of which are suppressed by additional exponential factors of 
e"™""^. The diagrams with this extra exponential suppression result from two-loop diagrams 
where intermediate states in both loops are purely off-shell and hence "going around the box" . 

We have focussed on the exponential corrections to phase shifts close to threshold. One can 
easily extend this work to include the exponential volume dependence of the phase shifts at 
higher energies. Alternatively, one can access non zero momenta by using twisted or partially 
twisted boundary conditions to probe the low-momentum dependence of the scattering ampli- 
tude ESlllMlIlSlIlSSllinillSIlllSHlIlSn]- TWs method boosts the entire two-particle system, 
however, requiring the extraction of scattering parameters in a boosted frame |4H {44, 260 . 
Our methods generalize trivially to this case. These methods can also be extended to other 
interesting two-hadron systems |43| I42j , where these exponential volume effects may be more 
significant. 



shown that the these two lattice artifacts were largely suppressed, and almost non-existent for the mass tunings 
used in Ref. |47|. However, as shown in Ref. |137| . there are still partial quenching effects and in particular, 
in the t- and u-channel diagrams the hairpin contributions can be significantly more sensitive to the boundary 
effects. For 1 — 2 these effects are only exponential, and for the pion masses and box sizes used in Ref. (47| . 
we have found they are the same order of magnitude as the corrections of this paper, and thus not-significant 
to the work of Beane et. al. 
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Chapter 6 

TWO-HADRON INTERACTIONS: 
GINSPARG- WILSON PIONS SCATTERING IN A SEA OF STAGGERED 

QUARKS 

The work in this Chapter is based upon Ref. jl37j . 
6.1 Introduction 

Lattice simulations with staggered fermions [IH] can at present reach significantly lighter quark 
masses than other fermion discretizations^ and have proven extremely successful in accurately 
reproducing experimentally measurable quantities [2651 126b'] . Staggered fermions, however, have 
the disadvantage that each quark flavor comes in four tastes. Because these species are degen- 
erate in the continuum, one can formally remove them by taking the fourth root of the quark 
determinant. In practice, however, the fourth root must be taken before the continuum limit; 
thus it is an open theoretical question whether or not this fourth-rooted theory becomes QCD 
in the continuum limit. ^ Even if one assumes the validity of the fourth-root trick, which we do 
in the rest of this chapter, staggered fermions have other drawbacks. On the lattice, the four 
tastes of each quark flavor are no longer degenerate, and this taste symmetry breaking is numer- 
ically significant in current simulations |266j . Thus one must use staggered chiral perturbation 
theory (SxPT), which accounts for taste-breaking discretization effects, to extrapolate correctly 
staggered lattice calculations to the continuum |275| 1276^ I277| l278] . Fits of S^PT expressions 
for meson masses and decay constants have been remarkably successful. Nevertheless, the large 
number of operators in the next-to-leading order (NLO) staggered chiral Lagrangian [278J and 
the complicated form of the kaon B-parameter in S^PT 279 both show that SxPT expressions 
for many physical quantities will contain a daunting number of undetermined fit parameters. 
Another practical hindrance to the use of staggered fermions as valence quarks is the construc- 
tion of lattice interpolating fields. Although the construction of a staggered interpolating field 
is straightforward for mesons since they are spin objects [28n[ 1281] . this is not in general the 
case for vector mesons, baryons or multi-hadron states since the lattice rotation operators mix 
the spin, angular momentum and taste of a given interpolating field |282[ I283[ 1284] . 

The use of Ginsparg- Wilson (GW) fermions [55'| evades both the practical and theoretical 
issues associated with staggered fermions. Because GW fermions are tasteless, one can simply 
construct interpolating operators with the right quantum numbers for the desired meson or 



We note that recently there have been significant improvements in the last few years for developing numeri- 
cally fast Wilson fermions ,261 262 263 264 . in addition to the twisted mass discretization method discussed 
briefly in Chapter |21 

^See Ref. |^ 1^ 1^ l?fCTl ETTl l?f^ 1^ for a recent discussion of staggered fermions and the fourth- 
root trick. 
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baryon. Moreover, massless GW fermions possess an exact chiral symmetry on the lattice |285j 
which protects expressions in xPT from becoming unwieldy.^ Unfortunately, simulations with 
dynamical GW quarks are approximately 10 to 100 times slower than those with staggered 
quarks |161j and thus are not presently practical for realizing light quark masses. 

A practical compromise is therefore the use of GW valence quarks and staggered sea quarks. 
This so-called "mixed action" theory is particularly appealing because the MILC improved 
staggered field configurations are publicly available. Thus one only needs to calculate correlation 
functions on top of these background configurations, making the numerical cost the same as that 
of quenched GW simulations. Several lattice calculations using domain- wall or overlap valence 
quarks with the MILC configurations are underway t292l 12951 including a determination 
of the isospin 2 (/ = 2) vrvr scattering length 07]. Although this is not the first 1 = 2 vrvr 
scattering lattice simulation j236| I237| I242| I244| I246j , it is the only one with pions light enough 
to be in the chiral regime [511 152j . Its precision is limited, however, without the appropriate 
mixed action xPT expression for use in continuum and chiral extrapolation of the lattice data. 
With this motivation we calculate the 1 = 2 vrvr scattering length in chiral perturbation theory 
for a mixed action theory with GW valence quarks and staggered sea quarks. 

Mixed action chiral perturbation theory (MA^PT) was first introduced in Refs. ISHU ITTH 
I178j and was extended to include GW valence quarks on staggered sea quarks for both mesons 
and baryons in Refs. j249j and jl93j . respectively, vrvr scattering is well understood in con- 
tinuum, infinite- volume xPT [Ml EB ED EH OHSl EHHl EHZl, and is the simplest two-hadron 
process that one can study numerically with LQCD. We extend the NLO calculations of 
Refs. |7()( I71j to MA^PT. A mixed action simulation necessarily involves partially quenched 
QCD (PQQCD) [UHlIIiniinZllMlIIlHlinil, in which the valence and sea quarks are treated 
differently. Consequently, we provide the PQxPT tttt scattering amplitude by taking an appro- 
priate limit of our MAxPT expressions. In all of our computations, we work in the isospin limit 
both in the sea and valence sectors. 

This chapter is organized as follows. We first comment on the determination of infinite 
volume scattering parameters from lattice simulations in Section [6.21 focusing on the applica- 
bility of Liischer's method |38| I39j to mixed action lattice simulations. We then review mixed 
action LQCD and MA^PT in Section [6.31 In Section [6.41 we calculate the 1 = 2 tttt scattering 
amplitude in MAxPT, first by reviewing vrvr scattering in continuum SU{2) xPT and then 
by extending to partially quenched mixed action theories with Nf = 2 and Nf = 2 -|- 1 sea 
quarks. We discuss the role of the double poles in this process |298j and parameterize the 
partial quenching effects in a particularly useful way for taking various interesting and impor- 
tant limits. Next, in section 16.51 we present results for the pion scattering length in both 2 
and 2 + 1 flavor MA^PT. These expressions show that it is advantageous to fit to partially 
quenched lattice data using the lattice pion mass and pion decay constant measured on the 
lattice rather than the LO parameters in the chiral Lagrangian. We also give expressions for 
the corresponding continuum PQxPT scattering amplitudes, which do not already appear in 
the literature. Finally, in Section 16.61 we briefly discuss how to use our MAxPT formulae to 
determine the physical scattering length in QCD from mixed action lattice data and conclude. 



■^In practice, the degree of chiral symmetry is limited by how well the domain-wall fermion |286l 12871 1288| is 
realized or the overlap operator |289l I29(jl I291| is approximated. 
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6.2 Determination of Scattering Parameters from Mixed Action Lattice Simula- 
tions 

Lattice QCD calculations are performed in Euclidean spacetime, thereby precluding the ex- 
traction of S-matrix elements from infinite volume |40j. Liischer, however, developed a method 
to extract the scattering phase shifts of two particle scattering states in quantum field the- 
ory by studying the volume dependence of two-point correlation functions in Euclidean space- 
time [381 139j ■ In particular, for two particles of equal mass m in an s-wave state with zero total 
3-momentum in a finite volume, the difference between the energy of the two particles and twice 
their rest mass is related to the s-wave scattering length:^ 



1 + 



(6.1) 



In the above expression, oq is the scattering length (not to be confused with the lattice spacing, 
a), L is the length of one side of the spatially symmetric lattice, and ci and C2 are known 
geometric coefficients.^ Thus, even though one cannot directly calculate scattering amplitudes 
with lattice simulations, Eq. 1)6. If) , which we will refer to as Liischer's formula, allows one 
to determine the infinite volume scattering length. One can then use the expression for the 
scattering length computed in infinite volume xPT to extrapolate the lattice data to the physical 
quark masses. 

Because Liischer's method requires the extraction of energy levels, it relies upon the existence 
of a Hamiltonian for the theory being studied. This has not been demonstrated (and is likely 
false) for partially quenched and mixed action QCD, both of which are nonunitary. Nevertheless, 
one can calculate the ratio of the two-pion correlator to the square of the single-pion correlator 
in lattice simulations of these theories and extract the coefficient of the term which is linear 
in time, which becomes the energy shift in the QCD (and continuum) limit. We claim that 
in certain scattering channels, despite the inherent sicknesses of partially quenched and mixed 
action QCD, this quantity is still related to the infinite volume scattering length via Eq. (|6.H) . 
i.e. the volume dependence is identical to Eq. (|6.1() up to exponentially suppressed corrections.^ 
This is what we mean by "Liischer's method" for nonunitary theories. We will expand upon 
this point in the following paragraphs. 

It is well known that Liischer's formula does not hold for many scattering channels in 
quenched theories because unitarity-violating diagrams give rise to enhanced finite volume 
effects |299j . For certain scattering channels, however, quenched xPT calculations in finite 
volume show that, at 1-loop order, the volume dependence is identical in form to Liischer's for- 
mula |299| I3f)n[ 130 Ij . Chiral perturbation theory calculations additionally show that the same 
sicknesses that generate enhanced finite volume effects in quenched QCD also do so in partially 
quenched and mixed action theories ^,128, 134, 294'", TtT, "302 , 249, 303 . It then follows that 

■^Here we use the "particle physics" definition of the scattering length which is opposite in sign to the "nuclear 
physics" definition. 

^This expression generalizes to scattering parameters of higher partial waves and non-stationary particles |38l 

Emu EH. 

^Here, and in the following discussion, we restrict ourselves to a perturbative analysis. 
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if a given scattering channel has the same volume dependence as Eq. (|6.H) in quenched QCD, 
the corresponding partially quenched (and mixed action) two-particle process will also obey 
Eq. ()6.1() . Correspondingly, scattering channels which have enhanced volume dependence in 
quenched QCD also have enhanced volume dependence in partially quenched and mixed action 
theories. We now proceed to discuss in some detail why Liischer's formula does or does not 
hold for various 2— >2 scattering channels. 

Finite volume effects in lattice simulations come from the ability of particles to propagate 
over long distances and feel the finite extent of the box through boundary conditions. Generi- 
cally, they are proportional either to inverse powers of L or to exp(-mL), but Liischer's formula 
neglects exponentially suppressed corrections. Calculations of scattering processes in effective 
field theories at finite volume show that the power-law corrections only arise from s-channel 
diagrams |299[ 13041 13021 Il3] . This is because all of the intermediate particles can go 
on-shell simultaneously, and thus are most sensitive to boundary effects. Consequently, when 
there are no unitarity-violating effects in the s-channel diagrams for a particular scattering 
process, the volume dependence will be identical to Eq. 1)6. 1() . up to exponential corrections. 
Unitarity-violating hairpin propagators in s-channel diagrams, however, give rise to enhanced 
volume corrections because they contain double poles which are more sensitive to boundary 
effects |299j .^ Thus all violations of Liischer's formula come from on-shell hairpins in the s- 
channel. 

Let us now consider 1 = 2 vrvr scattering in the mixed action theory. All intermediate states 
must have isospin 2 and s > Attt?. If one cuts an arbitrary graph connecting the incoming and 
outgoing pions, there is only enough energy for two of the internal pions to be on-shell, and, by 
conservation of isospin, they must be valence vr+'s.^ Thus no hairpin diagrams ever go on-shell 
in the s-channel, and the structure of the integrals which contribute to the power-law volume 
dependence in the partially quenched and mixed action theories is identical to that in continuum 
xPT. This insures that Liischer's formula is correctly reproduced to all orders in 1/L with the 
correct ratios between coefficients of the various terms. Moreover, this holds to all orders in xPT, 
PQxPT, MA^PT, and even quenched xPT- The sicknesses of the partially quenched and mixed 
action theories only alter the exponential volume dependence of the 1 = 2 scattering amplitude.^ 
This is in contrast to the / = tttt amplitude, which suffers from enhanced volume corrections 
away from the QCD limit. In general, the argument which protects Liischer's formula from 
enhanced power-like volume corrections holds for all "maximally-stretched" states at threshold 
in the meson sector, i.e. those with the maximal values of all conserved quantum numbers; 
other examples include K~^K~^ and K~^7r~^ scattering. We expect that a similar argument will 



^We note that, while the enhanced volume corrections in quenched QCD invalidate the extraction of scattering 
parameters from certain scattering channels, e.g. 1 = 299 301 , this is not the case in principle for partially 
quenched QCD, since QCD is a subset of the theory. Because the enhanced volume contributions must vanish in 
the QCD limit, they provide a "handle" on the enhanced volume terms. In practice, however, these enhanced 
volume terms may dominate the correlation function, making the extraction of the desired (non-enhanced) 
volume dependence impractical. 

*We restrict the incoming pions to be below the inelastic threshold; this is necessary for the validity of 
Liischer's formula even in QCD. 

^In fact, hairpin propagators will give larger exponential dependence than standard propagators because they 
are more chirally sensitive. 
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hold for certain scattering channels in the baryon sector. 

Therefore the s-wave 1 = 2 vrvr scattering length can be extracted from mixed action lattice 
simulations using Liischer's formula and then extrapolated to the physical quark masses and to 
the continuum using the infinite volume MAxPT expression for the scattering length.^'' 

6.3 Mixed Action Lagrangian and Partial Quenching 

Mixed action theories use different discretization techniques in the valence and sea sectors and 
are therefore a natural extension of partially quenched theories. We consider a theory with 
Nf staggered sea quarks and valence quarks (with corresponding ghost quarks) which 
satisfy the Ginsparg- Wilson relation j45l I285j . In particular we are interested in theories with 
two light dynamical quarks (A'^^ = 2) and with three dynamical quarks where the two light 
quarks are degenerate (commonly referred to as Nf = 2 + 1). To construct the continuum 
effective Lagrangian which includes lattice artifacts one follows the two step procedure outlined 
in Ref. jl66j . First one constructs the Symanzik continuum effective Lagrangian at the quark 
level |18H llSDj up to a given order in the lattice spacing, a: 

Csym = C + a£(5) + a2^(6) + . . . ^ (6.2) 

where contains higher dimensional operators of dimension 4 + n. Next one uses the 

method of spurion analysis to map the Symanzik action onto a chiral Lagrangian, in terms of 
pseudo-Goldstone mesons, which now incorporates the lattice spacing effects. This has been 
done in detail for a mixed GW-staggered theory in Ref. |249j : here we only describe the results. 
The leading quark level Lagrangian is given by 

4Nf+2Ny 
a,b=l 

where the quark fields are collected in the vectors 

g^/=2 = ( , ii , j2 , h ,kM,h,h, k , u,^ d) , (6.4) 

valence sea ghost 

QNf =2+1 ^ (n, d, g , ji , j2 ,h,j4,li, h ,h ,h,ri,r2,rs, n , u, d, sf (6.5) 

valence sea ghost 

for the two theories. There are 4 tastes for each flavor of sea quark, j, Z,r.^^We work in the 
isospin limit in both the valence and sea sectors so the quark mass matrix in the 2+1 sea flavor 
theory is given by 

mq = diag(mu, m^, ms, mj , mj , mj , mj , rrij , mj , rnj , ruj , rur, mr,mr,mr,mu, niu, mg). (6.6) 

valence sea ghost 

^"For a related discussion, see Ref. |255| 

^^Note that we use different labels for the valence and sea quarks than Ref. |249| . Instead we use the labeling 
convention consistent with Ref. |193| . 
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The quark mass matrix in the two flavor theory is analogous but without strange valence, sea and 
ghost quark masses. The leading order mixed action Lagrangian, Eq. (jHIHI), has an approximate 
graded chiral symmetry, SU{ANf + Ny\N^)L (g) SU{ANf + Ny\N^)R, which is exact in the 
massless hmit. In analogy to QCD, we assume that the vacuum spontaneously breaks this 
symmetry down to its vector subgroup, SU{4:Nf + A^„|A''^)y, giving rise to {ANj + 2N^Y ~ 1 
pseudo-Goldstone mesons. These mesons are contained in the field 



exp 



~T 



M 

X M 



(6.7) 



The matrices M and M contain bosonic mesons while x &iid contain fermionic mesons. 
Specifically, 
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(6.8) 



In Eq. 1)6. 8() we only explicitly show the mesons needed in the two flavor theory. The ellipses 
indicate mesons containing strange quarks in the 2+1 theory. The upper Ny x block of M 
contains the usual mesons composed of a valence quark and anti-quark. The fields composed 
of one valence quark and one sea anti-quark, such as (puj, are 1x4 matrices of fields where we 
have suppressed the taste index on the sea quarks. Likewise, the sea-sea mesons such as (pji are 
4x4 matrix- fields. Under chiral transformations, T, transforms as 



(6. 



In order to construct the chiral Lagrangian it is useful to first define a power-counting 
scheme. Continuum xPT is an expansion in powers of the pseudo-Goldstone meson momentum 
and mass squared [701 El ■ 

^pl/Kl^ml/Kl, (6.10) 

where oc mq and A-^ is the cutoff of xPT. In a mixed theory (or any theory which incorpo- 
rates lattice spacing artifacts) one must also include the lattice spacing in the power counting. 
Both the chiral symmetry of the Ginsparg- Wilson valence quarks and the remnant U{1)a sym- 
metry of the staggered sea quarks forbid operators of dimension five; therefore the leading lattice 



This is a "fake" symmetry of PQQCD. However, it gives the correct Ward identities and thus can be used 
to understand the symmetries and symmetry breaking of PQQCD 



118 



spacing correction for this mixed action theory arises at 0{a'^). Moreover, current staggered 
lattice simulations indicate that taste-breaking effects (which are of 0(0^)) are numerically of 
the same size as the lightest staggered meson mass |266| . We therefore adopt the following 
power-counting scheme: 

~ ~ ^q/^QCB ~ a^AqcD • (6-11) 

The leading order (LO), 0{e^), La grangian is then given in Minkowski space by j249j 

C = ^str (a^S a^St) + ^str (Sm^ + mgSt) - {Us +U's+ Uy) , (6.12) 

where we use the normalization / ~ 132 MeV and have already integrated out the taste singlet 
<I>o field, which is proportional to str(<I>) |128j . Us and Ug are the well-known taste breaking 
potential arising from the staggered sea quarks j275[ 1276] . The staggered potential only enters 
into our calculation through an additive shift to the sea-sea meson masses; we therefore do 
not write out its explicit form. The enhanced chiral properties of the mixed action theory are 
illustrated by the fact that only one new potential term arises at this order: 

Uv = -CMix str (TaSTsSt) , (6.13) 

where 

T^ = Vs-Vv = diag(-/v, h -Iv)- (6.14) 

The projectors, Vs and Vy, project onto the sea and valence-ghost sectors of the theory, ly 
and Is are the valence and sea flavor identities, and It is the taste identity matrix. From 
this Lagrangian, one can compute the LO masses of the various pseudo-Goldstone mesons in 
Eq. (|f).8j) . For mesons composed of only valence (ghost) quarks of flavors a and b, 

mlf^ = Binia + mb). (6.15) 

This is identical to the continuum LO meson mass because the chiral properties of Ginsparg- 
Wilson quarks protect mesons composed of only valence (ghost) quarks from receiving mass 
corrections proportional to the lattice spacing. However, mesons composed of only sea quarks 
of flavors si and S2 and taste t, or mixed mesons with one valence (v) and one sea (s) quark 
both receive lattice spacing mass shifts. Their LO masses are given by 

rnis,,t = B{m,, + m,,) + a'/^iit), (6.16) 
ml^ = B{m^ + nis) + a^Anjix- (6-17) 

From now on we use tildes to indicate masses that include lattice spacing shifts. The only sea- 
sea mesons that enter vrvr scattering to the order at which we are working are the taste-singlet 
mesons (this is because the valence- valence pions that are being scattered are tasteless), which 
are the heaviest; we therefore drop the taste label, t. The splittings between meson masses of 
different tastes have been determined numerically on the MILC configurations |2(i6j . so A(^/) 
should be considered an input rather than a fit parameter. The mixed mesons all receive the 
same shift given by 

IGCmiX Ir. 1 

AMix = J2 ' {^.\^) 
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which has yet to be determined numerically. 

After integrating out the <&o held, the two point correlation functions for the flavor-neutral 
states deviate from the simple single pole form. The momentum space propagator between two 
flavor neutral states is found to be at leading order |128j 

a (n^) - : ilkil(P^ ~ ™l ' -I In -,q^ 



- m^^ + ie Nf (jp, _ j„2^ _|_ ig)(p2 _ _|_ W^L^{p^ - rh\, + ie) 



where 

+ 1 for a = valence or sea quarks 



1 for a = ghost quarks . (6.20) 

In Eq. (|6.19l) . k runs over the flavor neutral states {(pjj, (pu, and k' runs over the mass 
eigenstates of the sea sector. For vrvr scattering, it will be useful to work with linear combinations 
of these ija fields. In particular we form the linear combinations 

for which the propagators are 

QAp^) = -o , (6.22) 



p2 _ j^2 _j_ ' 



p2 - m2 + ie Nf (p2 -rnl + ie)^ nfe'Li (P^ " + ie) 



Specifically, 



2 ~ 2 



2i {p^ - rnjj){p'^ - m^^) 



6.4 Calculation of the 1 = 2 Pion Scattering Amplitude 



for Nf = 2 + 1, (6.25) 



where m'i = 4(m^,- + 2m2j, 



Our goal in this work is to calculate the 1 = 2 vrvr scattering length in chiral perturbation 
theory for a partially quenched, mixed action theory with GW valence quarks and staggered 
sea quarks, in order to allow correct continuum and chiral extrapolation of mixed action lattice 
data. We begin, however, by reviewing the pion scattering amplitude in continuum SU (2) 
chiral perturbation theory. We next calculate the scattering amplitude in Nf = 2 PQxPT and 
MAxPT, and finally in Nf = 2 + 1 PQxPT and MA^PT. When renormalizing divergent 1-loop 
integrals, we use dimensional regularization and a modified minimal subtraction scheme {MS) 
where we consistently subtract all terms proportional to [70] : 

2 

- 7i? + log47r + 1, 

where d is the number of space-time dimensions. The scattering amplitude can be related to 
the scattering length and other scattering parameters, as we discuss in Section 16.51 
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6.4-1 Continuum SU (2) 

The tree-level 1 = 2 pion scattering amplitude at threshold is well known to be jlOOj 

Ami 



T 



f2 
J n 



(6.26) 



It is corrected at O(e^) by loop diagrams and also by tree level terms from the NLO (or Gasser- 
Leutwyler) chiral Lagrangian fT^P Thedia grams that contribute at one loop order are shown 
in Figure they lead to the following NLO expression for the scattering amplitude: 
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(6.27) 



where rriuu is the tree-level expression given in Eq. (|6.15|) and / is the LO pion decay constant 
which appears in Eq. 1)6. 7() . The coefficient l'^^ is a linear combination of low energy constants 
appearing in the Gasser-Leutwyler Lagrangian whose scale dependence exactly cancels the scale 
dependence of the logarithmic term. One can re-express the amplitude, however, in terms of 
the physical pion mass and decay constant using the NLO formulae for m,r and fj^ to find: 
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J TT 
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(6.28) 



where Ztttt is a different linear combination of low energy constants. The expression for It^tj 
can be found in Ref. |296j . We do not, however, include it here because we do not envision 
either using the known values of the Gasser-Leutwyler parameters in the the fit of the scattering 
length or using the fit to determine them. The simple expression 1)6. 28p has already been used in 
extrapolation of lattice data from mixed action simulations |47j , but it neglects lattice spacing 
effects from the staggered sea quarks which are known from other simulations to be of the 
same order as the leading order terms in the chiral expansion of some observables j266| . We 
therefore proceed to calculate the scattering amplitude in a partially quenched, mixed action 
theory relevant to simulations. 



6.4-2 Mixed GW-Staggered Theory with two Sea Quarks 

The scattering amplitude in the partially quenched theory differs from the unquenched theory 
in three important respects. First, more mesons propagate in the loop diagrams. Second, 
some of the mesons have more complicated propagators due to hairpin diagrams at the quark 
level 298, 128 . Third, there are additional terms in the NLO Lagrangian which arise from 
partial quenching jl32j , and lattice spacing effects |249l I278j . 

At the level of quark flow, there are diagrams such as Figure which route the valence 
quarks through the diagram in a way which has no ghostly counterpart. Consequently, the 
ghosts do not exactly cancel the valence quarks in loops. Of course, this is simply a reflection 
of the fact that the initial and final states — valence pions — are themselves not symmetric 

^^The continuum tttt scattering amplitude is known to two-loops |97l 12951 12961 1297| . 
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Figure 6.1: Example quark flow for a one-loop t-channel graph. This diagram illustrates the 
presence of meson loops composed of purely valence- valence mesons which are not canceled by 
valence-ghost loops. Different colors (shades of grey) represent different quark flavors. 



under the interchange of ghost and valence quarks, and therefore the graded symmetry between 
the valence and ghost pions has already been violated. This is well known in quenched and 
partially quenched heavy baryon xPT jl^-iSL II, SOL ll.Slj . This fact also partly explains the success 
of quenched tttt scattering in the 1 = 2 channel [2361 123?] : quenching does not eliminate all 
loop graphs like it does in many other processes, and in particular, the s-channel diagram is 
not modified by (partial) quenching effects. As a consequence, it is necessary to compute all 
the graphs contributing to this process in order to determine the scattering amplitude. 

Quark level disconnected (hairpin) diagrams lead to higher order poles in the propagator of 
any particle which has the quantum numbers of the vacuum 298 , 128]. In the isospin limit of 
the Nf = 2 partially quenched theory, conservation of isospin prevents the vr'^ from suffering 
any hairpin effects at leading order. At higher orders, the vr'' mass (as well as the other mesons) 
will receive hairpin contributions but the vr'' propagator will never develop double poles, the 
characteristic unitarity violating feature of partially quenched and mixed action theories. Hence 
only the fj acquires a disconnected propagator. Moreover, in the mo oo limit, the fj propagator 
(given for a general PQ theory in Eq. 1)6. 23(1 ') is given by the simple expression 



Gfiip 



(6.29) 



''^PQ 



2\2 ■ 



where the parameter 



^Iq = rh^j - ml (6.30) 

quantifies the partial quenching. (Recall that rhjj is the physical mass of a taste singlet sea-sea 
meson.) Notice that when Apg — > the propagator ()6.29|) also goes to zero; this is what we 
expect since, in the SU{2) theory, the only neutral propagating state is the vr". The propagator 
in Eq. 1)6. 29() can appear in loops, thereby producing new diagrams such as those in Fig. 16.21 ^^ 

^''We note that there are also similar contributions to the four particle vertex with a loop and to the mass 
correction. We do not show them, however, because they cancel against one another in the amplitude expressed 
in lattice-physical parameters, which we will show in the following pages. 



122 





Figure 6.2: Example hairpin diagrams contributing to pion scattering. The propagator with a 
cross through it indicates the quark-disconnected piece of the f] propagator, Eq. ()6.29|1 . 



After adding ah such hairpin diagrams, one finds that the contribution of the r] to the amphtude 
is 15 

4 ^PQ 



In addition to 1-loop contributions, the NLO scattering amphtude receives tree-level ana- 
lytic contributions from operators of O(e^) in the chiral Lagrangian. At this order, the mixed 
action Lagrangian contains the same O(p^), 0{p'^mg), and 0{in'^) operators as in the continuum 
partially quenched chiral Lagrangian, plus additional 0{a'^), 0{a^mq), and 0{a?p'^) operators 
arising from discretization effects. We can now enumerate the generic forms of analytic contri- 
butions from these NLO operators. Because of the chiral symmetry of the GW valence sector, 
all tree-level contributions to the scattering length must vanish in the limit of vanishing valence 
quark mass.^^ Thus there are only three possible forms, each of which must be multiplied by 
an undetermined coefficient: m^^, m'^^m'jj, and m'^^a'^. It may, at first, seem surprising that 
operators of O{o?mq), which come from taste-symmetry breaking and contain projectors onto 
the sea sector, can contribute at tree- level to a purely valence quantity. Nevertheless, this turns 
out to be the case. These 0{a^mq) mixed action operators can be determined by first starting 
with the NLO staggered chiral Lagrangian |278j . and then inserting a sea projector, Vs, next to 

every taste matrix. One example of such an operator is str ^Sm^^ str (T-'s^sS^sS"'') -|- p.c. , 
where, ^5 is the 75 matrix acting in taste-space and p.c. indicates parity-conjugate. This double- 
trace operator will contribute to the lattice pion mass, decay constant, and 4-point function 
at tree-level because one can place all of the valence pions inside the first supertrace, and the 
second supertrace containing the projector Vs will just reduce to the identity. 



^^We note that this contribution does not vanish in the Umit that — > with m^j 7^ 0. Similar effects 
have been observed in quenched computations of pion scattering amphtudes |300l 1295) . This non- vanishing 
contribution is the 7 = 2 remnant of the divergences that are known to occur in the / = amplitude at 
threshold. These divergences give rise to enhanced volume corrections to the I = Q amplitude with respect to 
the one-loop 1 — 2 amplitude and prevent the use of Liischer's formula. Moreover, it is known [1261 1531 that 
PQxPT is singular in the limit rriu —> with nonzero sea quark masses, so the behavior of the amplitude in 
this limit is meaningless. 

^® As we discussed in the previous footnote, this condition need not hold for loop contributions to the scattering 
amplitude. 
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Putting everything together, the total mixed action scattering amplitude to NLO is 
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The first fine of Eq. (|6.32() contains those terms which remain in the continuum and full QCD 
hmit, Eq. ()6.27|) . while the second and third lines account for the effects of partial quenching 
and of the nonzero lattice spacing. Note that, for consistency with the 1-loop terms, we chose 
to re-express the analytic contribution proportional to the sea quark mass as m^^Apg. In 
Eq. ()6.32|) we have multiplied every contribution from diagrams which contain a sea quark loop 
by 1/4, thus making our expression applicable to lattice simulations in which the fourth root 
of the staggered sea quark determinant is taken. 

It is useful, however, to re-express the scattering amplitude in terms of the quantities that 
one measures in a lattice simulation: and fj^. Throughout this chapter, we will refer to 
these renormalized measured quantities as the lattice-physical pion mass and decay constant. 
Because we are working consistently to second order in chiral perturbation theory, we can equate 
the lattice-physical pion mass to the 1-loop chiral perturbation theory expression for the pion 
mass, and likewise for the lattice-physical decay constant. Thus, in terms of lattice-physical 
parameters, the mixed action 1 = 2 vrvr scattering amplitude is 
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(6.33) 



where the first few terms are identical in form to the full QCD amplitude, Eq. ()6.28|) . This 
expression for the scattering amplitude is vastly simpler than the one in terms of the bare 
parameters. First, the hairpin contributions from all diagrams except those in Fig. 16.21 have 
exactly cancelled, removing the enhanced chiral logs and leaving the last term in Eq. (|6.33j) 
as the only explicit modification arising from the partial quenching and discretization effects. 
Second, all contributions from mixed valence-sea mesons in loops have cancelled, thereby re- 
moving the new mixed action parameter, Cmixi completely. Third, all tree-level contributions 
proportional to the sea quark mass have also cancelled from this expression. And finally, most 
striking is the fact that an explicit computation of the 0{a?mq) contributions to the amplitude 
arising from the NLO mixed action Lagrangian show that these effects exactly cancel when the 



Notice that once the lattice spacing a has been determined, the lattice-physical pion mass can be unambigu- 
ously determined by measuring the exponential decay of a pion-pion correlator. We assume that the lattice 
spacing a has been determined, for example, by studying the heavy quark potential or quarkonium spectrum. 

Another consequence of the exact cancellation of the loops with mixed valence-sea quarks is that one does 
not have to implement the "fourth-root trick" through this order. 
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amplitude is expressed in lattice-physical parameters. This result will be discussed in detail 
in Ref. Thus to reiterate, the only partial quenching and lattice spacing dependence 

in the amplitude comes from the hairpin diagrams of Fig. 16.21 which produce contributions 
proportional to ApQ = {rrijj + a'^A{^j) — m^)^, where m^^- + a'^A{^i) is the mass-squared of 
the taste-singlet sea-sea meson. Moreover, we presume that anyone performing a mixed action 
lattice simulation will separately measure the taste-singlet sea-sea meson mass and use it as an 
input to fits of other quantities such as the vrvr scattering length. Thus we do not consider it to 
be an undetermined parameter. 

It is now clear that one should fit tttt scattering lattice data in terms of the lattice-physical 
pion mass and decay constant rather than in terms of the LO pion mass and LO decay constant. 
By doing this, one eliminates three undetermined fit parameters: Cmix, I'pq, and l'^2, as well as 
the enhanced chiral logs. 



6.4-3 Mixed GW-Staggered Theory with 2 + 1 Sea Quarks 

The 2 -|- 1 flavor theory has three additional quarks - the strange valence and ghost and strange 
sea quarks - which can lead to new contributions to the scattering amplitude. Because we 
only consider the scattering of valence pions, however, strange valence quarks cannot appear 
in this process. Thus all new contributions to the scattering amplitude necessarily come only 
from the sea strange quark, r. Because the r quark is heavier than the other sea quarks 
there is SU{3) symmetry breaking in the sea. This symmetry breaking only affects the pion 
scattering amplitude, expressed in lattice-physical quantities, through the graphs with internal 
fj propagators because the masses of the mixed valence-sea mesons cancel in the final amplitude 
as they did in the earlier two flavor case. In addition, the only signature of partial quenching 
in the amplitude comes from these same diagrams. It is therefore worthwhile to investigate the 
physics of the neutral meson propagators further. 

There are more hairpin graphs in the 2-1-1 flavor theory since the rjg may propagate as well as 
the rju and the r/^. Because these mesons mix with one another, the flavor basis is not the most 
convenient basis for the computation. Rather, a useful basis of states is vr*^, fj = {rju + r\i)l\pl 
and r\s. Since we work in the isospin limit, the 'k^ cannot mix with f\ or t^^; in addition, there 
is no vertex between the r]a and tt^tt' at this order, so we never encounter a propagating r]s- 
Thus all the PQ effects are absorbed into the f] propagator, which is given by 



^fi^P ) ~ l2 Z3 ~ ^772 _2\2 ^„2 (QM) 



In SU{3) chiral perturbation theory, the neutral mesons are the vr'^ and the t]s- Therefore, in 
the PQ theory, we know that there will be a contribution from the fj graphs that does not result 
from partial quenching or SU{3) symmetry breaking. Therefore the extra PQ graphs arising 
from the internal 77 fields must not vanish in the Apg — > limit, in contrast to the two flavor 
case of Eq. (|6.31|1 . 

To make this clear, we can re-express the propagator of Eq. 1)6. 34() in terms of Apg as 
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(6.35) 
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This propagator has a single pole which is independent of Apq, as well as higher order poles 
that are at least quadratic in Apg. It is interesting to consider the large nir limit of this 
propagator. In this limit, ~ ^Brrir is also large. For momenta that are small compared to 



the second term of this equation goes to zero in the large m,. limit, and the 2 + 1 flavor 
propagator reduces to the 2 flavor propagator, Eq. H6.29() . as expected. 

While the above expression clarifies the Apg dependence of the propagator and the large 
rrir limit, it obscures the SU{3)sea limit. An equivalent form of the propagator is 
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where the quantity Aq 



m„ 



parametrizes the SU{3)sea breaking. When A3 = this 



propagator is similar in form to the corresponding 2 flavor propagator, Eq. ()6.29|1 . but it has 
an additional single pole due to the extra neutral meson in the SU{3) theory. 

Having considered the new physics of the hairpin propagator, we can now calculate the 
scattering amplitude. For our purposes here, it is most convenient to express the total 1 = 2 vrvr 
scattering amplitude in terms of Apg. Just as in the 2-flavor computation, the NLO analytic 
contributions due to partial quenching and finite lattice spacing effects exactly cancel when the 
amplitude is expressed in lattice-physical parameters. All sea quark mass and lattice spacing 
dependence comes from the hairpin diagrams, which produce terms proportional to powers of 
ApQ with known coefficients. The amplitude is 



MAxPT 



Pi=0 



Ami 



1 + 



mt 



(4vr^)2 
+ 



, mi\ 1 



1 



(4vrA)^ 



n=l ^ 




(6.37) 



where Apg 
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(6.38a) 
(6.38b) 
(6.38c) 
(6.38d) 



The functions J-'i have the property that J-i{x) ^ in the limit that x ^ 00. Therefore, when 
the strange sea quark mass is very large, i.e. m^/ml ^ 1, the 2 + 1 flavor amplitude reduces 
to the 2 flavor amplitude, Eq. (|6.33j) . with the exception of terms that can be absorbed into 
the analytic terms. The low energy constants have a scale dependence which exactly cancels 
the scale dependence in the logarithms. The coefficient ^tttt is the same linear combination of 
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Gasser-Leutwyler coefficients tliat appear in the SU{3) scattering amplitude expressed in terms 
of the physical pion mass and decay constant [HZl I297j . 

Because the functions J-i depend logarithmically on x, the 2 + 1 flavor scattering amplitude 
features enhanced chiral logarithms |126j that are absent from the 2 flavor amplitude. This is 
a useful observation, as we will now explain. Because there is a strange quark in nature and its 
mass is less than the QCD scale, Aqqd, lattice simulations must use 2 + 1 quark flavors. It is 
often practical to fix the strange quark mass at a constant value near its physical value in these 
simulations. This circumstance is helpful because, just as SU (2) chiral perturbation theory is 
useful to describe nature at scales smaller than the strange quark mass, the 2 flavor amplitude 
given in Eq. (|6.33|1 can be used to extrapolate 2 + 1 flavor lattice data at energy scales smaller 
than the strange sea quark mass used in the simulation (provided, of course, there are no strange 
valence quarks) |3n6j . This is valid because, at energy scales smaller than the strange quark 
mass (or actually twice the strange quark mass, since the purely pionic systems have no valence 
strange quarks), one can integrate out the strange quark. This is not an approximation, because 
all of the effects of the strange quark are absorbed into a renormalization of the parameters of 
the chiral Lagrangian. Moreover, since the 2 flavor amplitude does not exhibit enhanced chiral 
logarithms, signatures of partial quenching can be reduced by extrapolating lattice data with 
the 2 flavor, rather than the 2 + 1 flavor, expression. We note that in this case the effects of 
the strange quark are absorbed in the coefficients of the analytic terms appearing in Eq. (|6.33jl . 
and thus they are not constant, but rather depend logarithmically upon the strange sea quark 
mass. 



6.5 1 = 2 Pion Scattering Length Results 

In this section we present our results for the s-wave 1 = 2 tttt scattering length in the two 
theories most relevant to current mixed action lattice simulations: those with GW valence 
quarks and either Nf = 2 or Nf = 2 + 1 staggered sea quarks. We only present results for 
the scattering length expressed in lattice-physical parameters. The s-wave scattering length is 
trivially related to the full scattering amplitude at threshold by an overall prefactor: 
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6.5.1 Scattering Length with 2 Sea Quarks 

The 1 = 2 TTTT s-wave scattering length in a MA^PT theory with 2 sea quarks is given by 

3 In (4) - 1 + - TTW^I' (6-40) 

V / J (4vr/^)^6m4J 

where Apg = m|j + a^A(^/) — m^. The first two terms are the result one obtains in SU{2) 
xPT |296j and the last term is the only new effect arising from the partial quenching and mixed 
action. All other possible partial quenching terms, enhanced chiral logs and additional linear 
combinations of the 0{p^) Gasser-Leutwyler coefficients, exactly cancel when the scattering 
length is expressed in terms of lattice-physical parameters. And, most strikingly, the pion 
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mass, decay constant and the 4-point function all receive 0{a^mq) corrections from the lattice, 
but they exactly cancel in the scattering length expressed in terms of the lattice -physical 
parameters |3U5j . It is remarkable that the only artifact of the nonzero lattice spacing, m^- + 
a^A/, can be separately determined simply by measuring the exponential fall-off of the taste- 
singlet sea-sea meson 2-point function. Thus there are no undetermined fit parameters in the 
mixed action scattering length expression from either partial quenching or lattice discretization 
effects; there is only the unknown continuum coefficient, 1-,^^^. 

One can trivially deduce the continuum PQ scattering length from Eq. (|6.4nj) : simply let 
a ^ 0, reducing rhjj rrijj = 2Bmj in Apg, resulting in 
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6.5.2 Scattering Length with 2+1 Sea Quarks 

The 1 = 2 vrvr s-wave scattering length in a MAxPT theory with 2+1 sea quarks is given by 
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where the functions J-'i are defined in Eq. 1)6. 38() . As in the 2-flavor MA^PT expression, 
Eq. ()6.4U|1 . the only undetermined parameter is the linear combination of Gasser-Leutwyler 
coefficients, Ztttd which also appears in the continuum expression. 

We note as an aside that this suppression of lattice spacing counterterms is in contrast to the 
larger number of terms that one would need in order to correctly fit data from simulations with 
Wilson valence quarks on Wilson sea quarks. Because the Wilson action breaks chiral symmetry 
at 0{a), even for massless quarks, there will be terms proportional to all powers of the lattice 
spacing in the expression for the scattering length in Wilson xPT |3()5| 117(11 13()7j . Moreover, 
such lattice spacing corrections begin at 0{a), rather than 0{a'^). If one uses 0{a) improved 
Wilson quarks, then the leading discretization effects are of ©(a^), as for staggered quarks; 
however, this does not remove the additional chiral symmetry-breaking operators. Another 
practical issue is whether or not one can perform simulations with Wilson sea quarks that are 
light enough to be in the chiral regime. 



6.6 Discussion 



Considerable progress has recently been made in fully dynamical simulations of pion scattering 
in the 1 = 2 channel (2441 147j . We have considered 1 = 2 scattering of pions composed of 
Ginsparg- Wilson quarks on a staggered sea. We have calculated the scattering length in both 
this mixed action theory and in continuum PQxPT for theories with either 2 or 2 + 1 dynamical 
quarks. These expressions are necessary for the correct continuum and chiral extrapolation of 
PQ and mixed action lattice data to the physical pion mass. 
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Our formulae, Eqs. (|6.4n|) . (|6.42|) . not only provide the form for the mixed action scatter- 
ing length, but also contain two predictions relevant to the recent work of Ref. Beane 
et. al. calculated the 1 = 2 s— wave vrvr scattering length using domain wall valence quarks and 
staggered sea quarks, but used the continuum xPT expression to extrapolate to the physical 
quark masses. In Figure 2 of Ref. (see Fig. 16. 3j) . which plots m-n-a'^^ versus mj^/ f^^, the 
fit of the xPT expression to the lattice data does a remarkably good job. However, Eq. 1)6. 4U() 
predicts a known, positive shift to mTrOg^^ of size ApQ/(768/^7r^). Accounting for this positive 
shift is equivalent to lowering the entire curve. In Ref. |17], the valence and sea quark masses 
are tuned to be equal, so — cl^ 

A/ ~ (446 MeV)2 [SMl- Despite the large value of Apg, 
the predicted shift is insignificant, being an order of magnitude less than the statistical error. 
In Table Hm we collect the predicted shifts to m-^a^^ at the three pion masses used in Ref. [17| . 
We also list the magnitude of the ratio of these predicted shifts to the leading contribution to 
the scattering length, which turn out to be small, lending confidence to the power counting 
we have used, Eq. ()6.11|) . The other more important prediction is that there are no unknown 
corrections to the formula for the scattering length arising from lattice spacing corrections 
or partial quenching through the order 0{m'^), 0{a?"mq) and 0{a^). Therefore, to within sta- 
tistical and systematic errors, the continuum expression used by Beane et. al. to fit their 
numerical vrvr scattering data Jl?! receives no corrections through the 1-loop level. 

The central result of this chapter is that the appropriate way to extrapolate lattice vrvr 
scattering data is in terms of the lattice-physical pion mass and decay constant rather than in 
terms of the LO parameters which appear in the chiral Lagrangian. When expressed in terms 
of the LO parameters, the scattering length depends upon 4 undetermined parameters, Z^^, 
I'pQ, I'^i and Cmix- In contrast, the scattering length expressed in terms of the lattice-physical 
parameters depends upon only one unknown parameter, 1-,^^^, the same linear combination of 
Gasser-Leutwyler coefficients which contributes to the scattering length in continuum xPT. 
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Table 6.1: Predicted shifts to the scattering length computed in Ref. [17j arising from finite 
lattice spacing effects in the mixed action theory. The first two rows show the approximate 
values of and /jr while the third shows plus the statistical error calculated in In 

the fourth row, we give the predicted shifts in the scattering length (times m-,^) and, in the fifth 
row, we give the ratio of the predicted shift to the leading order contribution to the scattering 
length. 
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Figure 6.3: This is an updated plot of the 1 = 2 vrvr scattering length determined by the 
NPLQCD collaboration [37j. The scattering length is plotted as a function of m-,^/ as mea- 
sured in the lattice simulation, or in terms of the lattice parameters. The dashed line is the 
tree level prediction. The experimental point was not included in the fit. 
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Chapter 7 

MIXED ACTION LATTICE QCD WITH GINSPARG-WILSON 

VALENCE QUARKS 

This chapter is based upon work that is in progress |3U51 l3U8j . and which provides an 
understanding of the results found in the previous chapter. This chapter makes extensive use of 
Refs. [SI II2H1 HSOl CSIl ini Cni EB, ^93, 219, 278, ,2941 ISni- We discuss in some generahty 
the features of a mixed action theory with Ginsparg- Wilson valence quarks, and arbitrary sea 
quarks. We wih use staggered sea quarks at times for specific examples, as this combination of 
valence and sea quarks is the most relevant for current lattice simulations. 

7.1 Mixed Action Quark Level Lagrangian 

An understanding of mixed action lattice QCD must involve the effective theory for the mixed 
action simulation, but to construct this effective theory we must first understand the mixed 
action quark-level Lagrangian. As we have discussed several times in this work, the proce- 
dure is to first construct the effective continuum theory of the lattice action, the Symanzik 
Lagrangian |1811 ll8Uj , and then to construct the low-energy effective field theory which is con- 
sistent with the lattice symmetries jlHHj . The Symanzik Lagrangian for a mixed action theory 
has the form |171j 

C = CpQQCD + aC^''^ + a2£(6) + . . . (7.I) 

We see that the leading Lagrangian is simply the partially quenched Lagrangian which contains 
QCD as a limit. We comment again that a mixed action theory is automatically partially 
quenched, even if the bare quark mass parameters in the valence and sea sector are set equal. 
This is simply a consequence of the different lattice spacing modifications from the sea and 
valence Dirac operators, which force one to break the partially quenched symmetry and treat 
these quarks differently. However, we will show that it is useful to think of a mixed action 
theory with Ginsparg- Wilson valence quarks not as a mixed action theory, but rather as a full 
QCD or partially quenched theory of Ginsparg- Wilson quarks with a perturbative breaking 
of the symmetry between the valence and sea sectors, which can be seen by Eq. (|7.H) . To 
determine the mixed action effects, we must understand the operators in C^^'^ and beyond. We 
will in fact begin with C^^\ as typically, Wilson sea quarks are generally implemented with 0{a) 
improvement and staggered quark discretization effects begin at 0{a'^). To construct C^^\ it 
is convenient to group the different operators into a few categories. Firstly, there are those 
operators which consist of only sea-quarks, there are those of only valence quarks and then 
there are operators which mix the valence and sea quarks. We can further decompose these 
three sets into operators which break chiral symmetry, operators which do not break chiral 
symmetry and those operators which respect the lattice symmetry but break the continuum 
Lorentz symmetry. The mixed action symmetry forbids quark bi-linear operators which mix 
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valence and sea quarks and therefore the only type of operators which can mix the quarks are 
products of bi-linears. The Lorentz breaking operators are of the form^ 

QsJ^.D^Df^D'^Qs , Qyj,,D^Df'Df'Qv, (7.2) 

in which there is a sum over the Lorentz indices. In fact these operators can be distinguished in 
the dispersion relation for the various hadrons. For the mesons, these effects will be of 0{p'^a'^) 
and thus beyond the order we are interested, but for the baryons, these operators can enter at 
C(fl^) |lfl3j . and for a zero momentum projection they act as a simple mass shift. 

There are additional quark bi-linears which do not break chiral symmetry, and are of the 
form 

QsipD^Qs , QviPD^Qv- (7.3) 

These operators provide some difficulty from an effective field theory point of view. Under chiral 
transformations, these operators transform identically to the (PQ)QCD quark kinetic operator, 
and thus there is no way to distinguish these operators by their symmetry transformations. 
Thus, these operators lead to an O(a^) shift to all the LECs of the continuum chiral Lagrangian. 
For example, these operators will contribute a perturbative correction to the chiral condensate 
which drives chiral symmetry breaking |171j . Thus to disentangle these effects from the physics 
of interest (the continuum LECs), one must vary the lattice spacing in combination with the 
quark masses. However, unlike with the continuum LECs, we do not necessarily need or want to 
know the precise values of these LECs associated with finite lattice spacing, as they are lattice 
artifacts and not physical. Rather we want to be able parameterize their total contributions to 
a given correlation function or observable, such that there effects can be removed and we can 
extract the physics of interest. 

From a mixed action perspective, the most interesting operators are the 4-quark operators 
which allow an interaction of the valence and sea quarks through products of quark bi-linears. 
These 4-quark operators will have the general form 

{Qs^a,bQs){Qs^a,bQs) 1 {Qv^a,bQv) {Qv^a,bQv) , 

{QsrA,BQs){Qv'rA,BQv) , (7.4) 

where the matrix Ta_,b acts on the internal quantum numbers of the quark fields; the spin and 
for example in the case of staggered quarks, the taste. Perhaps a more illuminating form to 
write these operators is with the valence and sea projectors, 

Qs'^A,BQs = Q'PsrA,BrsQ 

Qv^a,bQv = QVvTa,bVvQ. (7.5) 

For the unmixed operators, one has all the possibilities discussed in Ref. |171j . The Ginsparg- 
Wilson symmetry of the valence quarks prevents the addition of any operator which breaks chiral 
symmetry except for operators involving the quark mass matrix. Thus, for the dimension-6 
operators with Ginsparg- Wilson valence quarks, there are only two mixed action operators, 

^Mix = C (Ql^Vv Q) (Qi^Vs Q) + C {Q^^j^Vv Q) {Qi^i^Vs Q) ■ (7.6) 



^In Euclidean spacetime, as with a lattice simulation, it is an 0(4) rotational symmetry which is broken. 
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These operators give rise to the mixed term in Eq. (|6.13|) for the example of staggered sea quarks, 
and in the following we will be interested in the next higher order effective terms arising from 
these mixed quark level operators. 

Why go through all this effort to construct the Lagrangian which includes all these lattice 
spacing corrections? Why not simply treat all the continuum LECs as polynomial functions of 
the lattice spacing? There are several answers to this question. Firstly, not all of the lattice 
spacing dependence can be captured in this way. For example, with a mixed action theory, there 
will be hairpin contributions such as those in Eq. (|6.31|) which can not be parameterized with the 
continuum LECs. For lattice actions which give rise to a splitting amongst hadron masses which 
are degenerate in the continuum limit, such as twisted mass LQCD or staggered LQCD, these 
mass splittings can only be understood with the appropriate effective theories jl58l 11751 1276j , 
and again can not be parameterized with a simple polynomial dependence of the continuum 
LECs. More importantly, with lattice simulations to date and for the foreseeable future, these 
lattice spacing contributions are not negligible and therefore there effects must be understood. 
Additionally, as a great example, the form of the 1 = 2 vrvr scattering length worked out 
in the previous chapter shows that there was a remarkable simplification of the scattering 
length expressed in terms of the lattice pion mass and lattice decay constant, in which all the 
explicit lattice spacing dependence cancelled from the expression up to a known shift from the 
hairpin interactions, Eq. (|6.40|) . We contrast this with the amplitude expressed in terms of the 
LO xPT parameters, Eq. (|6.32|) in which all of the lattice spacing dependence and partially 
quenched effects which one would expect, contribute to the amplitude. We then see that the 
naive guess of the form of the amplitude would add 3 extra undetermined fit parameters to 
the 1 = 2 scattering length, but that in reality, through the one-loop order there are no extra 
undetermined parameters. We now move on to discuss the mixed action effective theory which 
will allow us to understand this simple continuum like behavior of the vrvr scattering length and 
for which other processes this behavior holds, and more importantly how it breaks down. 

7.2 Mixed Action Effective Field Theory 

In this section we will show that for a mixed action theory with Ginsparg- Wilson valence quarks, 
most of the lattice spacing dependence for purely valence quantities can be parameterized as 
a multiplicative renormalization of the continuum xPT LECs, through the leading loop order. 
This is in contrast to a theory with Wilson valence quarks, in which every LEC receives additive 
as well multiplicative renormalizations. The benefits of this are obvious, and in general the 
additive corrections tend to be significant for Wilson and staggered quarks. Of course, for a 
dynamical simulation with GW valence and sea quarks, all of the lattice spacing dependence, 
aside from the 0(4) rotational breaking operators as in Eq. (|7.2|) . can be parameterized as 
multiplicative renormalizations of the continuum LECs. It is in this sense that it is useful to 
think of a mixed action theory with GW valence quarks as the full GW theory with perturbative 
breakings of symmetry between the valence and sea quarks. 

As has been shown in Refs. |17HI249] . the LO Lagrangian is given by the partially quenched 
Lagrangian, and thus is invariant under the full partially quenched symmetry. The first correc- 
tion arising from the mixed action is given in Eq. (|6.13|) and leads to an O(a^) shift to all the 
mesons composed of one valence and one sea quark. This does not depend upon the type of sea 
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quark used, as this operator does not break the taste symmetry of the staggered quarks. The 
value of Cmix will of course depend upon the type of sea quarks however. Through the leading 
loop order, the effect of this operator is to simply transform all valence-sea meson masses to 
their mixed action form 



Bo{mv + ms) Bo{'mv + ms) + a^A 



Mix 



(7.7) 



not just for the meson masses, but also for all interaction-vertices in tttt scattering and for the 
contributions to the meson decay constants. 

At the next order, things get more interesting, but as we will show, the continuum like 
behavior still holds. It is instructive to first consider what happens to the Gasser-Leutwyler 
operators. For partially quenched xPT the Gasser-Leutwyler operators retain their exact form 
with the replacement of tr — > str, and the LECs retain their exact values as in QCD. There is 
one additional operator which must be added to the mesonic chiral Lagrangian, but it has been 
shown that it does not contribute to physical observables until one order higher |132j , and so we 
shall ignore this operator in the following discussion. When considering a mixed action theory 
with Ginsparg- Wilson valence quarks, we find that at 0{mg) the standard Gasser-Leutwyler 
Lagrangian splits into more operators. For example, single super-trace operators split into 
two operators, one for the valence sector and one for the sea sector. The double super-trace 
operators split into three or four operators, depending upon whether the two super-traces are 
over the same or different sets of fields. To be precise, we give here two examples, 



2Bo L5 str (df^Ed'^T.'^ (mq^^ + Sm^ ) PQ MA 



2B^ C^^L5 str (pya^Sa^St Vv (rngSt + S 



+ 2^0^ Ci^L5 str (Vsd^^df'J:^ Vs (rngSt + S 



4BiLe 



str ( rngS^ + Sm 



PQ MA 



str (VvB^ (mgSt 



Q 



+ 8Le C^^ str (VyB^ (mgSt + Sm^)) str (VsB^ (mgSt + Sm^ 



+ 4L6 Ci^ str ( VsB^ [ mgST + Sm, 



(7. 



where in the above equation we have distinguished between the valence and sea quark conden- 
sates, Bq^ = Bq +0(0^). We have written the Lagrangian in this way because in the continuum 
limit, we know these operators have to recombine to form the Gasser-Leutwyler Lagrangian, as 
these operators are not distinct under the (PQ)QCD symmetries. In fact, we can write each of 
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the new coefficients as 




1 + 0{a^Kl) , (7.9) 



wliere i is for all relevant Gasser-Leutwyler LEG. We see also that assuming the lattice spacing 
corrections to these LECs are perturbatively small, these mixed action effects which break the 
PQ symmetry are of O(p^a^) and thus do not enter until a higher order than we considered in 
Chapter Understanding how these different operators contribute to m-„^, and %rn provides 
valuable insight into how all the lattice spacing corrections enter mixed action theories. We 
now track how the sea quark contributions from the operators in Eq. 1)7. 8() . as well as the 
equivalent operators related to the Gasser-Leutwyler operator with coefficient L4, contribute to 
these observables and then show how all the lattice spacing dependence enters. 

Computing the contributions to the above mentioned observables at O(p^) is a straight- 
forward exercise as this only involves tree level graphs. To be specific, we consider a 2-flavor 
theory and find that 

6ml = -CrL.ml ^ str(mQ) + L,ml ^ str(mQ) , 

<5/. = /CrL4-^str(mQ), 

5Z^ = -CrL,^str{mQ) , 

ST,i7^ = CFL4^str(n.Q) + ^ ^^^6^ str(rnQ) , (7.10) 

and when we add all these effects to the 1 = 2 tttt scattering length expressed in terms of the 
lattice parameters, we find 

6Ti,=^ = 0, (7.11) 

as claimed in Chapter^l Operationally, we understand this result as coming from a factorization 
of sea-quark effects, which was made explicit by the mixed action breaking of the Gasser- 
Leutwyler operators. Of course this same analysis explains why there is also no contribution 
at this order arising from partial quenching, aside from the hairpin contributions in the t- and 
ti-channel graphs of Fig. 16.21 The real power of this analysis comes when one then constructs 
the operators at this order which encode the lattice spacing effects. Regardless of the type 
of sea-quarks used, when building the mixed action effective Lagrangian with Ginsparg- Wilson 
valence quarks, one can determine that all of the operators which can give rise to lattice spacing 
effects must have those effects come from the sea quarks, and thus be of the form of the above 
mentioned operators with coefficients C^^ . In other words, the lattice spacing effects factorize 
in the same fashion as the sea quark mass effects shown above. 

One then realizes that all the lattice spacing operators one can construct which will con- 
tribute to quantities composed of Ginsparg- Wilson valence quarks, must all behave as mul- 
tiplicative renormalizations to the continuum operators, aside from the lattice spacing 
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Figure 7.1: This graph represents corrections to the pion mass from the Lagrangian, Eq. (|7.8)1 . 
For staggered sea quarks, this contribution is taste breaking. 



corrections which enter in the hairpin interactions. Not only are the lattice spacing correc- 
tions then multiplicative renormalizations to the continuum expressions, in contrast to all other 
lattice actions to date (aside from a full Ginsparg- Wilson theory), this also implies that for suf- 
ficiently small lattice spacings, these corrections are also small. There is now growing numerical 
evidence in support of these ideas given by the recent mixed action LQCD simulations which 
employ Ginsparg- Wilson valence quarks with staggered sea quarks [13 EHl 1^ IHHl IMl inHj • We 
would then like to understand how this simple picture breaks down, which we address in the 
next section. 



7.3 Taste Breaking in Mixed Action Lattice QCD 

In this section, we describe how this continuum like picture of the mixed action theories with 
Ginsparg- Wilson quarks breaks down. We will use the specific example of staggered sea quarks, 
both because currently, mixed action simulations are predominantly performed with staggered 
sea quarks, and because the oddities of the staggered theory help to highlight the points we 
are making. In fact the analysis of the previous section directly leads to an understanding of 
how the taste breaking of the staggered potential gives rise to contributions which can not be 
absorbed by multiplicative renormalizations of the continuum Gasser-Leutwyler LECs. To be 
very specific, we shall focus on the operator 

£ = SC^'^Le str (Vy (Smg + mQ^t) ) str (Vs (smQ + mQ^t) ) , (7.12) 

and work to 0{p^). We shall compute the pion mass and consider the contribution from the 
Figure 17^1 in which the pion loop comes from the second str in Eq. 1)7. 12() . This operator will 
give a contribution to the pion mass, 

2 _ 32C^'^L6 ml rit ABojmg + m^/) 

where the sum is over the various staggered taste mesons with nt counting the weighting of 
the i*'' taste meson for this loop, whose masses are given by Eq. (|6.16jl . Here we then see a 
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contribution to the meson masses which can not be treated as a multiphcative renormahzation 
of the continuum LECs. We note as a final comment that this feature wih also enter in the 
same manner in the baryon spectrum j3U8j . 

We then conclude this chapter by re-emphasizing that for quantities which are protected 
by chiral symmetry, eg. the pion mass, the scattering length, etc., the good chiral properties 
of Ginsparg- Wilson quarks, even when only used in the valence sector, protect these quantities 
from receiving chiral symmetry breaking corrections from the lattice spacing through the one 
loop order, irrespective of the sea-quark type. 
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Chapter 8 
CONCLUDING REMARKS 

Presently, lattice QCD simulations are performed with lattice spacings which are not negli- 
gibly small, lattice volumes which are not significantly larger than the typical hadron systems 
of interest and quark masses which are larger than those of nature. In order to make a rigorous 
connection between these lattice QCD simulations and nature, we need to understand how these 
lattice artifacts modify the computed correlation functions and thus the observable quantities 
of interest. The tool to do this is effective field theory. In this work, we have presented devel- 
opments and applications of effective field theory techniques to various observable quantities 
which are necessary for the lattice QCD simulations of today and tomorrow. 

In particular, in Chapter |2l we have extended the partially quenched heavy baryon La- 
grangian necessary to compute the masses of the lightest spin-^ and spin-| baryons to 0{mq)^ . 
This is necessary as the convergence of the effective field theory expansion for baryons is not 
as good as with the mesons, and so must be pushed to higher orders to test the convergence of 
the theory and hopefully provide more accurate theoretical knowledge of the baryon masses. In 
Chapter |H1 we showed how to include baryons into twisted mass chiral perturbation theory and 
then use this Lagrangian to determine the effects the twisted mass discretization technique has 
on the spectrum of the nucleons and deltas, which we were able to qualitatively compare to the 
available quenched twisted mass lattice QCD data of the baryon spectrum. We also explored 
the very interesting effects of simultaneously including both the isospin breaking of the quark 
mass parameter and the twisted mass parameter which leads to a mixing of quark flavors and 
obscures the definition of isospin. 

In Chapter m we showed how one can use background electromagnetic fields to determine 
the electromagnetic and spin polarisabilities of hadrons. This is particularly interesting because 
currently, we can only measure two of the four spin polarisabilities of the nucleon and therefore 
we are in a position to make a prediction of these polarisabilities with lattice QCD. From a 
theoretical standpoint, the polarisabilities are very interesting because of there sensitivity to 
the chiral physics. Unlike many baryon quantities, the leading contribution to the nucleon 
polarisabilities comes from the pion cloud and not a local point interaction. Therefore these 
quantities are particularly sensitive to finite volume modifications, and even for the present pion 
masses used in determinations of these quantities, m,r ~ 500 MeV, the leading finite volume 
modifications can be on the order of 10% corrections. 

In Chapters [5] and IHl we explored both the finite volume as well as the lattice spacing 
corrections to the two pion system. We first determined the exponential volume dependence 
of the 1 = 2 TTTT system and showed that these effects will become particularly important in 
the next generation of lattice simulations as the pion masses are brought significantly down 
from 300 MeV. We then determined the partial quenching and lattice spacing effects on the 
1 = 2 TTTT scattering length which arise from a mixed action simulation with Ginsparg- Wilson 
valence quarks and staggered sea quarks. We showed that if one expresses the scattering 
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length in terms of the lattice parameters measured from correlation functions, then the form 
of the answer is identical to continuum chiral perturbation theory up to a computably small 
shift arising from hairpin interactions of the flavor neutral mesons; all other possible partial 
quenching and lattice spacing corrections to the scattering length cancel in the final answer in 
terms of these lattice parameters. 

In Chapter d we showed that lattice QCD simulations which utilize mixed actions with 
Ginsparg- Wilson valence quarks exhibit continuum like behavior for quantities protected by 
chiral symmetry through the leading loop order, such that one can treat all the lattice spacing 
dependence as a multiplicative renormalization of the continuum operators, which is con- 
trast to all other types of fermion discretization methods presently used (aside from dynamical 
Ginsparg- Wilson sea quarks with Ginsparg- Wilson valence quarks). We then showed how this 
behavior breaks down at the next order, using the example of staggered sea quarks to show 
how the taste breaking of the staggered mesons disrupts this continuum like behavior of mixed 
action lattice QCD. 
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Appendix A 
TWISTED MASS APPENDICES 



A.l Absence of additional dimension five symmetry breaking operators induced 
by the mass splitting 



In this appendix, we show that the mass sphtting does not induce any symmetry breaking 
terms in the effective continuum Lagrangian at the quark level at quadratic order. The (mass) 
dimension six operators in the Symanzik Lagrangian we can drop for the same reason given in 
Ref. 169 , since they are either of too high order (cubic or higher in our expansion) or they do 
not break the symmetries further than those of lower dimensions. For dimension five operators, 
we will show that the only allowable terms by the symmetries of the lattice theory are those that 
either vanish by the equations of motion, or can be removed by suitable 0{a) redefinitions of 
the parameters in Cq, the effective Lagrangian in the continuum limit (the lowest order effective 
Lagrangian) . 

In the mass-degenerate case |169| . the only dimension five operator that appears is the Pauli 
term. Since in the limit of vanishing mass splitting (the isospin limit) the mass non-degenerate 
theory must be the same as the mass-degenerate theory, any additional operators induced by 
the mass splitting must be proportional to the mass splitting. These can only be of the form 



eli^OoiJ : Oo = To\{l} , Lq = {1 , , 75 , Ts^fc} , 
eq^Oi^ : Oi = {pFo , m Lq , /x Lq }\{^ rg , mra} , 



k = l,2, 3, 



dim[Oo] 
dim [Oi] 



0, 
1, 



(A.l) 



where the notation means "the set P excluding the set Q". The quantities Oq and 

Oi are all the possible independent structures with the correct dimension, which do not lead 
to dimension five operators vanishing by the equations of motion, or are not removable by 
redefinitions of parameters in Cq. However, none of these operators are allowed under the 
symmetries of the lattice theory. Specifically, they are forbidden by charge conjugation (C) and 
the pseudo-parity transformations that combine the ordinary parity transformation (V) with a 
parameter sign change 



or a flavor exchange or both 



Vi X 



F 1 



(A.2) 



(A.3) 



where 



V 



2,3 



Xp 
k = 



Uoix) ^ Uoixp) 
Uk{x) ^ Ulixp) 

^{x) -i-i/'(xp)7or2,3 



= (-x,t) 
1, 2, 3 



(A.4) 
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Structure 


C 






' F,e, 




T2, 75T2 


n, 75, 75T3 


75n 


T3 


Oi 


^75 X {1, Ti, T3} 
mT2, mj5T2 

ifiT2, i/i75r2 


Pn, lpT2 

"1 X {ti, 75, 75X3} 

i/i X {ri, 75 , 75 T3} 


Plhr2 





Table A.l: The structures of the dimension five operators that are non- vanishing by the 
equations of motion and non-removable by parameter redefinitions. They are classified by the 
symmetries that forbid them. 



and are the lattice link fields. Note that we have displayed the symmetries of the lattice 
theory jl6,S[ I179j in the form which applies to the effective continuum theory. 

In Table 1, we show explicitly which symmetry forbids each of the possible structures of Oq 
and Oi listed in iTOTl .^ We group the operators in columns according to the symmetry under 
which they are forbidden. 

The conclusion of the above discussion is that the mass splitting does not induce any ad- 
ditional operators that do not vanish by the equations of motion, or can not be removed by 
redefinitions of the parameters in the theory. Thus beyond Cq, the effective continuum La- 
grangian contains only the Pauli term to the order we work, exactly as in the mass-degenerate 
case. 

A. 2 Diagonalization of the delta mass matrix 

Here we diagonalize the tree level mass matrix for the delta states. We reiterate that the 
difference between first diagonalizing the tree level mass contributions, then calculating loop 
effects, versus calculating the loop contributions then diagonalizing, is of higher order than we 
work. To proceed, first we list all the independent operators to 0{e'^) that have tree level mass 
contributions, 

Oie^): (T^M+T^), (T^T^) tr(A^+) , (T^T^) tr(W^") 
0(am) : {T,M+T^) tr(W^") , (T^T^) tT{Wf) tr(A^+) 

0(a2) : (T^T^) tr(W^"') tr(W^"') , (T^T^) tr(Wi"'Wi"') , r'^J^W'Ty'' (>Vi"')^'^'r;'^'^ 

(A.5) 



^Most of what we show can be readily inferred from |179| . What is new here is the need for and the use 
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The tree level delta mass matrix at the order we work, Ma, is then given by 

(-A + C \ 

\{A-2B) + C 

-\{A + 2B) + C 













B 



VA , 



= U^|Cl4x4 + Ka jvA , 

where the vectors and va are vectors of the (QCD) delta basis states. 



(A.6) 



VA 



(A.7) 



and 



Ka 



( 


-A 




B 











\ 






B 
V3 












-1{A + 2B) 

B 
V3 



\ 



B 

A 



The entries in Ma are given by 

A = 2e, L/ + tr-^cosH 



D jW~ 2 A 3 ■ 2 / \ 

B = t2 a Aqcd sm (a;) , 



(Ai 



C = 2mq {-fM — — 4crvi/a Aq(-7£,cos(6<j) +2mq (^t^^^ + 2t^^+^ — _^CD ^^^^^^ 



+ a'KlcD {t + U)+ a'K^QCD ( 4*7+ cos^(t^) - 2t';- sin^(u;)) . 



.2a3 



W+ 



(A.9) 



Note that to the accuracy we work, oo can be either lvq or the non-perturbatively determined 
twist angle. 

Except for the operators 

which contribute to Ka, all other operators listed in ()A.5|) above have trivial flavor struc- 
ture, and so contribute to the identity part of Ma- Hence, to diagonalize Ma, we need only 
diagonalize Ka- The orthogonal matrix that accomplishes this is 



/ {2A-B+2X. 

%/3S 



>l/2 



{~2A+B+2X- 
\/3B 



-,1/2 



2Xy^(2yl-S+2X_)i/2 2Jfy^(-2A+B+2X_)i/2 








(2A+_B+2X+) 



1/2 






[-~2A-B+2X+ 



\ 



ll/2 



2X 



1/2 



V 











+ 



2x;' 

^/3B 



2X^''^(2A+B+2Jf+)-i/2 2X^/^(-2A-B+2X+)-i/2 y 

(A.IO) 
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where X± = ± AB + S^, and each column of S is a normalized eigenvector of Ma (and 
hence Xa also). The diagonal matrix one obtains after diagonalizing Ma is then 

(-A -B- 2X_ + 3C\ 



V = S-^ ■ Ma • S 



1 



■ diag 



V 



-A-B + 2X_ + 3C 
A-B- 2X+ + 3C 
A-B + 2X++3CJ 



(A.ll) 



where each entry in D is an eigenvalue of Ma- 

Now ii €q ^ 0, A ^ 0. Hence, since in our power counting A 
may expand X± in the ratio of -B/^ ~ 0(6^) <C 1 as 

X± = A 



0{e^) and B ~ 0{e^), we 



^ B B^ 
'^A^A^ 



A 



2 A 8A^ ^ ■ 



(A.12) 



from which it follows that 



S 



( 1 

4A 



V 



1 + 



V3B 
4A 
B 



4A 











1 

V3B 
AA 



1 






V3B 
AA 

- 

AA/ 
-6^ 



V = diag 



v 



-A + C ' 
— — — + c 

'¥ ~ ~ + 
A + C 



(A.13) 



up to corrections of 0{e ) for S and 0(e ) for V. 

If eg = 0, i.e. in the isospin limit, A = and X± = B. In this case, one can not find S 
and T> in the isospin limit by taking the limit A — > in ()A.13|) . since expansion in the ratio of 
B/A IS clearly not valid. Instead, one has to go back to Eq. HA.10|) and Eq. HA.11|) . which in 
the isospin limit reduce to 





V3 
2 





o\ 


2 


1 
2 














2 


1 

2 


\o 





1 

2 


V3 
2 ' 



V = diag 



(-B + C\ 

f + C 
B + C 

B 



(A.14) 



and the eigenvalues contained in T) given in Eq. 



(IA.14I) are the masses of the deltas at tree level 



in the isospin limit given in Eq. ()3.38|) and Eq. 1)3.39(1 . Note that as discussed in the text, in the 
isospin limit, we need not perform any mass matrix diagonalization at all, since we can simply 
rotate from the outset to the basis where the twist is implemented by the diagonal t^,. 
The new delta basis states are defined by 



v'a 



S ^ • TJA 



■ s. 



-"A ~ "^A ■ ' (A. 15) 

in which the delta mass matrix is diagonal to the order we work. By writing the old (unprimed) 
delta basis states in terms of the new (primed) basis states using the defining relations given 
above, i.e. 

^^a = 5'-ua, v^ = v'^-S~^, (A.16) 

the Lagrangian in the new delta basis can be obtained. Note that in the case where ^ 0, the 
new basis states contained in v'^ are "perturbatively close" to those contained in the va, i-e. 
the difference is 0{e ) as can be easily seen from Eq. (|03l) . This is of course not true if we 
are in a region where i? ~ j4, or ~ o'^Aq^j-j. 
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Appendix B 

COMPTON SCATTERING APPENDICES 



B.l Quenched chiral Lagrangian 

In this Appendix, we display the relevant pieces of the quenched chiral Lagrangian in the meson 
and baryon sectors and note particular pathologies of the quenched theory. In a quenched two 
flavor theory, we have valence (u, d) and ghost (n, d) quarks with masses contained in the 
matrix 

mq = diag(m„, m^, rriu, mj) , (B.l) 

where m- ^ = rriu^d to maintain the exact cancellation from the path-integral determinants 
arising from the valence and ghost quark sectors. The corresponding low-energy meson dynamics 
are described by the QxPT Lagrangian. At leading order, the form of this Lagrangian is the 
same as in Eq. (|1.57|) where the pseudo-Goldstone mesons are embedded non-linearly in E with 
the matrix $ now given by 

where 



M = ( ^ , M= ni I ] , x= T ) • (B.3) 



The matrix M contains the usual valence-valence mesons, while mesons in M are composed of 
ghost quarks and anti-quarks, and finally those in x of ghost-valence quark-anti-quark pairs. 
Unlike the partially quenched theory, there is no strong U(l)yi anomaly, and the flavor-singlet 
field, <I>o = str[$]/\/2 (along with its couplings rn-o and «#), must be retained in the theory. For 
the electromagnetic and spin polarisabilities in QQCD, no loop contributions from the singlet 
are needed to the order we work as flavor-neutral mesons are not present in loop diagrams at 
this order. Despite flavour neutral mesons being absent in loop graphs, the anomalous tree-level 
term couples the quenched singlet to the nucleon. Cancellations, however, lead to final results 
that are independent of mo and a$. 

For the quenched electric charge matrix of the valence and ghost quarks, we choose 

Q = diag{qu,qd,qu,qd) ■ (B.4) 

Notice the peculiarity that str Q = is unavoidable in the quenched theory. In the quenched 
theory, there are anomalous decays of fiavour neutral mesons into two photons. In terms of 
SU(2|2) QQCD quark fields, contributions to the anomaly from the valence and ghost sectors 
come weighted with squares of the quark charges, and we are thus not restricted to only the 
fiavor singlet current (as is the case for the strong U(l)yi anomaly). The relevant term of the 
anomalous quenched chiral Lagrangian is the same as has been detailed above in Sec. 14.5.11 
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In SU(2|2) HBxPT, the nucleons (those composed of three valence quarks) enter as part of 
a 20-dimensional representation described by a three index flavour-tensor, B. The quenched 
A-isobar is contained in the totahy symmetric three index flavour-tensor transforming in 
the 12-dimensional representation of SU(2|2). The leading-order Lagrangian describing these 
baryons and their interactions with Goldstone mesons is 

rg)^ = i{Bv VB) + 2a (BS^'BAf,) + 2/3 (BS''Af,B) + 27 {BS^'B) str^^ 

+ 27' str^^ + [ij"A,B) + (BA.T'')] . (B.5) 

In contrast to partially-quenched and unquenched chiral perturbation theory, there are two 
additional axial couplings 7 and 7' due to the presence of the flavour-singlet field. One should 
keep in mind that although we use the same notation for simplicity, all of the coefficients in 
the quenched Lagrangian have distinct numerical values from those of the partially-quenched 
Lagrangian. In the large Nc limit, the coefficients of the two theories are related |3U9| . 

Again the photon is minimally coupled in the above Lagrangian with fixed coefficients. At 
the next order in the expansion, the relevant terms that appear are 

+V|/^T^F,, [{BS>^Q^^T-) + {T^S'^Q^+B)] . (B.6) 

The PQxPT term with coefficient //^ is absent in the quenched theory. This only affects the 
Born terms of the Compton amplitude, which are essentially unknown because they depend on 
the quenched magnetic moment. Finally, the leading two-photon operators that give completely 
local contributions to the Compton scattering tensor appear in quenched chiral perturbation 
theory in essentially the same form as PQxPT. However, there are fewer operators per spin 
structure compared to the partially quenched case because of the super-tracelessness of the 
electric charge matrix. Our computation is unchanged since these terms do not contribute at 
the order we work. 



B.2 Finite volume functions 



The sums required in the evaluation of the polarisabilities at finite volume are {k 
n a triplet of integers) 



MM) 

MM) 
IC^{M) 
Cp{M) 



1 



k [|^|2 + M2] 

J^_i(M)-M%(M), 

lp-2{M) - 2M%_i(M) + M%(M) , 

T^_3(M) - 3M2t^„2(M) + 3M%_i(M) - M%{M) 



with 



(B.7) 

(B.8) 
(B.9) 
(B.IO) 
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At infinite volume these can be simplified using, 



1 r(/?-|) 1 



(47r)i r(/3) (M2)^-f 



(B.ll) 



for /3 > 3/2. 

In numerically evaluating these sums, it is useful to note that 



5^ (|n|2+x2)/3 5Z 



r(/3) 



+ 



7r2 



r(/3) Jo 



I"! 

e-— + 1 



(B.12) 



where En{x) is the exponential integral function. This form is valid for /3 > |, x G IR and the 
remaining sums converge exponentially fast in |n|. 
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